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ABSTRACT 

Motivated by a new wave of kinematical tracers in the outer regions of early-type galaxies (ellipticals 
and lenticulars), we re-examine the role of angular momentum in galaxies of all types. We present new 
methods for quantifying the specific angular momentum j, focusing mainly on the more challenging 
case of early-type galaxies, in or der to der ive firm empirical relations between stellar j\ and mass 
Mi, (thus extending the work bv lFaUlll983f ). We carry out detailed analyses of eight galaxies with 
kinematical data extending as far out as ten effective radii, and find that data at two effective radii 
are generally sufficient to estimate total reliably. Our results contravene suggestions that ellipticals 
could harbor large reservoirs of hidden in their outer regions owing to angular momentum transport 
in major mergers. We then carry out a comprehensive analysis of extended kinematic data from the 
literature for a sample of ~ 100 nearby bright galaxies of all types, placing them on a diagram of 
versus M*. The ellipticals and spirals form two parallel j\-M+ tracks, with log-slopes of ~ 0.6, 
which for the spirals is closely related to the Tully-Fisher relation, but for the ellipticals derives from a 
remarkable conspiracy between masses, sizes, and rotation velocities. We find that on average, the el- 
lipticals contain roughly 3-4 times less angular momentum than spirals of equal mass. We decompose 
the spirals into disks and bulges and find that these subcomponents follow similar j+-M+ trends to the 
overall ones for spirals and ellipticals. The lenticulars have an intermediate trend, and we propose that 
the morphological types of galaxies reflect disk and bulge subcomponents that follow separate, funda- 
mental j^—Mi, scaling relations. This provides a physical motivation for characterizing galaxies most 
basically with two parameters: mass and bulge-to-disk ratio. Next, in an approach complementary to 
numerical simulations, we construct idealized models of angular momentum content in a cosmological 
context, using estimates of dark matter halo spin and mass from theoretical and empirical studies. We 
find that the width of the halo spin distribution cannot account for the differences between spiral and 
elliptical but that the observations are reproduced well if these galaxies simply retained different 
fractions of their initial j complement (~ 60% and ~ 10%, respectively). We consider various physical 
mechanisms for the simultaneous evolution of j* and M* (including outflows, stripping, collapse bias, 
and merging), emphasizing that the vector sum of all such processes must produce the observed j+—M* 
relations. We suggest that a combination of early collapse and multiple mergers (major or minor) 
may account naturally for the trend for ellipticals. More generally, the observed variations in angular 
momentum represent simple but fundamental constraints for any model of galaxy formation. 



1. INTRODUCTION 

Many schemes for classifying galaxies have been pre- 
sented over the years, focusing on somewhat ephemeral 
properties such as morphology and color. Alternatively, 
one may consider three fundamental physical parame- 
ters: mass M, energy E, and angular momentum J. 
Qualitatively, these are related to the amount of ma- 
terial in a galaxy, to the linear size, and to the rotation 
velocity. 

An important advantage of these parameters is that 
they may be related back to the earlier states of galax- 
ies without having to unravel all of the messy interven- 
ing details such as baryonic dissipation, star formation, 
and morphological transformation. As an example, the 
simple assumption that J is approximately conserved 
during the collapse of gas within hierarchically-forming 
dark matter halos naturally exp lains the observed basic 
scaling relations of disk galaxies (TFall fc Efstathioull 19801 : 
IDalcanton et al.lll997HMo et al.lfl998l) . 

Here "conserved" means that the initial J* is re- 
tained at a factor of ~ 2 level, unlike E, which can 



be readily lost by factors of ~ 10 through dissipative 
collapse and radiation. Note that the "weak" conser- 
vation of total J is less restrictive and more plausible 
than the "strong" conservation of the internal distribu- 
tion of J with radius, which could be readily altered by 
secular pro cesses within disks while still pres erving to- 
tal J ("e.g.. iKormendv fc Kennicuttl 12004 see iFalU [20021 



Bland-Hawthorn 2002 for further discus- 



and lFreeman 
sion) . 

In this vein, Fall (1983], hereafter F83) introduced a 
general diagram of j* versus stellar mass M*, where 
= is the stellar specific angular momentum. 

This diagram has the important advantages that it deals 
with conservable physical quantities, and that the axes 
represent independent variables. The M+ axis embod- 
ies a mass scale, while the j* axis represents a length 
scale times a rotation-velocity scale. On the contrary, 
the standard relations between M + and circular veloc- 
ity v c (e.g.. iTullv fc Fisherl [l977t iDutton et all [20101: 
iTruiillo-Gomez et al.l 120111 ) involve correlated variables, 
since v c may be directly connected to M t . Another re- 
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Fig. 1. — The total intrinsic stellar specific angular momentum of 
galaxies p lotted against their total stellar mass, reproduced from 
IFalll Q1983T ). with corrections from a Hubble constant of h = 0.5 to 
0.7. The symbols show galaxy types according to the legend at the 
upper left; for the ellipticals (E), open circles show galaxies with 
an upper-limit estimate of j^. The dotted line shows a trend of 
2/3 

jir oc . The logarithms plotted here and used throughout the 

paper are in base 10. These j*— M* scaling relations ar e th e focus 
of this paper, and will eventually be updated in Figure [T4l 

lated parameter is the spin (A), which is useful for char- 
acterizing dark matter halo rotation, and which we will 
discuss later in this paper. 

The simple j'^-M* diagram is still charged with useful 
information for understanding galaxies, and to orient the 
remainder of our discussion, we begin by reproducing 
the original version from F83 here in Figure [TJ The only 
change is to rescale the data for a Hubble constant of ft, = 
0.7 rather than h = 0.5. These data were for late-type 
spirals (Sb and Sc) based on extended optical rotation 
curves, and for elliptical galaxies based on observations 
from their inner half-light radii, as feasible in that era. 

The first key feature to note from Figure [1] is that 
the spirals follow a fairly tig ht scaling relation of 7+ oc 
M", where a ~ 0-7 (see alsolTakase fc Kinoshitalll967t 
lHeidemannlll969t lFreemadH970t lNordsiecklll973[ ). which 
is a phenomenology that is now understood to provide a 
remarkable link between visible galaxies and their invisi- 
ble dark matter halos. F83 provided a simple theoretical 
framework in which the gaseous baryons of galaxies are 
initially mixed with the dark matter and share in the 
same j. The baryons then cool and decouple from the 
dark matter, collapsing into star-forming disks. If the 
baryonic j is approximately conserved in this process, 
both the zeropoint and the slope of the observed spiral- 
galaxy relation are reproduced. 

The formation of disk galaxies can thus be explained 
at a basic level through this long-standing picture of 
(weak) j conservation. To provide further understand- 
ing, hydrodynamical simulations of galaxy formation 
have been pursued for decades, with the j'*-M* obser- 



vational diagram from F83 as a key benchmark for the- 
ory. Attaining that benchmark has turned out to be a 
major challe nge, with early stud i es finding catastrophic 
7 loss (e.g.. iKatz fe Gunnl 119911: iNavarro fc Benzl 119911: 
Nav arro et al JU995i|N avarro fc Steinmetzl [19971 ) . 

This angular momentum "catastrophe" can be at- 
tributed partially to numerical limitations, and partially 
to uncertainties in modeling baryonic processes suc h as 
feedba ck following star formation, as reviewed by IFalll 
(2003). Over the years, the simulations have improved 
and can now come close to reproducing the j+-M+ obser- 
vations Ce.g.. iGovernato et all [20071 : lAgertz et all 120111 : 
IGuedes et alJl2011ft . although much work still remains in 
understanding both the numerics and the physics. 

Besides the angular momentum benchmark from F83 
which has become a standard ingredient in modeling the 
formation of disk galaxies, there is another aspect of the 
original j+-M+ diagram that has received relatively little 
attention: the inclusion of elliptical galaxies along with 
the spirals. The diagram thereby provides a fundamental 
diagnostic of scaling relations for all galaxies, which is 
important because there is st ill not a f ull ex planation for 
such a basic property as the iHubblel (|1926f ) sequence of 
galaxy morphologies. 

Star formation considerations aside, there is an ob- 
vious dynamical distinction between galaxy disks and 
spheroids, which are characterized by cold, ordered rota- 
tion versus random motions with fairly low net rotation, 
respectively. Differences in the conservation and distri- 
bution of j may very well be pivotal to explaining these 
differences and to governing the fates of galaxies. 

As shown in Figure [TJ F83 found that ellipticals 
followed a j'*— M* trend roughly parallel to the spi- 
rals, but lower by a fact or of ~ 6, and with mor e 
apparent scatter (see also iBertola fc Capacciolil 11975ft . 
There are several potential explanations for such a dif- 
ference between spirals and ellipticals, but the most 
plausible one is traced to a violent, clumpy gene- 
sis for spheroids. For example, mergers could nat- 
urally redistribute angular momentum from the cen- 
tral regions of a galaxy to its outer parts by dy nam- 
ical friction Ce.g.. TAarseth fc Falll [l980t iGerhardl fl98lt 
| Barnes fc Efstathioul 119871; iZurek et aJJ 119881; iBarnes 
1992tlHernauistlll99alNavarro fc Whitdll994tlHevl et al. 
[1991 IrTOnghia fc Navarro! 120071; iZavala et all 12008ft . 
Thus, j should be basically conserved but inconveniently 
locked up in unobservable components such as the dark 
halo and the faint outer stars. 

With this theoretical sketch in hand, the j+ disparity 
between spirals and ellipticals has received little further 
attention over the years. However, the scenario of angu- 
lar momentum redistribution has not yet been directly 
tested by observations — a situation that may now finally 
be remedied via the advent of new techniques for optical 
spectroscopy in galaxy halo s (with preliminary resu lts 
along these lines reported in iRomanows kvet al.l l2004). 

In this paper we re-open various questions about an- 
gular momentum in all types of bright galaxies, follow- 
ing and extending the treatment of F83. Are the j*- 
M* slopes, zeropoints, and scatter in Figure [TJ supported 
upon re-examination? Does the "missing" in ellipticals 
emerge in large-radius data? Can the j* variations be as- 
sociated with the natural dispersion in spin expected for 
standard dark matter halos, or is it necessary to invoke 
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Fig. 2. — Physically-motivated classification diagram of galaxies, 
using the parameter space of stellar mass and specific angular mo- 
mentum. The solid blue and red lines show parallel scaling relations 
for disks and bulges, which are based loosely on our observational 
results to be presented in Section \E\ Approximate positions are 
also shown for different galaxy types: Sc, Sb, Sa, SO, fE, and sE 
(the latter two being fast and slow-rotating ellipticals). 

additional baryonic j evolution? 

F83 also proposed that the Hubble sequence may be 
understood as a systematic variation in at a fixed M* 
(or equivalently, variation in M* at fixed j'*), but could 
not test this idea owing to the lack of adequate data 
for the crucial, intermediate cases of Sa and SO galaxies. 
Here we will pursue this theme, and advance a framework 
where every galaxy can be considered basically as a lin- 
ear combination of a disk and a bulge, with each of these 
components following a characteristic j t -My, scaling rela- 
tion. In this idealized model, the j'^-M* parameter space 
maps uniquely to a space of M+ and bulge fraction B/T. 

Figure [2] provides a schematic overview of this frame- 
work, showing decompositions of the Hubble sequence 
in parameter space. One of our goals in this pa- 

per will be to include observational results for Sa and SO 
galaxies in this diagram for the first time, to see if such 
systems fill in the gap (if any) between earlier and later 
types, and if bulges and disks are homologous enough to 
explain the j^—M* trends as primarily reflecting a B/T 
sequence. 

The ji,-Mi, diagram does not simply provide a basic de- 
scription of galaxies and their subcomponents, but also 
permits a novel approach to modeling the evolution of 
galaxies which is complementary to numerical simula- 
tions. As mentioned previously, there are simple models 
for the formation of disk galaxies that relate their and 
M* values to the initial conditions of their host halos. 
More generally, any stage in the evolution of a galaxy 
will involve a vector of change in the j-M diagram that 
is not arbitrary, since in real physical processes, changes 
in j and M will be linked in characteristic ways. There- 



fore the empirical offsets between the j^-M* sequences 
of different galaxy types, and of their subcomponents in- 
cluding bulges, disks, and dark matter halos, can reveal 
the evolutionary connections among them. 

We set out to explore the preceding questions and is- 
sues as follows. In Section[2]we present a methodology for 
careful estimation of in various types of galaxies and 
observations, with most of the details of its derivation 
given in Appendix lAl Section [3] uses detailed models of a 
handful of real galaxies to examine a simplified procedure 
for ji, estimation. Our updated analysis of the observed 
ji, trends in a large sample of galaxies follows, with the 
observational ingredients and their inter-correlations de- 
scribed in Section [U and the full results presented in 
Section [5] including a definitive confirmation of the large 
offset between spirals and ellipticals. These empirical j* 
trends can be considered as fundamental, enduring tools 
for constraining theories of galaxy evolution. In Section[6] 
we go on to connect the observations to generalized the- 
oretical predictions for angular momentum in a modern 
cosmological context. We summarize in Section [7] 

In addition, Appendix [X] is an important part of this 
paper, providing an extended presentation of new con- 
tent relating to the derivation of , which has been split 
off from the main text for the sake of readability. Appen- 
dices lEHD] provide data tables of and other properties 
of observed galaxies, along with detailed discussion of the 
observations and data analysis for a subsample of these 
galaxies. 

The reader looking for immediate answers to the ques- 
tions above may wish to skip ahead to the results of 
Sections I5.2I and onwards. 

2. BASIC FORMULAE: DISKS AND SPHEROIDS 

The foundation for this paper is a revised, general ob- 
servational analysis of specific stellar angular momentum 
j* for bright galaxies in the nearby universe. This quan- 
tity is most generally calculated by the following expres- 
sion: 

. = Jt_ J r rx vpd 3 r 

h ~ Ah J r p<Pr ' { ) 

where the subscript "t" denotes the "true" angular mo- 
mentum in three-dimensional space, r and v(r) are the 
position and mean-velocity vectors (with respect to the 
center of mass of the galaxy), and p(r) is the three- 
dimensional density of the population under study (gen- 
erally assumed to be stars in this project). 

For spiral galaxies, we approximate the density dis- 
tribution as an infinitely-thin, axisymmetric disk with 
an exponential surface density profile. Assuming also a 
radially-constant rotation curve, Equation ([T]) yields the 
simple expression: 

j t =2v c R d , (2) 

where v c is the intrinsic circular rotation velocity, and i?d 
is the intrinsic exponential-disk scale length. These de- 
projected quantities are relatively easy to infer from ob- 
servations because it is straightforward to estimate disk 
galaxy inclinations. Equation ^ is widely used in the 
literature (including in F83), but we will demonstrate 
explicitly that it provides an excellent approximation to 
real galaxies whose rotation curves vary with radius. 
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Fig. 3. — Fraction of enclosed cumulative quantities vs. cylindri- 
cal galactocentric radius (normalized by the effective radius R c ) for 
model galaxies with an exponential profile (n = 1 disk, top) and 
a de Vaucouleurs profile (n = 4 spheroid, bottom). A constant, 
cylindrical rotation field is assumed. The quantities are projected 
stellar mass JW* {dotted curve), angular momentum J* {dashed), 
and specific angular momentum j* {solid). The latter quantity is 
computed using the cumulative values of both J* and M* within 
the radius R. The vertical dashed line marks 1 R e . To capture 
half of j± , the observations must extend to ~ 1 Re in a disk galaxy, 
and to ~ (4-5) R c in a spheroid. 

For more general cases including elliptical galaxies Q 
there is no established recipe equivalent to Equation ([2]). 
For multiple reasons, estimating j t for these galaxies is 
much harder than for spirals. Not only are their incli- 
nations and intrinsic shapes uncertain, but large-radius 
rotation measurements are both more difficult and more 
critical. 

We illustrate the last point with some basic galaxy 
models. Adopting the simple assumption of an axisym- 
metric system with cylindrical rotation that is constant 
with respect to the intrinsic radius R, we consider both 
a disk galaxy with an exponential surfac e density profile, 
and an elliptical galaxy with a standard Ide Vaucouleurs! 
(1948) i? 1 / 4 profile. Although ellipticals are in general 
triaxial systems, the axisymmetric model is sufficiently 
accurate for our purposes. 

Figure [3] then shows the cumulative distribution of 
angular momentum (both total and specific) with ra- 
dius. For the disk galaxy, the specific angular momentum 
reaches roughly half of its total value at the effective ra- 
dius R D that encloses half of the stellar light. This implies 
that observational estimates of jt will be relatively easy 

1 We use the term "spheroid" to mean a pressure-dominated 
stellar system (which may also rotate). A "bulge" is the spheroidal 
component of a spiral galaxy. An "elliptical" is a galaxy with only a 
spheroidal component, although many galaxies commonly classified 
as ellipticals probably have embedded disklike components, similar 
to those in lcnticulars but less obvious. We consider jointly the 
ellipticals and lenticulars under the general rubric of "early-type" 
galaxies. 



for disk galaxies. 

For the elliptical galaxy on the other hand, the halfway 
mark for j t is reached at 4.5 R e . This is because ellipti- 
cals contain a fairly large fraction of their light in their 
outer regions where the radius lever-arm in r x v is large. 
The implication is that observations of elliptical galaxies 
need to extend to much larger radii than for spirals, in 
order to be confident of capturing the total jt ■ 

Typical stellar kinematics observations in 1983 ex- 
tended to ~ 1 R c , and even today, only a small handful of 
galaxies have been observed kinematically out to ~ 5 R e , 
which means the positions of the ellipticals in the orig- 
inal j*-M* diagram (Figure [T]) were highly uncertain, 
and continue to be challenging to determine with surety. 
Fortunately, after a great deal of experimentation, which 
we will discuss below, we find that there is a heuristic 
approach where observations around ~ 2 R c can be used 
to estimate the total jt of ellipticals with reasonable ac- 
curacy. 

Returning to a general framework for estimating jt 
from observations, there is not only the challenge of ex- 
tending the data to large radii, but also of having only 
three of the six phase-space quantities in Equation ((TJ) 
accessible (i.e., the projected positions and line-of-sight 
velocity). Even the projection of jt on the sky involves 
unobservable velocity components tangential to the line 
of sight, and requires additional modeling assumptions. 

To cope with these issues, we will model the observed 
rotation and luminosity profiles of galaxies and convert 
these to jt estimates using approximate deprojection fac- 
tors. Although these factors are based on highly simpli- 
fied models, the dominant source of uncertainty is still 
the limited extent of the data to large radii. 

We derive in Appendix [X] two alternative expressions 
for estimating jt from observations, both of them based 
again on the simplifying assumption of cylindrical rota- 
tion. The first expression starts with a detailed calcula- 
tion of a "projected" specific angular momentum proxy 
that can be estimated directly from observations: 

^ p jT,(x)xdx 

Here ^rot,p(a^) is the observed profile of rotation velocity 
along the projected semi-major axis x, and S(x) is the 
surface density profile, again along the semi-major axis. 

The quantity j p is related to jt through a "deprojec- 
tion" factor Cc 

3t = C t j P . (4) 

Therefore the problem of estimating jt separates into two 
parts: the calculation of j p from observations, and the 
factor Cj which can be calibrated from theoretical mod- 
els. 

As we describe in Appendix [A] this latter factor has 
some dependence on the detailed density-velocity struc- 
ture of the galaxy, but is primarily a function of the 
inclination i relative to the line of sight. For thin-disk 
galaxies, it is simply Q = (sini) -1 . With spheroidal 
galaxies, there is an additional dilution effect that comes 
from the line-of-sight intersecting the rotation field at 
non-tangent points. In principle, this effect is dependent 
on the detailed shape of the rotation profile, but we have 
found with simplified test models that such variations 
can be neglected in practice. We also find that as long 
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as the major-axis radius x, rather than a circularized ra- 
dius R, is used in Equation then C, is insensitive to 
galaxy flattening. 

A general approximation to Ci as a function of incli- 
nation is provided by Equation (|A29|) . It is normally 
difficult to determine i for spheroidal galaxies, and we 
will when needed adopt inclination-averaged values. 

Equation ([3]) yields accurate results that are commen- 
surate with the quality of modern observations, but in- 
volves numerical integration, and careful compilation of 
E(x) and v TO t lP (x) profiles along with extrapolation be- 
yond the bounds of the data. 

We could in principle simplify the problem further by 
using parametric models for v I0 t, p (x) and £(x). Unfortu- 
nately, the diversity of observed rotation profiles (when 
non-spiral galaxies are considered) defies parametriza- 
tion. We can at leas t adopt for the surface density the 
general ISersid (|1968l ) law which accurately represents a 
wide range of galaxy types: 



£(x) cx exp —b n (x/a e ) i '' r ' 



(5) 



where a e is the effective radius along the semi-major axis, 
and the shape index n determines the steepness of the 
outer density profile (higher values are shallower: e.g., an 
exponential disk profile has n — 1 and the de Vaucouleurs 
law for ellipticals has n = 4), while b n is a numerical 
function of n [Equation (IA27[) ] . 

We use this E(x) simplification in practice when deriv- 
ing j p from a detailed v TO t tP (x) profile in expression ([3]). 
We also generally base our £(x) profiles on observations 
of stellar surface brightness profiles I(x), assuming for 
simplicity that there are no variations of stellar mass-to- 
light ratio with radius (e.g., due to dust). 

Our second method is a quick-and-dirty shortcut for 
estimating j t , as needed to generate an initial overview 
of the trends for a large sample of galaxies. We simply 
calculate the following linear scalar expression [derived 
in Appendix [X] from Equation ([5])]: 



jp k n 



(6) 



where j p means an approximation for j p , v s is the ob- 
served rotation velocity at some arbitrary measurement 
location x s , and k n ~ 1-5 is a numerical coefficient that 
depends on the Sersic index n of the galaxy [see Equa- 
tion flXSTJ) ]. As in Equation (J4j) , j p is multiplied by C* 
to provide an approximate jt- Here the basic idea is that 
a galaxy can be represented by a characteristic observed 
rotation velocity scale v s , a length scale a e , and a fac- 
tor k n that relates to the moment of inertia (discussed 
further below). 

The heuristic approximation that we make here is to 
select v s at x s ~ 2 a c for all galaxies. We will show in 
the next section that this choice allows us to estimate j p 
with an accuracy of ~ ± 0.1 dex, which is good enough 
to start making some interesting inferences about trends 
in it- 

For n — 4 spheroids, the expression equivalent to Equa- 
tion for spirals is: 

j t =3.03v s R c , (7) 

for a median, unknown inclination (Equation (|A32[0 . An 
important concept with the more general expression ([6]) 



is that k n increases strongly with n; for fixed galaxy size 
and rotation velocity, a more extended luminosity profile 
implies a higher j p owing to the large fraction of mass 
residing at large radii. This also means that a spheroidal 
(n ~ 4) galaxy with the same observed rotation v s and 
size a e as a spiral has a larger specific angular momen- 
tum. Late-type and early-type galaxies near the L* char- 
acteristic luminosi ty do have simila r sizes for the same 
stellar mass (e.g.. IShen et~alll200l . Therefore we can 
already make the basic prediction that if j p at a fixed 
mass is independent of morphology, then the early-types 
should have v s values relative to late-types of ~ fci/fei, 
i.e., lower by a factor of ~ 2. 

The formalism that we have outlined here represents 
a modest extension of the simpler methods in F83. The 
improvements introduced here include allowance for a 
range of luminosity profiles (not only n — 1 and n = A), 
and better treatment of elliptical galaxies where rota- 
tion at large radii is critically important. It also be- 
comes more straightforward to understand the interplay 
between observations and uncertainties in the esti- 
mates, as explored in the next section. 

3. OBSERVATIONS: ANALYSIS METHODS 

Before we move on to j'^-M* analyses of a large sample 
of galaxies, we examine a small sample in more detail. 
The goals here are to illustrate the nature of the available 
data, to demonstrate that the simplified Equations ([2]) 
and ^ are good approximations to a full treatment with 
Equation , and to understand some systematic effects 
in the and M* determinations. 

Because this paper is concerned with the angular mo- 
mentum bound up in the stellar components of galaxies, 
the preferred kinematic tracer comes from integrated- 
light absorption-line spectroscopy. In many cases, such 
data do not extend to large enough radii, so we make use 
of additional tracers as proxies for the field stars: cold 
and warm gas, planetary nebulae (PNe), and metal-rich 
globular clusters (GCs). 

We consider disk- and bulge-dominated galaxies in Sec- 
tions 13.11 and 13.21 respectively. We evaluate our simpli- 
fied jp estimate © in Section 13.31 describe our mass 
estimates in Section 13. A\ and then consider systematic 
uncertainties in Section \~ 



3.1. Disk- dominated galaxies 

The most straightforward galaxies for estimating angu- 
lar momentum are the gas-rich spirals, since the stel- 
lar rotation profile, which cannot always be measured 
directly, follows the gas rotation profile to a good ap- 
proximation. Also, the observed rotation can easily be 
corrected for projection effe cts in order to recover the in- 
trinsic value (see Appendix IA.2j) . The detailed analysis 
below is overkill for these galaxies, whose jt can be read- 
ily estimated through Equation |(3J), but we wish to illus- 
trate how our more general treatment works for them, 
before moving on to the spheroids. 

We consider two real galaxies: NGC 3054 and 
NGC 3200, which are well-studied disk-dominated spi- 
rals from th e class ic optical rotation curve analyses of 
iRubin et all (fl982T ). These cases are chosen to bracket 
the typical range of inner rotation profile shapes for spi- 
rals (slowly and rapidly rising, respectively). 
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Fig. 4. — Observed rotation- velocity profiles of two spiral galaxies 
(NGC 3054 and NGC 3200) vs. semi-major axis radius (renormal- 
ized by the effective radius). Each galaxy is labeled with its Hub- 
ble ty pe. The data are ionized gas velocities from IPizzella et al.l 
1120040 . The solid curves with shaded regions show power-law fits 
(with uncertainties) used to extrapolate the rotation velocity to 
larger radii. See main text and Appendices \K\ and [B] for further 
details. Dotted horizontal lines show the characteristic rotation 
velocity v s for each galaxy; the approximate intersection with the 
corresponding rotation-velocity profil e is m arked with a $ symbol 
and defines the radius x 3 (see Section 13.31 . 

We take the long-slit major- a xis ionized-gas kinematics 
data from Piz zella et al.l (|2004f) , shown in Figure 0] after 
a modest amount of re-binning. These rotation profiles 
have high-frequency bumps and wiggles that are presum- 
ably caused by local perturbations such as spiral arms. 
Fortunately, these features tend to average out when cal- 
culating a cumulative j and are not important in this 
context. 

To calculate the projected specific angular momentum 
jp, we carry out a piecewise integration of Equation ([3]), 
using the major-axis rotation- velocity data v TO t :P {x) up 
to ~ 2 a c , along with simple power-law extrapolations at 
larger radii, as shown in Figure |4] For we use an 

exponential model [n = 1 in Equation J5J], with the disk 
scale-lengths i?d taken from r-band photometry as we 
will discuss in the next section. Note that a c — 1.68i?d 
for a pure exponential disk. 

The resulting cumulative j p {< x) profiles with radius 
for these galaxies are shown in Figure EJ Here it would be 
trivial to convert j p (< x) immediately to jt(< R) using 
the known inclinations of these galaxies, but our general 
strategy is to focus first on the direct modeling of the 
observations for all galaxies, and later apply the depro- 
jection factors C», which involve different systematics. 

It can be seen that j p hardly changes outside ~ 3 fl c , 
and that the large-radius extrapolations make very little 
difference: the regions outside ~ 2-2.5 a c (~ 3-4 i?d) 
contain only ~ 8%-15% of the total luminosity, and con- 
tribute only ~ 15%-25% of the total j p (half of j p is 



Fig. 5. — The cumulative projected specific angular momentum, 
j P {< x), of several nearby galaxies as a function of semi-major 
axis radius (with log axes), based on modeling of kinematic ob- 
servations. Solid curves show the best-fit models, with shaded 
regions illustrating the uncertainties (including those due to ex- 
trapolations at large radii). See Table H3lfor the distances and a e 
values adopted. For most of the galaxies, j p has nearly reached its 
asymptotic value by x ~ 5 a c . 

enclosed within ~ 1.2 a c ~ 2i?d! Figure Given rea- 
sonable extrapolations of the data, the total j p for these 
two galaxies, using our basic modeling assumptions, is 
constrained to ~ 5% (~ 0.02 dex). 

Thus the kinematics is not a major source of uncer- 
tainty for jt estimation in disk-dominated galaxies. Ad- 
ditional complications that we have not considered here 
are deviations of the disk surface density profile from a 
simple constant mass-to-light ratio exponential model, 
and inclusion of a bulge (to be discussed later). We will 
examine more general systematic uncertainties in Sec- 
tion r~ 



3.2. Bulge-dominated galaxies 

We now turn to the novel component of this paper, 
which is the careful treatment of j t in early-type, bulge- 
dominated galaxies. Figure [3] demonstrated that tradi- 
tional observations within 1 a c provide little assurance 
about the total angular momentum content of these sys- 
tems, while even current cutting-edge observations out 
to ~ 5 a c might in principle not be adequate. 

Here we analyze a sample of eight real galaxies in de- 
tail in order to characterize the accuracy of jt estima- 
tions. Seven of these galaxies were chosen because of the 
availability of high-quality extended kinematic data us- 
ing inte grated stellar light spectroscopy from two recent 
papers (|Coccato et~al] l2009t fProctor et all 12009( 1. Both 
papers represent the first installments of systematic sur- 
veys of early-type galaxies in the local universe, and there 
is no obvious selection bias for the seven galaxies. Five 
of them are "ordinary" near-L* early-types with central 
"fast-rotator" kinematics as is typical for such galaxies 
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Fig. 6. — Rotation-velocity profiles for eight early-type galaxies. See Figure[4]for further details, including an explanation of the shaded 
uncertainty regions. For ease of inter-comparisons, the vertical axis of each panel has been scaled according to the velocity dispersion of 
the galaxy at at 2 a e , which is marked in each panel by a * symbol. Note the dashed lines at zero rotation velocity in some cases. The 
galaxies show a diversity of rotation-velocity trends with radius. 

full stellar populat i on in the regions of overlap (e.g. , 
Coccato et al.l^OOl iMcNeil et aljboiOHDas et al.ll2011t 



(|Kormendv fc Bender| [l996: Em sellem et"aTll2011ft . The 
other two (NGC 1407 and NGC 4374 = M84) are exam- 
ples of round, bright "slow rotators" that are comm on in 
high-density environments (jCappellari et al.ll2011bf ). 

Five of these galaxies also have PN or GC kinematics 
data available (jCoccato et al.l 120091 : iRomanowskv et all 
2009) , which we incorporate into our analysis in order to 
extend the range of galactocentric radii probed. We in- 
clude an eighth galaxy in our sample, NGC 5128 (Cen A), 
because it has the most extended (PN) kinematics data o f 
any early- type galaxy in the literature (jPeng et al.ll2004ft . 
It may also be the re mnant of a recent major merger (e.g., 
iBekki fc Penal 12006ft . which as discussed in Section Q] is 
expected to generally transfer angular momentum into 
the outer regions. Analysis of this galaxy thus provides 
a golden opportunity to search for the "missing" angular 
momentum, and to see if any clear jt difference emerges 
with respect to the other galaxies in the sample. 

The use of PNe and GCs to provide proxies for 
stellar kinematics may seem risky, given the consider- 
able uncertainties that remain about the parent stel- 
lar populations of these tracers. However, in most 
galaxies studied to date, both the density and kine- 
matical profiles of PN and metal-rich GC systems 
have been found to correspond well to those of the 



McNeil-Movlan et all l2012t iCortesil 120121 : iPota et al l 
2012ft . We have also verified that this is generally the 
case for the galaxies in our sample. 

Further details of the observations as well as of the 
kinematical modeling are provided in Appendix [B] along 
with the resulting rotation and angular momentum pro- 
files. It should be emphasized that the careful, homoge- 
neous construction of these profiles is laborious, which is 
why the current sample of galaxies that we consider in 
detail is relatively small. 

The rotation- velocity profiles of these eight galaxies are 
summarized in Figure [6] Unlike the spirals (Figure [4]) , 
the early-types show great diversity in the characteris- 
tic shapes of their profiles. Some are fairly constant 
with radius, others plummet rapidly from a central high 
value, and one continues increasing to the limits of the 
data. This diversity is not simply a matter of inclination, 
as can be seen by the divergent cases of NGC 821 and 
NGC 2768, which are both highly flattened and proba- 
bly close to edge-on. We thus find that the central rota- 
tion properties of early-type galaxies cannot be used to 
reliably estimate the total angular momentum content, 
and there is probably no simple function that universally 
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characterizes their full rotation-velocity profiles. 

As with the spirals, we fit power laws to the outer 
regions of the rotation data in order to extrapolate to 
larger radii (see Appendix [B] for further details). We 
then use Equation ([3J to calculate profiles of cumulative 
jp with radius, which we plot in Figure [5] Even though 
the data do not reach the total asymptotic value for j p , 
the requirement of a smooth power-law extrapolation for 
the rotation-velocity profile does in most cases strongly 
limit the total j p , which is typically determined at the 
± 15% level (± 0.06 dex). The radius enclosing half of 
the total jp varies from galaxy to galaxy depending on 
the shape of its rotation-velocity profile: 0.7-3 a e (for 
the two spirals, it is 1 a e ). 

The exceptions to these findings are the two bright, 
round ellipticals NGC 1407 and NGC 4374. Figure G2 
shows that much of the angular momentum in these 
galaxies is found at very large radii (half of j p within 
9 o e and 4 a c , respectively), as expected from their fairly 
high Sersic indices of n ~ 4-8 (the ordinary early-types 
have n ~ 2-4). However, beyond the usual uncertainties 
introduced by extrapolating the rotation velocity, there 
are a couple of other practical considerations. 

One issue is that although these particular galaxies 
have relatively well studied surface brightness profiles, 
many such massive ellipticals do not, with their n and a 
values poorly known. This situation produces "double 
jeopardy" for angular momentum estimation, since both 
the luminosity and the rotation-velocity profiles at very 
large radii are important yet poorly constrained. 

The other issue demonstrated with NGC 4374 is that 
its cumulative j p has not yet converged at the (estimated 
total) virial radius of ~ 35 a c , so it is not clear how its 
angular momentum should even be defined. This class of 
high-n galaxies is clearly problematic, and we will con- 
sider any jt results on them to be tentative for now. 

Figure [S] also reveals a first glimpse of the basic result 
of this paper. For most of the early-types in the sam- 
ple, there is relatively little angular momentum hidden 
beyond ~ 1-2 a c , and their total values of j p are lower 
than those of the spirals. We will make more detailed 
comparisons later in this paper. 

3.3. Simple J/M approximations 

We now arrive at a question that is critical for the 
wider survey of angular momentum in the rest of this 
paper: how accurate is the simplified Equation ([6])? As 
a reminder, this j p -estimator would replace the detailed 
calculations based on Equation ([3]) that we have carried 
out in the preceding subsections, but which are time- 
consuming to carry out for a larger sample of galaxies, 
and are not even possible for cases without very extended 
kinematic data. 

In Appendix [2] we have motivated the construction of 
Equation © via toy models of galaxies, and calculated 
the corresponding coefficient k n . We will now apply this 
formula to the set of ten real galaxies just discussed (both 
late- and early- types) , and find an optimum radial loca- 
tion x s for measuring the characteristic rotation velocity 
v s . 

For each galaxy, it is straightforward to find the con- 
stant value of v s which when substituted in Equation ([3J 
yields the same j p as with the full observed rotation- 




log j [km s -1 kpc] 

Fig. 7. — Comparison of a simple projected specific angular mo- 
mentum estimate [j p ; Equation (O] with the more accurate value 
(jp). Results are shown for ten different galaxies, each with a 
choice of three reference radii: x s /a e = 1 (red crosses), 2 (green 
filled circles), and 3 (purple open circles). Some of the points are 
given a 0.02 dex horizontal offset for visibility. The dashed and 
dotted lines mark the one-to-one relation with a ± 0.1 dex scatter. 
The optimal choice here for x s is 2 a c . 

velocity profile. These results are listed in Table E[3]aiid 
shown in Figures [4] and [6j where the intersection of v s 
with the rotation-velocity profile determines the char- 
acteristic measurement radius x s . As an example, for 
NGC 821, it is clear that x s ~ 2a c . For NGC 4494 on 
the other hand, a broad range of choices for x s would 
work, owing to its nearly constant rotation- velocity pro- 
file. 

Considering this issue in more detail, we calculate 
j p using Equation (|6]) with an arbitrary choice for x s 
(which in turn determines a guess for v s from the ob- 
served rotation velocity at this radius). The results for 
x s /a c — (1, 2, 3) are shown in Figure[7J plotted against j p 
calculated in full from Equation ([3]). It can be seen that 
x s /a c — 2 provides a reasonably good match between j p 
and j p for all of the galaxies in this sample. The other 
radius choices fare worse, owing to galaxies like NGC 821 
that have rotation- velocity profiles with a distinct tran- 
sition between the inner and outer regions near 2 a e , and 
thus v s measurement elsewhere would be biased. 

Now to home in more finely on a choice for x s , in Fig- 
ure [5] we present the ratio of estimated and "correct" j p , 
as a function of the chosen x s , for each galaxy. Some 
of the galaxies permit a broad range of choices for x Sl 
while others do not. Especially noteworthy again are the 
galaxies like NGC 821 and NGC 3377 which have sharp 
drops in their rotation-velocity profiles, so v s measured 
at small radii would overestimate j p by factors of ~ 2-3. 

We do not find a strong correlation between n and 
optimal X s clS expected from the simple models we con- 
structed in Appendix lA.4[ the dominant effect on x s with 
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Fig. 8. — The logarithmic ratio between simple estimates of pro- 
jected specific angular momentum [Equation (0] and more accu- 
rate values [Equation (3)]) vs. the rotation-measurement radius x s 
in units of the effective radius. Each point indicates a sample ratio 
for an individual galaxy, with error bars indicating the kinematics- 
driven uncertainties in total j p from the detailed models. Results 
are plotted for 10 galaxies: two spirals (orange profiles with open 
circles), six ordinary early- types (blue profiles with filled circles), 
and two giant ellipticals (red profiles with filled squares). As in Fig- 
ure [7] x s ~ 2a c provides a good measurement location, resulting 
in minimal scatter and bias for the angular momentum estimates. 



the real galaxy sample is the scatter in the shapes of the 
rotation-velocity profiles. Future detailed analyses of a 
larger sample of galaxies may reveal systematic trends 
with n that motivate improved j p estimation methods, 
but for now we stick with our simple j p approach. 

Because the real galaxies so far do not show strongly 
rising outer rotation- velocity profiles, and if anything the 
reverse, x s ~ 2 a c appears to be a good overall choice 
for the rotation- velocity measurement radius. This min- 
imizes the galaxy-to-galaxy scatter in the j p approxima- 
tion (~ ± 0.1 dex), and appears to produce little system- 
atic bias (<~ 0.1 dex). Such "errors" are comparable to 
the uncertainties from carrying out the full j p calcula- 
tions, and are therefore acceptable for our purposes in 
this paper. 

One caveat here is that this sample of galaxies is still 
small, and we cannot yet be sure of the universal valid- 
ity of our approximation, e.g., for the larger sample of 
galaxies that we will study in the remainder of this pa- 
per. However, we will show that there is no apparent 
systematic bias, i.e., the overall scientific conclusions are 
consistent with the subset of detailed j p profiles. 

3.4. Stellar mass estimates 

So far we have focused on estimating j*, but the other 
key component in constructing the j^—M* diagram is of 
course the stellar mass M*. Assuming that we have a 
well-determined surface brightness profile I(x) or total 



luminosity, we then need to know the stellar mass-to- 
light ratio T*. In this paper, we assume for simplicity 
that T* is constant throughout each galaxy, which also 
means that its value is not relevant in our calculations. 

Estimating in galaxies is a classic and not fully 
resolved problem. One standard approach is to use theo- 
retical models for stellar populations in combination with 
observations of the stellar light (e.g., broad-band colors, 
or spectroscopic line indices). Although there are well- 
known degeneracies between the ages and metallicities 
inferred for the stars, fort unately can be estimated 
with more certainty (e.g., iTortora et al.l l2009h . modulo 
the initial mass function (IMF) of the stellar populations. 

The IMF affects the overall normalization of T* via 
the mass contributions of late-type dwarf stars or com- 
pact stellar remnants, which are observationally difficult 
to tally. If all galaxies have the same IMF, then our 
analyses of the relative differences between galaxies in 
the j*-M* plane will be secure. There are also recent, 
indirect clai ms for poss i ble galaxy-to-galaxy IMF vari- 
ations (e.g.. lDavell2008t iTreu et al.l 120101: ITortora et all 
2010t Ivan Dokkum fc Conrovl I20llf IDutton et all 120121: 
Ferreras et al.ll2012t iSmith et al.l 120121 ). However, even 
in this case we do not expect a major impact on our 
conclusions. 

As a n example, the recent analysis of iCappellari et all 
(2012) implies that strong IMF variations tend to occur 
in only the most massive, and relatively rare, early- type 
galaxies, which would have log (M*/Mq) > 11.3 in our 
plots (based on a standard IMF). Such galaxies might 
have masses larger than our estimates by factors of ~ 2, 
but given the relatively small numbers of such galaxies 
and the weak constraints on their values, they will 
have little effect on our estimated j*-M+ trends. 

It is outside the scope of this paper to estimate T* 
for each galaxy in detail. Instead, we adopt the sim- 
plification that all galaxies have Y* t K = 1.0. The 
near infrared (NIR) A"-band is only mildly affected by 
internal and foreground extinction, is thought to be 
fairly insensitive to variations in stellar populations, and 
has uniform photometry av ailable from the 2MASS sur- 
vey (|Skrutskie et al.l 120061 ). The systematic variation 
in T^ t K across our entire sample of late and early-type 
galaxies is conventionally expected to be no more than 
-0.1 dex, based on B - V colors (jBell et al.ll2003h . al- 
though there are re cent suggestions of v ariations at the 
level of ~ 0.4 dex (jZibetti et al.l l2009t iBershadv et all 
120111 ). Our adopted value of T*^ = 1.0 is motivated 
by d ynamical estim a tes in both spirals and lenticulars 
from I Williams et all (|2009l figure 9 ) , and corres ponds to 
an IM F midway between iKroupal (|2001f ) and iSalpeterl 



(1955) 

Our calculations of M* also require estimates of to- 
tal luminosity, Lx- However, we do not simply adopt 
the total magnitudes provided by the 2MASS archive. 
These values are not reliable f or ea rl y-type galaxies 
(e.g.. iNoordermeer fc~V erheiicn 120071; Devereux et all 
2009; Will iams et al.ll2009HSchombertll2011[ ). particularly 
the variety with extended high-n envelopes, where the 
2MASS values could be too faint by as much as 1 mag. 

Instead, we construct our own "aperture corrections" . 
We adopt the 2MASS magnitudes within the 20 th -mag 
isophote, A20, and use the best available optical pho- 
tometry for each galaxy along with a Sersic model fit to 
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estimate the fraction of the galaxy light residing beyond 

K 20 - 

This procedure neglects any bandpass-dependence in 
the light profiles I(x), which ar e often more radially 
extended in bluer bands (e.g.. Ide Vaucouleursl 119611 : 
iPeletier et all [l990t iRoediger et al.| 1201 1|) Such differ- 
ences imply Y* variations with radius ([Tortora et al.l 
l20ll . which is a reminder of the limitations of our 
constant- Y* approximation. Given our reliance on op- 
tical profiles I(x) to derive S(x) and estimate j p , as 
in Equation ([3]), for consistency we do need to use the 
optical data to extrapolate the -fT-band photometry in 
estimating M*. However, the scale- lengths a c of the 
stellar mass distributions are probably smaller on av- 
erage than the a e values that we use based on opti- 
cal luminosity distributions, leading us to overestimate 
both j p and M*. Improvement on this point could be 
made in the future by analysis of deep I(x) data at 
NIR wavelengths. NIR spectroscopy would then also 
be needed fo r full consistency of both j u and M* es- 
timates (e.g.. iSilge fc Gebhardtl [200l iSilva et al.l [20081: 
I Vanderbeke et al.ll2011lT ~ 

3.5. The j+-M* diagram 

Here we focus on the j*-M+ plane, our ultimate desti- 
nation in this paper, but for now considering the pro- 
jected specific angular momentum j p rather than the 
true jt in order to isolate various effects that are disjoint 
from inclination uncertainties. Figure [3] shows our de- 
tailed galaxy sample where cumulative j p (< R) — J P {< 
i?)/M*(< R) is plotted not as a function of radius (as in 
Figure [5]) but of enclosed projected stellar mass, M*. 

For reference, we show dashed lines corresponding to 
jp cx M", with a — 2/3. This value for a is motivated 
by previous observations (Section [1]), and by theoretical 
predictions for jt-M+, given constant values of an initial 
halo spin parameter A, as we will see in Section [6.11 We 
are most concerned with the locations of galaxies rela- 
tive to these tracks, and with any systematic effects that 
could shift the data in a direction perpendicular to them. 

The shaded regions of the curves in Figure [9] indi- 
cate the uncertainties due to the kinematic data, includ- 
ing the extrapolations to large radii. For most of the 
galaxies, the asymptotic position in the j p -Af* diagram 
is relatively well determined. The main exceptions are 
NGC 1407 and NGC 4374, which as discussed before are 
extended giant ellipticals whose total j p is very difficult 
to determine. The early-type galaxy NGC 2768 is also 
a concern even though the formal j p uncertainties are 
small, since there are large contributions to the total j p 
estimate from the region of extrapolation. 

An offset in total j p between the late-types and most 
of the early- types as in Figure [5] is also apparent in Fig- 
ure [9] However, the mass dimension brings the relative 
positions into sharper focus. For example, NGC 4374 
and NGC 5128 have similar j p values to NGC 3054, but 
also have larger stellar masses, which means that their 
inferred halo spins will be lower (considering distances 
perpendicular to the dashed tracks). 

We next consider some systematic uncertainties that 
apply even if the rotation- velocity profiles are perfectly 
measured. First, there is a typical distance uncertainty 
of ~ 10%. This affects j p linearly and M+ quadratically, 
moving the position of the data by a very small amount 
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Fig. 9. — The cumulative projected specific angular momentum 
of nearby galaxies (as in Figure^, now plotted vs. cumulative pro- 
jected stellar mass. The curves are solid where constrained by the 
data, and dotted for extrapolations. Circles show intervals of 1 a e , 
up to 4 a c . Error bars at the end of the NGC 3054 curve illus- 
trate the effects of systematic uncertainties (see text for details): 
diagonal for the distance, vertical for scale-length, and horizontal 

2/3 

for T*. Diagonal dashed lines show tracks of j p <x MJ , which 
represent constant halo spin. 

nearly parallel to the A tracks (see sample error bars 
marked for NGC 3054 in the Figure). 

Next we consider an uncertainty of ~ 30% (~ 0.11 dex) 
in the scale lengths a e , which translates into a similar 
uncertainty in j p [see Equation J2J|]0 Also, in some cases 
the surface brightness profile is well constrained and the 
associated j p uncertainty is very small (e.g., ~ 5% or 
~ 0.02 dex in the case of the n ~ 3 elliptical NGC 4494). 

Finally there is T*, which may be uncertain by a fac- 
tor of ~ 50% (~ 0.2 dex) and would affect M* by the 
same amount. For spiral galaxies, this is probably the 
limiting factor for inferring their A values. For the early- 
types, the inclination is generally unknown and may be 
a significant source of uncertainty for estimating j t , even 
when j p is well constrained. We will return to this theme 
in Section [5TT1 

4. OBSERVATIONS: SCALING RELATIONS AND 
DERIVATIONS OF J/M FOR THE FULL SAMPLE 

Having carried out detailed analyses of for a hand- 
ful of galaxies in the previous section, we now derive j'* 
for a much larger galaxy sample, using simpler methods. 
Besides these derivations, in this section we also examine 
some basic scaling relations for galaxies, in order to un- 
derstand the observational underpinnings of the j+—M+ 
results in the next section, and to verify that our results 
are consistent with some well-known properties of galax- 
ies. We also introduce a novel, generalized version of the 

2 In practice, the a e uncertainty is correlated with an uncer- 
tainty in the galaxy luminosity and thus in M*, but this is a rela- 
tively weak effect. 
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Tully-Fisher relation for galaxies of all types. Those who 
are keen to get straight to the angular momentum results 
may wish to skip to Section [5.21 

In order to populate the observational j^—M* diagram, 
we will use the j p approximation of Equation ([5]) which 
we have found to be generally accurate at the ~ 0.1 dex 
(~ 25%) level. The basic parameters that we then need 
for all of the galaxies are: the total stellar mass (M*) 
and its scale- length (i?d or a c ), the Sersic index n, and 
the characteristic rotation velocity v s . 

The distances to the galaxies are estimated from red- 
shifts and surface brightness fluctuations. As discussed 
in Section 13.41 M* is derived from aperture-corrected 
2MASS magnitudes rriK , assuming T*,if = 1.0. 

The other parameters are derived differently for the 
late-type and early-type samples, as we will discuss in 
Sections 14. II and 14.21 respectively. Section |4~31 brings the 
data together in an examination of basic scaling rela- 
tions, before proceeding to the final j+-M+ analyses of 
Section [5] 

4.1. Late-types 

Because spiral galaxies are dominated by their disk 
components, whose photometric and kinematic proper- 
ties are relatively straightforward to measure, past stud- 
ies of their angular momenta have generally treated them 
as pure disks, e.g., using Equation ([2]) to calculate j t . 
However, this approximation may be inadequate for the 
spirals with relatively large bulges (Sa and some Sb), 
and it is one of the goals of this paper to consider these 
components. 

With Equation © in mind, we could use values for the 
parameters n, a c , and v s that characterize the compos- 
ite bulge-disk systems (e.g., with an overall n somewhat 
larger than 1). However, the required stellar photometry 
and kinematic data are not available for a large sample of 
galaxies. Instead, we analyze disk and bulge components 
separately, make some simple assumptions for the bulges 
to compensate for the missing data, and then combine 
the disks and bulges into global j* analyses. 

We focus on the classic s piral galaxy data set assem- 
bled by Kent! (|1986l 119871 119881 ). comprising 64 galax- 
ies from type Sa to Sm, at distances ranging from 1 to 
100 Mpc. These data include r-band CCD photometry 
along with bulge/disk decompositions, and inclination- 
corrected gas-disk rotation curves from both o ptical 
emission-lines (e.g., iRubin et al.l I1980L I1982L 119851 ) and 
HI radio emission (based on various sources in the lit- 
erature). Most of Kent's sample comes from the Ru- 
bin et al. surveys, which selected for spiral galaxies with 
high inclinations, spanning a wide range of luminosi- 
ties, scale-lengths, and Hubble types, and without strong 
bars. Despite advances in observational resources in the 
intervening decades, we know of no comparable, publicly- 
available sample that includes both rotation curves and 
photometry with detailed bulge/disk decompositions for 
a wide range of disk-galaxy types. 

We estimate the disk and bulge scale-lengths (Rd and 
Qte,h) by modeling the nonparametric Kent decomposi- 
tions with simple exponential and de Vaucouleurs pro- 
files (n = 1 and n — 4, respectively). Our mod- 
els thereby treat all bulges as "classical", with n ~ 
4, neglecting some variations in their detailed proper- 
ties, such as the n ~ 1-2 indices of "pseudo" bulges 



(|Kormendv fc Kennicuttl 12004ft . The latter bulges tend 
to be much less massive, and make only minor contri- 
butions to the total for spirals, which is insensitive 
to the details of the adopted bulge density and rotation 
profiles 

For 34 of these sample galaxies (type Sb to Sc), in- 
dependent d ecompositions were car ri ed ou t on the same 
data set by Andredaki s fc Sanders! (|1994D . using para- 
metric fits to the raw surface brightness profiles. Our i?d 
values agree with theirs at the ~ 10% level, while the 
bulge results are highly variable, both between our anal- 
yses and theirs, and between different model fits by these 
authors. Most of these galaxies are very disk dominated 
(B/T < 0.1), so it is not surprising that the bulge param- 
eters would be very uncertain. Fortunately the bulges in 
such cases turn out to be only very minor contributors 
to the total of their host galaxies. Other parameters 
and their sources are listed in Table 031 

For v s of the stellar disk components of these galaxies, 
we assume that they rotate with the same velocities as 
their gas disks. We derive v c based on the rotation curves 
over the range (2-3) i?d, re-projecting this intrinsic value 
to the observed v s according to the inclination (v s — 
v c sin i). 

The final and most challenging parameter to estimate 
is the characteristic rotation velocity v s for the bulges. 
Direct estimates of bulge rotation- velocity profiles over a 
large range in radius require extensive spectroscopic data 
combined with careful bulge-disk kinematic decomposi- 
tion. As far as we kn ow, this has only bee n done for one 
spiral galaxy to date (jDorman et al.ll2012|) . Thus we are 
much worse off with estimating for spiral bulges than 
for early-type galaxies, and must make even stronger 
simplifying assumptions than in the original F83 anal- 
ysis of ellipticals. Fortunately, because the spirals are 
disk-dominated, we will find that their total estimates 
are only mildly sensitive to the assumptions about bulge 
kinematics. 

Our strategy for the bulge v s values is to estimate 
these indirectly, based on other observables: the ellip- 
ticity e = 1 — q and the central velocity dispersion a®. 
These three parameters may be related together through 
the following model: 



«• - G)' 



00 



1 -e 



1/2 



(8) 



where (v/cr)* is a parameter describing the relative dy- 
namical importance of rotation and pressure. In an edge- 
on galaxy, (v/cr)* ~ 1 represents an oblate isotropic sys- 
tem where the observed ellipticity is supported by rota- 
tion, and thi s model also turn s out to work well at other 
inclinations (|Korm cndv 1983). 

The standard lore is that spiral bulges and 
low-luminosity ellipticals ar e near oblate-isotropic, 
with t ypical (v/a) * ~ . 9 dKormendv fc Illingworthl 



19821: iDavies et al 



Binnev fc Tremaine 



19831: IBinnev fc Merrifieldl Il998fc 
20081 ). However, some con- 
cerns about these conclusions were raise d early on 
(|Whitmore et al.l 119841: iFillmore et all I1986D and mod- 

3 More extensive observations and modeling in the future could 
be used to establish the j*—Mj, trends for morphologically differ- 
ent bulges, and thereby provide physically-based information as to 
whether or not there are genuinely distinct subtypes. 
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Fig. 10. — Relation between bulge rotation velocity and velocity 
dispersion as a function of ellipticity. The points show data for 26 
spiral galaxies from the literature, with symbol shapes and colors 
corresponding to different Hubble types as in the legend. The 
curves show Equation ((SJl with (v/a)* = 1 and (v/a)* = 0.7 for 
the dotted and solid curves, respectively. We adopt (v/a)* = 0.7 
as our default model. 

ern integral-field analysis of early-types has revealed that 
their rotation velocities tend to be si gnificantly lower 
than in the oblate i s otrop ic model (Cap pellari et al.l 
120071: lEmsellem et al.l 1201 11 1. The rotation of spiral 
bulges, on the other hand, has not seen systematic inves- 
tigation in decades (s ome new work has just appeared in 
iFabricius et al.|[20T2l) . and here we attempt only a quick 
look at the implications of recent papers that have re- 
ported bulge kinematics for a handful of cases. 

We ta ke results on (v / a) and e from lLaurikainen et al 
2007), iMorelli et all (|2008f >. and iMacArthur et al 



2009) , and plot them in Figure \W\ We see that the 



oblate isotropic model is not a good representation of 
most of the data, nor is any other simple value of (v/a)* . 
However, in order to have a simplified framework for 
bulge rotation, we characterize this data set as having 
(v/a)* = 0.7 ± 0.4 (median and 68% scatter). 

We therefore adopt the following procedure for esti- 
mating bulge j*. We use the observational values for e 
and (To, and then estimate v s using Equation ([8]) with 
(v/a)* = 0.7 representing a typical value for bulges. We 
test the impact of the latter assumption on the results by 
also using (v/a)* = 0.3 and 1.1 to bracket the possible 
range of average bulge rotation. We thereby explore the 
systematic uncertainty in bulge rotation but not the in- 
trinsic scatter, keeping in mind also that this bulge model 
is based on the central regions and does not account for 
the uncertainties in extrapolating the rotation to large 
radii, as discussed in detail for the early-type galaxies. 

The e values are taken from the Kent derivations. We 
take the op measure ments in most c ases from Hyp e rLeda 
(iPaturel et al.ll2003ft and also from lCorsini et al.l (|1999Q 
and iPizzella et alT(|2004l ) . For some of the later- type 



galaxies, there are no Co measurements available, and for 
these we use an empirical relation (which we infer from 
other galaxies in these studies) that ao is approximately 
equal to the gas-disk rotation velocity. Such cases all 
have B/T < 0.15, so this approximation is not of major 
importance for the total J* estimates, but any inferences 
for these particular bulges will be relatively uncertain. 

We now have enough information to proceed with the 
specific angular momentum calculations for the spiral 
galaxies. Again, our basic approach is to estimate sep- 
arately the bulge and disk angular momenta jh and jd- 
Given a bulge stellar mass fraction quantified as /t,, we 
can then estimate the total specific angular momentum 
by: 

3 = /bib + (1 - h)jd. (9) 

In practice, we use the bulge-to-total r-band luminosity 
ratio B/T (from the series of Kent papers) as a proxy 
for f h . 

To calculate the projected values of jb and jd, we use 
Equation (J6j> . For the intrinsic values, we assume that 
both the bulge and the disk in a given galaxy have the 
same inclination i, which is estimated from the observed 
disk ellipticity. We then use the deprojection factor d to 
convert projected to intrinsic values [see Equation Q]. 
For the disk, this is a simple factor of (sini)^ 1 , and the 
calculation reduces to Equa tion Q. For the bulge, we 
calculate Ci from Equation (IA29[) . 

Using these procedures, we construct a catalog of spi- 
ral galaxies with characteristic masses, scale-lengths, and 
rotation velocities for both their bulge and disk compo- 
nents. We report these values in Table 03J along with 
the total galactic specific angular momenta (bulge and 
disk combined), both projected and intrinsic. When we 
vary the assumed bulge rotation systematically across 
the bracketing range, the total is changed by no more 
than ~ 0.03 dex (~ 7%) for the vast majority of the 
galaxies, and up to ~ 0.1 dex (~ 25%) for a few of the 
Sa-Sab galaxies. Therefore the details of the bulge mod- 
eling are of only very mild importance to the overall 
results for the spirals. These data will be used in later 
sections to examine various scaling relations for these 
galaxies and for their subcomponents. 

4.2. Early-types 

For the gas-poor early-type galaxies (lenticulars and 
ellipticals), the challenge is to assemble a large sample 
with all of the ingredients that we need to calculate 
(i.e., v s , a e , n). The information is scarcest for v s , and 
therefore we have scoured the literature for kinematic 
data sets extending to radii of at least ~2a c , assembling 
a sample that, although not exhaustive, is unprecedented 
in its size and scope. The sources include integrated- 
starlight absorption-line spectroscopy, and velocities of 
GCs and PNe. To estimate approximate values for v s , 
we simply read off the major-axis rotation velocity at 
2 a c (as explained in Section |3~3|) . We thereby assemble 
a total sample of 40 early-type galaxies, including the 8 
galaxies that we modeled in detail in Section [3] 

Table 03 provides a summary of our sample, along 
with the sources of kinematic data. Given that the data 
are drawn from a variety of literature sources with com- 
plex selection effects, it is important to check whether or 
not the sample is a fair representation of early-types in 



Angular momentum and galaxy formation revisited 



13 



the nearby universe. We have done so in Appendix [C] 
using the ATLAS 3D volume-limited sample of nearby 
galaxies as a reference, and focusing on the masses M+ 
and central rotation parameters (v/a)*. 

We find that the distribution of our sample galaxies in 
the (v / a)*-M* parameter space is fairly similar to that of 
an unbiased sample over a similar mass range. The me- 
dian galaxy mass in our sample is log (M*/Mq) — 10.8, 
whi ch is near the cha racteristic mass M* of nearby galax- 
ies (|Guo et al.ll201Q[ ). We thus conclude that our obser- 
vational results should be representative of low-redshift 
ordinary early-type galaxies. The only caveat here is 
that our sample is biased toward ellipticals at the ex- 
pense of lenticulars, which we must take into account 
later when drawing conclusions about the overall popu- 
lation of early-type galaxies. 

An alternative scheme for classifying early-types is 
as "fast rotators" (including almost all lenticulars) 
and "slow rotators" , base d on their central kinemat- 
ics (|Emsellem et al.l I2007D . The central rotation is 
known to correlate with many other galaxy p roperties 
(jDavies et alJll983t iKoTmendv fc Benderlll996[ ). and the 
fast and slow rotators have been interpreted as having 
different formation histories. Therefore it is important 
that we investigate to what extent the global specific an- 
gular momentum j± correlates with the central rotation 
classification. Our sample includes three slow-rotators, 
which is consistent with the fraction of su ch galaxies in 
the nearby universe (jEmsellem et al.ll2011[ ) , and will pro- 
vide a rough initial idea of any systematic differences 
between fast and slow rotators. 

Returning to the remaining observational parameters, 
for each early-type density profile, we need both the 
Sersic index n and the corresponding scale-length a e 
(which can differ significantly from the value obtained 
with a classic n = 4 fit, e.g., in the RC3 catalog of 
Ide Vaucouleurs et al.lll99l"l ). Unfortunately there is no 
comprehensive source available for such measurements, 
and we resort to a medley of literature data. 

For 34 of the galaxies in our sample, there are pub- 
lished Sersic fits, and we take the (a c ,n) values accord- 
ing to the following priority: detailed photometric anal- 
ysis i n ind ividual galaxy papers ( e.g., iNapolitano et al] 
120091) ; t he IKormendv et all (120091) tabulation for Virgo 



galaxies: IHopkins et al.l (|2009an : IHopkins et al.l (|2009 
ID'Onofriol (|2001h . 

For the remaining 6 galaxies, we have as a starting 
point the RC3 value for the effective radius. Then 
we use the well-established observation that there are 
strong correlations between early-type gala xy size and 
luminosity, and the Sersic index n (e.g.. iCaon et al 
1993HPrugniel fc Simierill997tlGraham fc Guzmanl l200 
Blanton et al.l 120031 : IKormendv et alJl2009f ). This allows 



us to estimate a most-probable n value for each galaxy 
(see Appendix [C] for details) . 

Note that if we were simply to approximate all of 
the early-types as n = 4 spheroids, the k n values in 
Equation (j6]) would be too high on average by ~ 30% 
(~ 0.15 dex, given a median index value of n ~ 2.5). This 
would translate to an equivalent systematic error on j*. 
We could adjust for this effect by adopting n = 2.5 in all 
cases, but n also has a systematic dependence on galaxy 
mass, and ignoring this fact would produce a spurious 
mass-dependent trend in j± of ~ 50% (~ 0.2 dex) over 



the full range in mass. 

In Table Cj5j we compile the observed parameters v s , 
a e , and n for our full early-type galaxy sample. We use 
these to calculate j p approximately from Equation ([6]), 
and tabulate these values as well. For some of the very 
extended galaxies like NGC 4374, the total luminosity 
and angular momentum (via the factor k n ) are integrated 
out only to the estimated virial radius. 

In order to convert projected j p to intrinsic jt for 
analysis in later sections, we must apply a deprojec- 
tion factor Cj which depends on the inclination i. Un- 
fortunately, the individual inclinations are not generally 
known, but neither are they completely random, because 
of an inclination-bias in galaxy classification. As dis- 
cussed in Appendix IA. 31 we therefore apply median de- 
projection factors of C mc d = 1-21 (+0.08 dex) to the 
lenticulars, and C mc d = 1-65 (+0.22 dex) to the ellipti- 
cals. 

Since one of our eventual goals will be to quantify the 
intrinsic scatter in the observed j+—M* relations, it is im- 
portant to be clear about the error budget in our analy- 
ses. Again, the basic parameters that go into our j* cal- 
culations are Ci, a e , n, and v s . For early- type galaxies 
with an assumed n = 4 profi le, the typical uncertain ties 
in a c are - 25% (~ 0.1 dex: iCappellari et al.|[20llab . If 
we allow for a more general n, which for some galaxies 
is measured directly and in other cases is derived sta- 
tistically (Appendix [Cj, then we estimate a combined 
uncertainty on j* from a e and n of ~ 40% (~ 0.15 dex). 
The uncertainty on v s from our simplified measurement 
and extrapolation approach is ~ 25% (~ 0.1 dex; Sec- 
tion HT3J)- 

Table [T] summarizes the uncertainties introduced by a 
number of different ingredients in the calculations. 
The separate uncertainties for and M* are mapped to 
the direction perpendicular to a cx M^ 3 trend, as dis- 
cussed in Section [3.5l This net uncertainty is designated 
AA, owing to the connection with spin-based theoretical 
models. 

The total uncertainty in A for late-type galaxies is typ- 
ically ~ 30% (~ 0.1 dex), and is driven by the estimate of 
M* (via T*) rather than j*. For the vast majority of the 
early-types (apart from the special class of massive, ex- 
tended ellipticals), the uncertainty is ~ 60% (~ 0.2 dex), 
and is driven by the four parameters mentioned above 
that enter into the j+ calculation. 

This full j*-M* dataset is assembled from a generally 
unbiased ~ M* galaxy sample that we can use to investi- 
gate differences in angular momentum not only between 
early-types and spirals, but also between ellipticals and 
lenticulars, and between fast and slow rotators. 

4.3. Size and rotation-velocity scaling relations 

Before considering specific angular momenta and their 
correlations in the next section, we examine some trends 
among the raw ingredients that go into these analyses, 
a c , v s , and M*. Doing so provides a check that our 
results are consistent with the familiar size-mass and 
mass-rotation velocity (Tully-Fisher) relations that have 
been established for nearby galaxies. We also introduce 
novel relations involving rotation, and explore some pre- 
liminary indications about angular momentum. 

We first consider the standard scaling relation of galaxy 
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TABLE 1 
Uncertainty budget 



Galaxy 








AA (dcx) 










type 


D 


c, 


v s 


Va 


n, a e 


bulge 


X* 


total 


Sb-Sm 


0.01 


0.01 


0.02 


0.03 


0.05 


0.03 


0.07 


0.09 


Sa-Sab 


0.01 


0.01 


0.02 


0.03 


0.05 


0.1 


0.07 


0.13 


SO 


0.01 


0.05 


0.06 


0.1 


0.15 





0.07 


0.18 


fE 


0.01 


0.15 


0.06 


0.1 


0.15 





0.07 


0.22 


sE 


0.01 


0.12 


0.35 


0.35 


0.2 





0.2 


0.40 



2004; 



2011 



Gadottill2009]lMaltbv et al.l l2010t iMcDonald etal 
Dutton et alj 120111) . Given the many different as 



sumptions and data sources that went into our sizes and 
masses, these parameters match the ATLAS 3D results re- 
markably well overall (with some nuances discussed fur- 
ther in Appendix [C|. This suggests that our size and 
mass data are representative and reliable at the ~ 0.1 dex 
level. 

We can also consider separately the spiral bulges, 
plotting their sizes and masses for our sample in Fig- 
ure [TTJ Although the full range of sizes is not vis- 
ible in this plot, the bulges follow a roughly parallel 
size-mass relation to the elliptical galaxies, but smaller 
on average by a factor of ~ 4 (~ 0.6 dex) and with 
a great deal of scatter (possibly because of the ap- 
proximate nature of these size measurements). Other 
studies have als o found that bulges are more compact 
than ellipticals (IGraham fc Worlevll2008t IGadottil 120091: 
lLaurikainen et all 120101: iDutton et al.1 1201 lft , "but the 
quantitative details vary considerably, and we therefore 
regard our bulge scaling relations as provisional. 

The next scaling relation that we consider is rotation 
velocity versus mass. For spiral galaxies, this is the Tully- 
Fishcr relation, but it has to our knowledge never been 
constructed previously for all galaxy types. We can al- 
ready generate a broad expectation for what we will find, 
given the observed size-mass relations along with the 
assumption that is independent of galaxy type. As 
mentioned in Section [2j we can then use Equation ([6j) to 
predict the ratio of characteristic rotation velocities for 
ellipticals and spirals: 



Vs,E fcl Qe.Sp 

VsSp k A a CjE ' 



(10) 



where we are approximating the spiral galaxy parameters 
as dominated by the disk component. With fci/fej = 0.5, 
and a e ,Sp/ie,E ~ 2 for our sample, we therefore predict 



10 



Note. — The uncertainties on and have been converted 
into equivalent uncertainties on A. The different galaxy types in- 
clude fast- and slow-rotating ellipticals (fE and sE). The listed 
sources of potential error are distance (D), corrections for projec- 
tion effects including inclination (C;), the rotation velocity scale 
calculated in detail (v s ), the alternative approximate rotation ve- 
locity scale (v B ), the stellar density profile Sersic index (n) and 
scale radius (a e ), the incorporation of bulge contributions, and the 
stellar mass-to-light ratio including IMF variations (T*). 



size versus mass, or a c versus M* in our notation, show- 
ing the results in Figure 111! where we again compare 
our results to the volume-limited ATLAS 3D sample as 
a baseline check. We find that in both samples, late- 
and early-type gal axies have rough l y the same sizes at 
a given mass (cf. IShen et aLl I20031 12007T) , but there is 
a clear systematic trend for the more bu lge-dominated 
galaxies to be more compact (see also Ide Jong et al 
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Fig. 11. — Relation between size and stellar mass for our galaxy 
sample. The former is the semi-major axis effective radius, and 
the latter is based on if-band total luminosities with an adopted 
mass-to-light ratio of M^/Lpc = 1 in solar units. Different symbols 
denote different galaxy types as shown in the legend; for the spirals, 
the disk and bulge ("B") components are shown separately. The 
range of the plot is restricted in order to better see the main trends 
in the data; the bulge data extend to radii as small as a c ~ 0.01 kpc 
(note also that the most compact elliptical shown is NGC 4486B, 
which is considered a rare, highly-stripped galaxy). For compari- 
son, diagonal lines show power-law model fits to the data from the 
ATLAS 3D survey (i.e., independent from our data set): lenticulars 
and fast-rotator ellipticals (dot-dashed), Sa-Sb spirals (dashed), 
and Sc-Irr spirals (dotted). For both data sets, the late- type galax- 
ies are systematically larger than the early-types at a given stellar 
mass. The absolute normalizations of the trends are similar be- 
tween the ATLAS 3D sample and ours, with some small differences 
as discussed in the text. 

v s ,e/vs,Sp ~ !■ Thus, ellipticals should rotate at roughly 
the same velocity as spirals if they have the same specific 
angular momenta at a given mass. 

Without proceeding any further, this scaling analy- 
sis already suggests that ellipticals have lower j+ than 
spirals, or else they would be extremely flattened by 
rotation, similarly to the spiral disks which have near- 
maximal rotational support (modulo possible differences 
in dynamical mass between spiral and elliptical galaxies 
at the same stellar mass). The same argument applies 
even more strongly to the spiral bulges, since they are 
far more compact than the disks at a given mass. If the 
bulges had the same j* as the disks, then they would 
have to rotate much faster, which is impossible. 

We now examine what our new collection of observa- 
tions tells us directly about the rotation scaling relations. 
The left-hand panel of Figure [T2l shows the characteristic 
rotation velocity v s for the elliptical and lenticular galax- 
ies, and the spiral disk and bulge subcomponents, in our 
sample. Here we are plotting the intrinsic rotation ve- 
locity, multiplying by the deprojection factor C%, which 
is jus t (sin i)^ 1 for disks (see Appendix I A. 2 [) . and Equa- 
tion (|A29[) for bulges. For the early-type galaxies, the 
inclinations are unknown, and we have adopted median 
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Fig. 12. — Relations between characteristic rotation velocity d v s , st ella r mass (left-hand panel), and size (right-hand panel) for our full 
galaxy sample, using the same data sources and symbols as in Figure 1111 For the spiral disks, v s is the outer gas-disk rotation velocity. 
For the lenticulars and ellipticals, v s is the stellar rotation velocity measured along the semi-major axis at 2 a e , except for the points with 
error bars, which are the eight cases studied in detail in Section [3] with v s derived from full modeli ng o f the rotation- velocity profiles. For 
the bulges, v a is estimated indirectly using flattening and velocity dispersion observations (Section 14.11 1. In all cases, the rotation velocity 
has been deprojected for both inclination and "dilution" effects, using the factor d (see text for details). In the left-hand panel, the dotted 
blue line shows a least-square fit to the Sb-Sc disks, a dashed red line shows a proposed inverse trend for a subset of the E/SOs, and the 
blue dot-dashed line shows the baryonic Tully-Fisher relation for late-type galaxies from Truiillo-Gomcz ct al. ( 2 011) for comparison. In 
the right-hand panel, the diagonal line shows a prediction for the spiral disks based on ACDM models (see Section l6.2|l . Overall, the spiral 
and elliptical galaxies follow mass— rotation velocity and size— rotation velocity trends that have remarkably opposite slopes. The trends for 
the lenticulars are between the spirals and ellipticals. 

ity and size (right-hand panel) that we will discuss later 
in this paper. 

The behavior of the lenticulars in the mass-rotation ve- 
locity diagram is difficult to discern in detail owing to the 
small sample size, but in general it appears intermediate 
to the other galaxy types. We also notice an interest- 
ing pattern when considering the lenticulars and ellipti- 
cals together: there may be a bimodal mass-rotation ve- 
locity relation^ with some galaxies following the trend 
for spirals, and others following a steep reverse relation, 
Ci v s oc M~ 3 . The implication is that there may be two 
distinct populations of early-type galaxies, one of which 
is closely related to spirals, and which are not equivalent 
to standard E and SO classifications. 

The bulge rotation velocities appear to follow a similar 
trend to the spirals, at about half the amplitude. Here it 
should be remembered that the bulge "data" points are 
indirect estimates constructed in order to provide plau- 
sible adjustments to the total angular momenta of the 
spiral galaxies (Section l4.ip . The results so far suggest 
that bulges are different from ellipticals in their mass- 

5 This pattern may be partially an artifact of inclination effects. 
In particular, some of the edge-on lenticulars were observed with 
long-slit spectroscopy directly along their embedded disks, which 
may not provide an accurate measurement of the overall rotation. 
However, for the ellipticals we find no correlation between apparent 
rotation velocity and ellipticity. An additional issue is that the 
occasional extremely low-inclination galaxy will not be treated well 
by our me dian-dcprojection method (cf. the right-hand panel of 
Figure |24) |. so in any fits to the data, we will discard outliers with 
very low v B or (e.g., NGC 1419). 



factors for C, as discussed in Section [ 

We see that the disks follow a fairly tight relation 
of approximately Ci v s cx Mj~ , with a residual trend 
for the later-type disks to rotate more slowly. This is 
equivalent to the familiar Tully-Fisher relation, and in 
the Fi gure we include a recent result from the liter- 
ature (Trujillo- Go mez et al.l l20"Tll ). which matches our 
data ver y well (cf. the type -dependence among spirals 
found bv lMasters et aLl feOOS) . We also show in the right- 
hand panel of Figure [12] the relation between size and ro- 
tation velocity, which are strongly correlated parameters 
for disk galaxies. 

The elliptical galaxies are completely different, show- 
ing an anti- correlation between rotation velocity and 
mass0 with Ci v s oc M~ 01 . This result also con- 
trasts markedly with standard relations for ellipticals 
involving the velocity dispersi on op or the dynami- 
cal mass (e.g., op oc M ° 25 ; iFaber fc Jackson! 119761 : 
iTrujillo-Gomez et al.ll201lT ). In galaxy disks, the rota- 
tion velocity traces the dynamical mass, so the Tully- 
Fisher relation is a measure of both mass and angular 
momentum. In elliptical galaxies, on the other hand, the 
mass and angular momentum relations are decoupled. 
We also find an anti-correlation between rotation veloc- 



4 This echoes a similar trend in the central rot atio n properties 
of early-type galaxies in general (shown in Figure |33). The e ight 
galaxies studied in detail (points with error bars in Figure I12H 
are consistent with this trend but do not include enough lower- 
luminosity ellipticals to distinguish between v B being constant or 
decreasing with mass. 
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size-rotation velocity relations, and we will see in the 
next section how their angular momenta compare. 

Since both the sizes and the rotation velocities of el- 
liptical galaxies are systematically lower than for spiral 
disks, we can already predict that the ellipticals will on 
average have much lower j*. Note that although this 
conclusion has already been widely adopted for decades, 
only now have the kinematic data reached large enough 
radii to confirm it with confidence. 

To see that the low characteristic rotation velocities 
for ellipticals are not a mathematical sleight of hand, 
one may consider the specific cases of NGC 821 and 
NGC 3377 in Figure |H1 The rotation- velocity profiles of 
these galaxies decline dramatically outside x ~ (1-2) a e , 
which may be contrasted with the spiral galaxies in Fig- 
ure |U Preliminary analysis of additional edge- on cases, 
where the deprojection uncertainties are minimized, in- 
dicates that such declines are a generic feature of ~ M* 
early- type galaxies (A. Romanowsky et al., in prepara- 
tion). 

This conclusion includes NGC 2768, which from the 
current data appears consistent with a constant or ris- 
ing outer rotation velocity, but which with more ex- 
tensive new PN data may have a declining outer pro- 
file. Even the cases of strongly rising rotation-velocity 
profiles out to x ~ 2 a c found by iRix et al.l (fT9991 ) 
appear upon closer inspection to turn over at larger 
radii. These results all contrast with early claims of 
high outer rotation in some early-types, which were 
recently overturned w i th improved observations (e.g. , 
Arnaboldi et all [l99l iKissler-Patig fc Gebhardtl \wm 
Romanowskvl 120061: iMcNeil et al.l l2010t iStrader et al I 
2011). 

We can also begin making some interesting inferences 
about the relations among other galaxy types, based on 
both size and rotation- velocity trends (Figures QT] and 
[T2"|) . As discussed, the lenticulars share similar properties 
to spirals in some cases, and to ellipticals in others. The 
distinction between "fast" and "slow" rotator ellipticals 
based on their inner regions does not appear to hold up 
when considering their global rotation properties. 

This overview of the observable scaling relations be- 
tween mass, size, and rotation velocity gives us a preview 
of some of our overall conclusions about angular momen- 
tum, and provides more confidence in the solidity of those 
conclusions. We construct a novel mass-rotation veloc- 
ity relation for ellipticals, which is the analogue of the 
Tully-Fisher relation for spirals, but with the remark- 
able difference of having a negative slope. The data also 
imply that both elliptical galaxies and spiral bulges must 
have lower specific angular momenta than spiral disks of 
the same mass. We address this issue more quantita- 
tively in the next section, incorporating the additional 
mass-dependent factor k n in calculating . 

5. OBSERVATIONS: ANGULAR MOMENTA OF THE FULL 
SAMPLE 

Having derived estimates of the j* and M* parameters 
for our full galaxy sample, we now examine the resulting 
observational trends, which constitute the key results of 
this paper. We begin by focusing on the late-type galax- 
ies in Section lSTTl and combine these with the early-types 
in Section [5.21 We discuss our proposed replacement for 
the Hubble sequence in Section 15.31 which we test by 
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Fig. 13. — The total (disk plus bulge) stellar specific angular 
momentum of nearby spiral galaxies plotted against total stellar 
mass. The top and bottom panels show estimates of projected 
and intrinsic jf* , respectively; the uncertainty in j+ for each galaxy 
is in almost all cases smaller than the plotted symbols. Different 
symbols denote galaxy sub-types as specified in the legends. The 
dotted lines show fits to the data in each panel, while the dashed 
lines show fits to the disk components alone (data not shown). 
The spiral galaxies follow a universal j*-M* relation, with some 
dependence on Hubble type. The projected relation is very sim- 
ilar to the intrinsic relation, but with a small offset, and slightly 
increased scatter, in j*. 

examining systematic residuals from the j*-M* trends 
in Section 15.41 We further convert the j*-M* data into 
one-dimensional histograms in Section [5.51 

5.1. Lessons from spirals 

Although the main novelty of this paper is our care- 
ful consideration of early type galaxies, we also include 
the oft-studied category of spirals in order to provide an 
integrated analysis of bright galaxies of all types. Fur- 
thermore, the well-constrained angular momenta of the 
spirals also permit us to better understand systematic 
issues such as inclination corrections that are trickier to 
handle for early-types. 

We plot the total (disk+bulge) j+-M+ data for the spi- 
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TABLE 2 

Mass-angular momentum fits to data 



Sample 



log JO 



Q 



0~lo 



All spirals, total, projected 3.11 ± 0.03 0.53 ± 0.05 0.22 

All spirals, total, intrinsic 3.18 ±0.03 0.52 ± 0.04 0.19 

Sa-Sab, total, projected 2.93 ± 0.05 0.60 ± 0.06 0.17 

Sa-Sab, total, intrinsic 3.02 ± 0.04 0.64 ± 0.07 0.12 

Sb-Sbc, total, projected 3.15 ±0.03 0.65 ±0.14 0.16 

Sb-Sbc, total, intrinsic 3.21 ± 0.03 0.68 ±0.13 0.15 

Sc-Sm, total, projected 3.25 ± 0.04 0.58 ± 0.06 0.20 

Sc-Sm, total, intrinsic 3.29 ± 0.04 0.55 ± 0.05 0.18 

All spirals, disks, projected 3.25 ± 0.02 0.62 ± 0.05 0.20 

All spirals, disks, intrinsic 3.31 ± 0.02 0.61 ± 0.04 0.17 

Sa-Sab, disks, projected 3.25 ± 0.05 0.76 ± 0.09 0.21 

Sa-Sab, disks, intrinsic 3.34 ± 0.04 0.82 ± 0.08 0.17 

Sb-Sbc, disks, projected 3.24 ± 0.03 0.71 ±0.14 0.16 

Sb-Sbc, disks, intrinsic 3.30 ± 0.03 0.75 ±0.12 0.13 

Sc-Sm, disks, projected 3.29 ± 0.05 0.61 ± 0.07 0.21 

Sc-Sm, disks, intrinsic 3.33 ± 0.05 0.57 ±0.05 0.19 

All spirals, bulges, projected 2.20 ± 0.31 0.69 ±0.11 0.58 

All spirals, bulges, intrinsic 2.32 ± 0.31 0.69 ±0.10 0.57 

Sa-Sab, bulges, projected 2.30 ± 0.32 0.99 ±0.15 0.47 

Sa-Sab, bulges, intrinsic 2.44 ± 0.32 0.99 ±0.15 0.46 

Sb-Sbc, bulges, projected 1.89 ± 0.34 0.34 ± 0.20 0.58 

Sb-Sbc, bulges, intrinsic 2.01 ± 0.33 0.34 ±0.19 0.56 

Sc-Sm, bulges, projected 2.21 ± 0.57 0.64 ± 0.27 0.60 

Sc-Sm, bulges, intrinsic 2.30 ± 0.58 0.63 ± 0.28 0.60 



Lenticulars, projected 
Lcnticulars, intrinsic 
Ellipticals, projected 
Ellipticals, intrinsic 



2.97 ±0.08 
3.05 ± 0.08 
2.52 ±0.05 
2.73 ±0.05 



0.80 ±0.14 
0.80 ±0.14 
0.60 ±0.09 
0.60 ±0.09 



0.29 
0.29 
0.24 
0.24 



Sb-Sm, intrinsic, fixed a = 2/3 3.28 ± 0.03 0.67 0.19 

Ellipticals, intrinsic, fixed a = 2/3 2.75 ± 0.05 0.67 0.24 

KCDM halos 2.50 0.67 0.23 



rals from Table CHI in Figure [13] In the top panel, we 
show the projected value, j p , and in the bottom panel, 
the intrinsic value, jt- These are related trivially by the 
disk inclination, but we wish to investigate how well the 
trends in projection reflect the intrinsic trends, since de- 
projection for the early-type galaxies will be more diffi- 
cult. 

Overall, the spiral galaxies appear to follow fairly tight 
j+-M* trends, with similar slopes, regardless of Hubble 
sub-type. In more detail, we carry out least-square fits 
to j+ as a function of M± in log-log space: 



log jmod = log jo + a [log(M*/M ) - 11] 



(11) 



with a residual rms scatter that we parameterize as 
ciogjj,- The uncertainties in the fit parameters jo and 
a are estimated by bootstrap resampling. 

Our fitting results for various spiral subsamples are re- 
ported in Table [2] For total j*, the systematic uncertain- 
ties from the bulge rotation (see Section l4~Tj) turn out to 
be smaller than or equal to the statistical fitting uncer- 
tainties, even for the Sa-Sab galaxies, and in the Table 
we have combined both uncertainties in quadrature. 

The data are basically consistent with a universal j*- 
M* slope for spiral galaxies of all types, with a ~ 0.6 
and an rms scatter of o~\ og j ~ 0.2 dex. There is also a 
clear residual trend with Hubble type: the Sa-Sab galax- 
ies have systematically lower j* than the Sb-Sm galax- 
ies. These conclusions hold for both j p and jt, although 
the uncertainties and the scatter are smaller for j t , as 



expected if there are genuine, underlying physical corre- 
lations that become clearer after deprojection. 
_ The multi-component nature of our model galaxies al- 
Jpws us to look further at disk and bulge properties sep- 
arately. We will take up this issue in Section [521 an d f° r 
now provide the fits to the jd-M^ and jb—M^ relations in 
Table [5] It should be remembered that the bulge results 
depend on model assumptions, although as discussed, we 
have plausibly bracketed their upper and lower limits for 
j*- 

As anticipated, the bulges turn out to have little im- 
"pact on the total j* trends for the Sb-Sm galaxies, which 
are dominated by the disk components. For the Sa-Sab 
galaxies, the bulges are responsible for the systematic 
offset with respect to the later types; this offset changes 
slightly but persists when adopting the upper or lower 
limits to the bulge rotation. The disks of all the galaxy 
types turn out to follow nearly the same j^-M* relations. 
_ This analysis demonstrates that inclination effects are 
not expected to have a major impact on our overall re- 
sults, since for both disks and bulges, the intrinsic and 
projected j+-M+ trends as well as their scatter are very 
similar. There is an overall offset between disk jt and j p 
of ~ 0.07 dex, which is comparable to the range of 0.04- 
0.06 dex that we would expect, given the median incli- 
nation i = 67° of our sample, and depending on whether 
The j^—M* trend represents a median or an average fit 
(see Appendix IA.2I for further discussion) . 

For our ensuing study of early- type galaxies, we will 
.therefore simply adopt median deprojection values for 
all of the galaxies, which we estimated in Section |4~21 to 
mean adding offsets of 0.08 dex and 0.22 dex to j p to 
-derive jt, for lenticulars and ellipticals, respectively. We 
can also in general drop the usage of j p in the rest of 
this paper, in favor of the more physically meaningful j t 
which we now adopt as our estimate for j*. 

5.2. Combined observational results 

We are now ready to include the early-type galaxies in 
our analysis, and thereby address most of the key science 
questions raised in Section [TJ As a reminder, our start- 
ing point is the j*— M ± diagram from F83 that we have 
reproduced in Figure [T] Do we find the same j+—M+ 
trends with an updated and expanded dataset, and more 
detailed analysis? Do ellipticals still appear to have sys- 
tematically low j± relative to spirals, or do we discover 
large reservoirs of additional j* at large galactocentric 
radii, using modern data? Do Sa and SO galaxies fill in 
any "gap" between spirals and ellipticals, and can we 
then connect the Hubble sequence to a sequence in j*? 
Can we characterize all galaxies as combinations of disks 
and bulges that follow universal scaling relations? (The 
main remaining question that connects to galaxy forma- 
tion theory will be pursued in the next section.) 

Taking our early-type galaxy j* and M+ estimates from 
Table C[5] (after statistically correcting projected to in- 
trinsic quantities; see Table [1] for an error analysis), we 
plot them in Figure [Ml (left), along with the spirals re- 
sults discussed in Section 15.11 This new Figure is the 
centerpiece of our paper. Focusing first on the elliptical 
galaxies, our basic finding is that they follow a j^-M* 
trend which is roughly parallel to the spirals but with a 
large systematic offset to lower j-,, . 

We thereby confirm the conclusions of F83, finding 
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Fig. 14. — Left-hand panel: The total intrinsic specific angular momentum of galaxies plotted against their total stellar mass. Symbols 
show galaxy types according to the legend at the upper left. The points with error bars shown are based on the more detailed j* estimator 
[Equation J3}]; for the remainder of the galaxies, the approximate estimator [Equation J6}] was used. The uncertainties are similar in 
both cases. The deprojection from observed j p to intrinsic jt was accomplished using individual inclinations for the spirals, and median 
deprojection factors for the lenticulars and ellipticals (see main text). The least massive early- type galaxy in the sample is the compact 
elliptical NGC 4486B, which is probably in the process of being tidally stripped by the giant galaxy M87; the other low-j* outlier is 
NGC 1419. Both are marked with black X symbols and excluded from all fits in this paper. Dotted lines show the best fits for the Sb-Sm 
and elliptical galaxies: these two galaxy types follow j^—M* trends that are parallel but separated in j+ by ~ 0.5 dex. Right-hand panel: 
As left-hand panel, but now plotting spiral disks and bulges alone, along with elliptical galaxies, as indicated by the legend. The upper 
line is now the fit to the disks (for all spiral types) rather than to the whole galaxies. Note that the slopes of the lines in this panel and the 
left-hand one should not be compared by eye, owing to the different axis ranges. The uncertainties in j+ for the disks are typically ~ 0.04 
dex, and for the bulges at least ~ 0.2 dex; the M* uncertainties are systematic (see main text). Many of the most massive spiral bulges 
appear to a follow a similar j*—M* relation to the ellipticals. 



from a new synthesis of modern photometric and kine- 
matic data that the "missing" angular momentum in el- 
lipticals does not emerge at large radii, as had been ex- 
pected from some theoretical studies. As discussed in 
Section l4~3l the new observations tend to show outer ro- 
tation profiles that decline rather than rise. Even the 
nearby galaxy NGC 5128 (Cen A), which is often con- 
sidered to be an elliptical formed through a recent major 
merger, shows a relatively low j* when compared to spi- 
rals of the same stellar mass. Whether or not these obser- 
vations pose a genuine problem to major-merger expla- 
nations for forming ellipticals will require renewed the- 
oretical analysis, but as discussed in Section 14.31 there 
seems to be a pattern in the literature of misdiagnoses of 
high outer rotation from early, sparse data - which led 
to premature claims of evidence for major mergersQ 

The specific angular momentum difference between 
spirals and ellipticals is also apparent from a simple, di- 
rect consideration of the data in Section 14.31 where the 
smaller sizes and rotation velocities for ellipticals sug- 



6 INorris et al.1 (120121 ') also recently noted an emerging trend 
for low rotation in elliptical-galaxy halos, at odds with major- 
merger expectations. One possible counter-example is the SO 
galaxy NGC 1316, which is generally thought to be a major- 
merg er remnant. Based o n the new PN kinematics results 
from McNcil-Movla n et al.l j20i"2T ). we confirm the finding of 
IArnabol*dT"eTlrtT!ri99*8l) that the j+— M* values for this galaxy are 
close to the mean trend for spirals. However, we caution that our 
photometric parameters and Y* value are particularly insecure for 
this galaxy. 



gested that they have lower j*. As an arbitrary bench- 
mark, we use the median j+ at the L* characteristic lu- 
minosity, which is log (L* K / Lk.q) ~ 11, corresponding 
to log (M*/M Q ) ~ 11. For ellipticals and Sb-Sm spi- 
rals, we find projected values of j p ~ 330 km s _1 kpc 
and ~ 1600 km s~ x kpc, respectively, and true values of 

= j t ~ 540 km s _1 kpc and ~ 1800 km s _1 kpc. 

In more detail, we report fits to the j+-M+ data to- 
ward the end of Table [2] The fitted slope for the ellip- 
ticals is consistent with that for the Sb-Sm spirals, but 
is significantly offset to lower by a factor of ~ 3.4 
(~ 0.5 dex). These findings are consistent with F83, ex- 
cept that the gap has narrowed from a factor of ~ 6 
(~ 0.8 dex)H Note that if the K-h&nd for the ellipti- 
cals were systematically higher than for the spirals by a 
factor of ~ 2 (perhaps owing to age or IMF differences; 
cf. Section [23]), then the offset would increase to a 

7 Our revised Sb-Sm relation is ~ 0.1 dex lower than in F83, 
partly owing to the inclusion of bulges, and partly to new estimates 
for disk sizes and mass-to-light ratios. Our revised ellipticals re- 
lation is ~ 0.2 dex higher than in F83; this difference appears 
to arise not so much from the rotation data (the extrapolations 
to large radius by F83 turn out very good on average), but from 
a refined treatment of the total angular momentum calculation 
for spheroids. Our slopes of a = 0.53 ± 0.04 and 0.60 ± 0.09 for 
the Sb-Sm and elliptical galaxies are shallower than the a = 0.75 
slope suggested by F83; for the Sb-Sm galaxies, this difference is 
driven mostly by our inclusion of bulges and of lower-mass galaxies 
[log (M* /Mo) ~ 9]; while for the ellipticals, a shallower slope was 
already apparent in F83. 
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factor of ~ 5 (~ 0.7 dex). 

The scatter of criogj, = 0.24 dex for the ellipticals is 
similar to the j p scatter for the spirals. We also note that 
the general trends for the ellipticals are supported by the 
small sample of galaxies that we modeled in detail (see 
points with error bars in Figure UM left). Although one 
might still have concerns that large formal uncertainties 
in ji, remain for most of the sample after extrapolating 
their rotation- velocity profiles beyond 2 i? G , in order to 
close the j* S a P between spirals and ellipticals, the ro- 
tation velocity would have to rise rapidly by a factor of 
~ 4 outside these radii, which seems implausible (cf. Fig- 
ure [5]). 

The parallel nature of the spiral and elliptical trends 
is an interesting and non-trivial result, since Figure [12] 
showed that the slopes of the rotation-velocity scaling 
relations for these galaxies have opposite signs. Some 
mass-dependent conspiracy of size, rotation velocity, and 
Sersic index must be at work in order for the j*-M* 
slopes to turn out the same. 

The few "slow rotator" ellipticals in our sample show 
no indication of deviating systematically from the over- 
all fa-M* trend for ellipticals, which disagrees with 
earlier findings of m uch lower for such galaxies 
(|Bender fc Nietol Il990l ). Although their outer regions, 
like their central parts, rotate slowly relative to most of 
the fast rotators (Figure [T2l , we find that this is com- 
pensated for by their larger scale radii and Sersic indices 
(keeping in mind that the results for these galaxies are 
the most uncertain). Thus the global j ± measurements 
suggest that the slow and fast rotators may have more 
in common than was previously suspected. 

Having confirmed the basic observational findings of 
F83, we now move on to fresh territory, beginning with 
the inclusion of Sa and SO galaxies in Figure [14] (left). 
F83 suggested that these would fill the gap in j+-M+ 
space between ellipticals and late-type spirals, which is 
confirmed by our sample. Both of these galaxy types 
are on average offset to lower from the Sb-Sm spirals 
trend by a factor of ~ 1.8 (~ 0.25 dex; we will discuss 
variations about the average in Section l5.4j) . 

One natural interpretation of this new finding is that 
the Hubble classifications are related to an underlying 
physical structure, where all galaxies are composed of 
some combination of two basic components: a disk and a 
spheroid (as illustrated schematically in Figure [2] of Sec- 
tion [TJ. These components would define two distinct se- 
quences in the j^-M* plane, which in combination would 
move the total values of galaxies to intermediate regions 
in this plane, depending on the bulge-to-total mass ra- 
tios, B/T. 

To explore this idea, we plot the j^—M* data separately 
for elliptical galaxies, and for spiral disk and bulge sub- 
components, in the right-hand panel of Figure 1141 The 
disks follow a similar relation to spiral galaxies overall, 
since these are dominated by their disks. More remark- 
ably, the j^-M* trend for bulges is fairly similar to the 
trend for ellipticals over the mass range where they over- 
lap0 This is a surprising result, because as shown in Fig- 

8 At lower bulge masses, the apparent tendency to relatively low 
j+ values should be viewed as speculative, since it is based on clas- 
sical bulges rather than the pseudo-bulges that may predominate 
in this regime. 



ure llli the bulge sizes are systematically smaller than the 
ellipticals, and thus their rotation velocities (Figure [12]) 
must be higher, in an apparent conspiracy to produce 
roughly the same j*. 

A similar analysis could in principle be carried out for 
the fast-rotator ellipticals, since they are widely consid- 
ered to host hidden, embedded disk-like components. Do 
the disk and bulge subcomponents of ellipticals follow the 
same j*-M+ relations as those of the spirals? We have 
investigated this question in Appendix [D] using decompo- 
sitions from the literature, but the results are somewhat 
ambiguous. Thus, although we have been able to address 
all of the major questions raised initially about empirical 
j*-M+ trends, we flag the trends for the subcomponents 
in ellipticals (and lenticulars) as an important aspect re- 
maining in need of clarification. 

5.3. Replacing the Hubble diagram 

The foregoing discussion brings us to the diagram that 
we have already introduced schematically with Figure [2] 
which constitutes our own, physically-motivated, substi- 
tute for the classic Hubble tuning fork, and which could 
provide the underlying explanation for the observational 
trends found in Figure HH In this scheme, all galaxies 
are composed of a disk and a bulge, each adhering to a 
distinct and parallel j±-M± scaling relation. If the disk 
and bulge relations are universal (which we will further 
test in Section I5.4|) , then the location of a galaxy in j ± - 
M+ space can immediately be used to infer its B/T value 
uniquely, and vice- versa (i.e., there is a coordinate trans- 
formation between the two parameter spaces). Elliptical 
galaxies would then be the cases with B/T ~ 1, and 
bulges could be thought of as mini-ellipticals. 

As with the original Hubble diagram, our j+-M+ di- 
agram provides a simple description of galaxies, along 
with the temptation to interpret it as some kind of evo- 
lutionary sequence. However, our diagram differs, since 
the parameters used are physical quantities that may in 
principle be conserved, and thus it is actually justified 
to begin using the diagram directly as a tool to motivate 
and test some evolutionary scenarios for galaxies. This 
will be the objective of Section [6] 

A key feature of our diagram is that it views galaxies as 
fundamentally populating a space of two parameters, an- 
gular momentum and mass, which are nearly equivalent 
to the more observationally accessible properties of bulge 
fraction and luminosity. In this framework, galaxies can- 
not be fruitfully reduced to a one-di mensional family co n- 
trolled by a single parameter (e.g., iDisnev et al1l2008[ ). 

Our diagram may also be contrasted with another cur- 
rently fashionable way to understand galaxies: as color- 
magnitude sequences th at are generally re l ated to star 
forma tion histories (e.g., iBaldrv et al.ll2004l : iFaber et"al~l 
2007). These properties are loosely related to j*-M* 
space if star formation generally occurs in high-j* disks. 
However, our framework is less astronomical and more 
astrophysical in nature, and we expect it to provide novel 
insights to galaxy formation that are complementary to 
other classifications, and perhaps more fundamental. 

Another recently-introduced classification for galaxies 
is also base d loosely on specific an gular momentum con- 
cepts: \r ([Emsellem et al.l I2007T ). which measures the 
rotational dominance in the central regions (typically in- 
side ~ R e /2) and is similar to a v/a metric. Applied 
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Fig. 15. — Specific angular momentum relative to the best-fitted trend for spiral disks. In the left-hand panel, these residuals are plotted 
vs. Hubble stage. For clarity, small random offsets have been added in the horizontal direction for the early-type galaxies. In the right-hand 
panel, the residuals are plotted vs. bulge-to-total mass ratio. The curved line shows a sample model prediction (not a fit to the data; see 
text for details). There are strong systematic trends of the j* residuals with respect to both Hubble type and to bulge- fraction, and the 
relative smoothness of this trend (particularly for the E/S0s) suggests that bulge-fraction is the more fundamental driving parameter. 

ble that the vertical displacements of in this diagram 
are somehow more directly related to Hubble morphol- 
ogy than to B/T (although one should keep in mind that 
B/T is one of the main factors in the morphological clas- 
sifications, along with spiral arm winding and dumpi- 
ness). 

To consider this point more clearly, and to better see 
the relative trends in the data, we flatten the j^-M* re- 
lations into one dimension, dividing by the mean trend 
for the spiral disks and thus generating the quantity: 



to early- type galaxies, a ho st of interesting patte rns and 
correlations have emerg ed (jEmsellem et aJJl2011f >. How- 
ever, this metric in practice is not only very scale depen- 
dent, but also misses exactly those scales that are most 
important for measuring true, physical angular momen- 
tum (recall Figure |3]). In fact, we have seen evidence 
that jn and the central Xr are disjoint properties: the 
slow rotators (low-Afl, galaxies) do not appear to deviate 
from the j+-M+ trend for fast rotators. 

A final related diagram to mention is j*— v c , where 
v c is the circular velocity, tracing the dynamical mass 
of a galaxy within some characteristic radi us (e.g., 
iNavarro fc Steinmetzl feOOCt iKassin et alJ l2012f >. There 
are complications with using this parameter space, since 
for spiral galaxies both and v c are normally based on 
the same rotation-velocity measurements, which causes 
a built-in correlation. Unlike M*, v c is not a physical 
quantity subject to straightforward conservation laws. 
In addition, a critical point for our goal of analyzing all 
types of galaxies in a unified manner is that it is very 
hard to estimate v c for a large sample of early-types since 
they rarely host extended gas disks. Instead, extensive 
data are required from other tracers such as stellar kine- 
matics (as needed for estimation), as well as gruel- 
ing dynamical modeling which even with the state-of-thc 
art techniques can stil l leave considerable uncertainties 
(|de Lorenzi et a l. 2009). Similar problems apply to a j+- 
M v i r (virial mass) diagram, where the masses can be es- 
timated on ly on a statistical rather than on an individual 
basis (e.g.. iDutton fc van d en Bosch 20l3). 



Alogj* = logj* - log j mo d(M*), 



(12) 



Our 



5.4. Examining the residuals 
bulge-disk framework, although rather com- 



pelling, is not a unique explanation for the systematic 
trends in the left-hand panel of Figure [14] It is possi- 



where j mo d is given by Equation (TTT]) . We plot Alogj* 
versus the Hubble stage parameter Tjjubbic in Figure [15] 
(left-hand panel) . There is clearly a strong positive corre- 
lation between Tnubbic and the j^-M* residuals. Among 
the spirals, this trend is clearest when considering the Sa- 
Sab versus Sb-Sc galaxies. The Scd-Sm galaxies appear 
to continue the trend, but they inhabit the lowest-mass 
area of the j^-M* diagram, where the mean relation is 
not defined well enough to be certain of the residuals. 

The SOs break the smooth trend of Alogj* decreas- 
ing for smaller Tnubbic- Many of them appear to have 
comparable specific angular momenta to typical Sb-Sc 
galaxies, which was foreshadowed by the rotation scaling 
relations of Figure [T^l The implication is that lenticu- 
lars and spirals are overall dynamically similar, differing 
more in their finer morphological features which may be 
related to star formation activity. We can thus think of 
these lenticulars as faded spirals, or of the spirals as reju- 
venated lenticulars, although they differ in average B/T 
values, and mo re nuanced compariso ns will require anal- 
ysis of Y* (cf. iWilliams et al.lT201~0h . As for the subset 
of lenticulars with low A log , they may either be very 
close to face-on, or else belong to a different family of 
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objects that are related to the ellipticals. 

Returning to our original hypothesis that B/T is the 
key parameter affecting the j^-M* trends, we consider its 
correlation with the residuals Alogj*. Since we do not 
actually have bulge/disk decompositions for the early- 
type galaxies in our sample, we introduce a novel tech- 
nique that uses the degree of central rotational support 
as a rough proxy for B/T . The idea here is that the bulge 
is to a first approximation non-rotating, so any observed 
rotation is from the disk: objects with higher (v/a) im- 
ply higher disk fractions and lower B/T. Appendix ID] 
describes our methods for early- type B/T estimation in 
more detail. For the late-types, we already have B/T 
estimates based on decompositions in the literature, as 
discussed earlier. 

We show the results in the right-hand panel of Fig- 
ure [TH The residuals do correlate clearly with B/T, in a 
fairly smooth trend that is followed equally well by all of 
the galaxy types, and which contrasts with the Tjjubbie 
trend. We have marked a simple expectation for the B/T 
trend with the curved line, given the summation of Equa- 
tion (|9|), along with an arbitrarily assumed jb = 0.1 x jd- 
This model mimics the data remarkably well, although 
it should be remembered that the agreement is some- 
what built-in already, since correlated rotational prop- 
erties were used both to estimate B/T and to calculate 
3*- 

Recalling that we also had to make strong modeling as- 
sumptions for the spiral bulges when calculating j*, the 
better connection of the residuals to B/T rather than 
^Hubble should be considered preliminary. It is also dif- 
ficult to tell how much of the scatter in at fixed B/T 
is due to observational error, and how much is due to 
intrinsic variations, i.e., with bulges and/or disks not 
following perfectly standardized j*—M+ relations. Defini- 
tive resolution of these issues will require more detailed 
bulge-disk decompositions of all ty pes of galaxies, in- 
cluding spectroscopic information (cf . Cortesi et al.ll2011t 
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We would however like to advance the proposition that 
bulge fraction is the fundamental driving parameter be- 
hind ji, variations, and is responsible for many of the 
observed variations in galaxy properties (see discussion 
in previous subsection). Not only does this make sense 
from a physical standpoint, but the agreements between 
ellipticals and spiral bulges in Figure IT4l (right) . and be- 
tween model and data in Figure ll"5l (right) . provide pro- 
visional but strongly suggestive observational support. 
The radially-declining rotation-velocity profiles of galax- 
ies like NGC 821 and NGC 3377 in Figure [6] could also 
be naturally explained by central disk components em- 
bedded in non-rotating bulges. Furthermore, we will 
see from consideration of a cosmological context in Sec- 
tion [(TU that the distribution of is more naturally rec- 
onciled with distinct disk and spheroid subpopulations 
than with a simple continuum of galaxy j*. 

5.5. Histograms of stellar j residuals 

Before moving on to theoretical analyses, we construct 
one more representation of the data whose relevance will 
become particularly clear in the next section. We com- 
press the preceding j+-M+ information into a histogram 
of residuals from the spiral disk relation, showing the 
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Fig. 16. — Histogram of specific angular momentum relative to 
the mean observed trend for spiral disks. In two of the panels, 
curves show example lognormal distributions for comparison to the 
data. In the uppe r pa nel, the red, green, and blue histograms show 
data from Figure [T5l for spirals, lenticulars, and ellipticals, respec- 
tively. The middle panel shows the bulge and disk subcomponents 
of spiral galaxies, with red and blue histograms, respectively. The 
lower panel is a summation of the data from the upper panel, after 
renormalizing each galaxy sub-type type by its frequency in the 
nearby universe (see main text). The specific angular momentum 
does not appear to have a simple lognormal distribution, and may 
even be bimodal. 

results in Figure [16] (upper panel) . Here it is apparent 
that the spiral galaxy data comprise a roughly lognor- 
mal distribution in Aj*, with an rms dispersion of ~ 0.2 
dex. The ellipticals have a less well-defined distribution 
that partially overlaps the spirals but is offset ^0.5 dex 
lower, while the small sample of lenticulars spans almost 
the full range of residuals. 

In the middle panel of Figure [T|)J we look instead at 
the disk and bulge subcomponents of the spiral galax- 
ies, where we have also ovcrplottcd a Gaussian with a 
width of uiog^ = 0.17 dex for reference. Given the un- 
certainties and possible selection bias in our analysis, we 
consider the disks to be reasonably consistent with a log- 
normal distribution. 

The A log j* distribution for the spiral bulges resembles 
that of the ellipticals in the sense that both are system- 
atically offset to lower values, as we have previously seen. 
The bulges apparently extend to much lower A log 
than the ellipticals, but as discussed in Section 15.21 this 
is not a secure result, given the uncertainties in the bulge 
calculations. 

Returning to the overall results, we would like to know 
whether or not galaxies follow a bimodal distribution in 
Alogj* as the top panel of Figure Il6l suggests. The 
complication here is possible bias in the galaxy sample: 
if we were to study all bright galaxies in a volume-limited 
sample, the Alogj* distribution might look very differ- 
ent. To investigate this issue, we must re- weight the 
distribution of j+ in our sample by galaxy type. 
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The simplest approach is to renormalize by frequency 
or number density. We use the ATLAS 3D results that 
70%, 22%, and 8% of the galaxies in the nearby uni- 
verse are spirals, lenticulars, and ellipticals (over a stel- 
lar mass range simila r to our observational sample; 
Cappellari ct al. 2011a). The fractions in our sample are 
63%, 14%, and 23%, demonstrating a strong bias toward 
ellipticals at the expense of lenticulars. 

We plot the re-weighted results in the lower panel of 
Figure 1161 showing also for reference a lognormal curve 
with <J\ gj t = 0.27 dex (a width that will be motivated 
in Section I6.2[) . The total distribution of log residuals 
appears slightly non-Gaussian, with a tail extending to 
low values. This feature may not be significant if one 
allows for systematic uncertainties in the selection effects, 
but the skewness will become clearer when compared to 
theory in Section ROl 

An alternative scheme would be to re-weight by the 
stellar mass density of the different galaxy types. This 
would bring us closer to a total distribution function for 
stellar j in the universe, rather than a distribution of 
galaxies with given . It is beyond the scope of this pa- 
per to carry out such an exercise in detail, but the basic 
outcome is clear. The high end of the mass distribution is 
dominated by early- types (cf. lower panel of Figure [33]), 
which means that the mass weighting would enhance the 
contributions of these galaxies relative to number weight- 
ing. The universal distribution of would then appear 
more non- Gaussian than in the lower panel of Figure [TBI 

We therefore find evidence that the residuals of the 
specific angular momenta of galaxies from the mean re- 
lation are not simply lognormal. The best match to a 
lognormal model is provided by the disk components of 
spirals, while the bulges and the ellipticals may comprise 
a distinct second population^ Again, a natural interpre- 
tation of this finding is that all galaxies are composed 
of some combination of high- and low-j* material, which 
may be identified with disks and bulges, respectively. 

Some implications of these results for galaxy forma- 
tion in a modern cosmological context will be discussed 
in the next section. It should be remembered, however, 
that our empirical findings — of specific, strong correla- 
tions between galactic angular momentum, mass, mor- 
phology, and bulge fraction — stand on their own and 
must be explicable by any successful theory of galaxy 
formation, whether now or in the future. 

6. CONNECTING TO THEORY 

We are now ready to present a fresh theoretical way of 
looking at galaxies, using the j^-M* diagram, which was 
introduced in F83, and which may now be reinvigorated 
by populating it with observational data for galaxies of 
all types. Our general approach is to take a step back 
from galactic details, whether these be spiral arms and 
dust lanes in observations, or unresolved gas physics and 
star formation recipes in simulations, and return to some 
simple physical parameters and conservation rules that 
may provide robust constraints and insights to galaxy 
formation. 

9 [Hernandez ct al. (2007) used a large photometric survey to 
estimate j* indirectly, with results that are less accurate than those 
presented here, but which similarly imply a bimodal distribution 
for ellipticals and spirals. 



We have shown in Sections 15.21 and 15.41 that the spe- 
cific stellar angular momenta of observed galaxies fol- 
low remarkably tight correlations with their masses and 
bulge fractions. Such patterns in Nature demand theo- 
retical explanations, as they could be tracing fundamen- 
tal physical processes. Indeed, the j*-M+ relation for 
spiral galaxies is well known in some circles, and pro- 
vides a crucial benchmark for models of galaxy forma- 
tion. However, the correlation for elliptical galaxies (al- 
ready shown in a preliminary version by F83) is less well 
known and addressed with theoretical models. Our goal 
is to advance a general, physical framework for integrat- 
ing these observational constraints into models of galaxy 
formation and evolution. 

Our approach here is different from, and complemen- 
tary to, the active field of hydrodynamical simulations of 
galaxy formation. Although such simulations have made 
notable progress toward the ultimate goal of reproducing 
realistic galaxies, they still have a long way to go, with re- 
cent work highlighting large differences in the basic prop- 
erties of simulated galaxies, depending on what code, res- 
olution, and physical rec ipes are used ( Scanna pieco et al~l 
12011 iTorrev et al.|[20lH ). 

Historically, such methods missed reproducing ob- 
served j* trends by factors of up to ~ 30, and even the 
most recent work shows variations at the factor of ~ 2 
level. The general concern is that many of the large- 
scale properties of galaxies could well depend strongly 
on transport processes at the scales of molecular clouds, 
which are not yet modeled satisfactorily in cosmological 
simulations. Therefore some caution is still needed in as- 
suming that the simulations are providing an adequate 
representation of reality. 

In this context, simplified "toy" models continue to 
play a key role in defining the broad but solid outlines 
of the galaxy formation theory that is required to match 
the observational constraints. These models may also 
prove useful in physical understanding of the output of 
numerical hydrodynamical simulations. 

We frame our analysis in the context of the current 
standard cosmological model for structure formation: 
cold dark matter wit h a cosmological constant (ACDM; 
iKomatsu et al.|[20Tl|) . This model makes specific, robust 
predictions for the angular momenta of DM halos. Be- 
cause the visible galaxies, consisting of stars and gas, 
are presumed to reside in these DM halos, we may then 
ask whether or not the observed stellar angular momenta 
bear any resemblance to the predictions for DM halos. 

We begin with the properties of ACDM halos as our 
"initial conditions" for galaxy formation, which we map 
to our observable space: j±-M+ for the stellar compo- 
nents of galaxies. We do this by parameterizing the re- 
tention of mass and angular momentum during galaxy 
formation, and then by introducing a menu of j^—M* 
vectors of change that correspond to plausible physical 
processes (outflows, mergers, etc.). 

We emphasize that the primary aim of this paper is not 
to concoct a new theory of galaxy formation, nor to weigh 
in on competing models by vetting specific simulation 
outputs against the j*-M* diagram. Instead, we wish to 
lay out a generalized framework that can both constrain 
and explain the models. The methodology and merits of 
this approach should become clearer as we develop the 
ideas throughout this section, and as we eventually work 
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through some practical examples. 

We develop general theoretical prediction s an d make 
basic inferences about j retention in Section 16.11 In Sec- 
tion [6T2] we investigate two possible explanations for the 
observed dichotomy between spirals and ellipticals. In 
Section[673]we consider coupling between changes in mass 
and angular momentum, and connect these to evolution- 
ary scenarios for galaxies. 

6.1. Basic constraints 

The overdense regions in an expanding universe are not 
spherically symmetric and exert tidal torques on each 
other, inducing a net angu lar momentum in each col- 
lapsing galaxy (Hoylc 1951). This rotational behavior is 
usually specified in terms of a dimensionless spin parame- 
ter that quantifies the dynamical importance of rotation, 
and is a combination of fundamental physical quantities: 



A = 



GM 5 / 2 ' 



(13) 



where J is the angular momentum, E is the energy (ki- 
netic and potential ), G is the g ravitational constant, 
and M is the mass |PeebIei[i969)E3 Whether analyzed 
through linear tidal torque theory, or through JV-body 
simulations of galaxy assembly, A is predicted to fol- 
low an almost lognormal distribution that is relatively 
insensitive to cosmological para meters, time, galaxy 
mass , and environment (e.g.. iBarnes fc Efstathioul 
1987b iZurek et al.lll98l iStemmetz fc Bartehnannl 1199a 
Cole fc LacevHl996t IMaccio et al.ll20r)7fi ~ 

The spin parameter provides a convenient way to char- 
acterize DM halos, but it is not straightforward to con- 
nect A to baryonic galaxies because it is not a physically 
conserved quantity (as energy is dissipated). We instead 
conduct our theoretical analysis in terms of the specific 
angular momentum parameter j, as we have done with 
the observations. Along with the mass M , j is a quantity 
that is potentially conserved at some approximate level 
during the evolutionary history of a galaxy. 

To re-cast A to j , we adopt a ACDM-based sphe rically- 
symmetric halo profile from iNavarro et all ((199 6) . trun- 
cated at the virial radius^ We then obtain: 



= 4.23 x 10 4 A 



\Q 12 Mr. 



2/3 



km s" 1 kpc. (14) 



We adopt a characteristic valueP^l of (A) = 0.035, along 
with a 1 a log dispersion of 0.23 dex, based on a study of 
relaxed halos in a cosmological simulation with WMAP5 

10 Recall that the parameters (J, E, M) can be translated 
roughly into a more observationally oriented basis set of rotation 
velocity, effective radius, and luminosity (v TO t, R e , L), where in ap- 
proximate terms: M oc L, E oc L 2 R^ , and J oc v I0 tLR,,. 

11 The virial radius is defined as bounding a region inside which 
the mean halo density is a factor of A v ; r times the critical den- 
sity p cr it = 3H 2 /(8nG). We adopt a WMAP5 cosmology, with 
H = 72 km s _1 Mpc -1 and A vir = 95.3 at 2 = IMaccio et al.l 
2008). T o calcu late E for this halo, we use an expression from 
Mo et aLI 11998) with a fixed concentration of c v ; r = 9.7; and we 



ignore variations due to concentration which affect A at the ~ 5% 
level. A related s pin-proxy parameter , A', is based on a singular 
isothermal sphere (Bullock et al. 2001), and is ~ 11% smaller than 
A. 

12 This is based on the average value of log A, but throughout 
this paper we use shorthand such as (A) and (j) for log-averages. 



parameters, by IMaccio et al.l ([2008). The log-averaged 
numerical coefficient in Equation (| 14[) then becomes 
1460 km s _1 kpc. Other recent studies are generally 
consistent with these results at the level of ~ 10%. The 
a = 2/3 expo nent is also an explic i t prediction of tidal 
torqu e theory (Sh ava fc Tullvl[l984l : IHeavens fc Peacockl 
|1988|) , and provides a reasonable approximation to the 
trends from direct calculations of .7 v i r and M v j r in A^-body 
simulations (jAntonuccio-Delogu et al.|[2010D . 

Equation (fT4"f can be considered as setting firm "ini- 
tial conditions" for galaxies, characterizing their angular 
momenta near the time of virialization. This is shown 
schematically in panel (a) of Figure 1171 which we have 
populated with toy-model "galaxies" consisting of pri- 
mordial halos of gas and DM. Their masses are drawn 
from a uniform logarithmic distribution, and their angu- 
lar momenta from a lognormal distribution using (j v ir) 
and cTiogj vir as above. 

We next consider a series of idealized evolutionary 
steps that allow us to parameterize evolution in the j-M 
diagram. We assume that the baryons consist initially of 
gas that is well mixed with the dark matter of its parent 
halo, and that does not collapse within the halo until 
after the linear and translinear regimes of tidal torque 
when most of the angular momentum is acquired. The 
gas may then be assumed to have the same value of j as 
the halo, which we show in panel (b) as a simple shift of 
the points to the left, according to a cosmol ogical baryon 
fraction of f b = 0.17 (jKomatsu et al.ll2011l ). 

In panel (c) we show what happens in a simple case 
where a fraction of the baryons form into stars, with 
a particular value of (/*) = 0.1, and a dispersion of 
o'log/* = 0.15 dex. Again, j is assumed to be conserved, 
and the galaxies shift to the left. It is also usually as- 
sumed, though not required by the diagram, that this 
process involves the formation of a thin stellar disk whose 
collapse was halted by the balance between gravity and 
centrifugal force. 

Our analysis does however assume that the baryon col- 
lapse extends all the way out to the halo virial radius. 
This conventional assumption is at some level implausi- 
ble since DM collapse and gas cooling are governed by 
different physical scales in space and time. A more gen- 
eralized approach where the baryon collapse radius is al- 
lowed to vary will be considered in Section [6.3.21 

Note that the /* parameter can take on a more gen- 
eral meaning of net stellar mass fraction relative to ini- 
tial gas mass, which allows for stars that are accreted 
by or ejected from the galaxy. We will shortly discuss a 
more refined model where /* varies systematically with 
mass, but for now we continue with our very simplified 
constant-/* model in order to consider its basic implica- 
tions. 

Our next model ingredient is an idealized process of 
angular momentum loss, with no concomitant change in 
mass, which we quantify by a fractional j net retention 
factor of fj. An example of such a process would be 
internal j transfer from the stars to the DM halo. Given 
the parameters and fj , we may then translate the j— 
M relation (fl"4|) for DM halos to an equivalent one for 
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Fig. 17. — Schematic evolution of galaxies in the space of specific angular momentum and mass. Each point shows a galaxy randomly 
selected from a simple model (see main text). Panel (a) shows the initial galactic halos of gas and DM. Panel (b) shows the gas component 
only, adopting a baryon fraction of /b = 0.17, with an arrow illustrating the direction that a single galaxy takes in this diagram. Panel (c) 
shows the stellar component after forming from the gas with an average relative fraction of (/*} = 0.1. Panels (d) and (e) show the stars of 
spiral and elliptical galaxies, respectively, after adopting more realistic variations of (f*) with mass. Panel (f) shows the effect of angular 



momentum loss, with a factor of {fj) =0.1. 



e.g., spiral galaxies probably do not exist at masses of M 

the stellar components of galaxies: 
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Note that these are simple, idealized models, and not every aspect should be taken literally; 
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ulation of galaxies can be determined by observations as 
a strict constraint on theory. 

We can immediately use Equation (|16[) in combination 
with the observational results for jo from Table [2] for 



where again using the prediction for (A), the numerical 
coefficient for (_?'*) becomes 1010 km s _1 kpc. 

This relation is identical to our parameterized fit to 
the observational data with Equation (fTTj) . modulo the 
numerical factors and the value for the exponent a. Since 
the observed j+-M+ relation can be approximated with 
a = 2/3 and a normalization jo, then we can express 
the difference between observation and theory through a 
simple combination of the parameters fj and 

jo = 1010 (fj f- 2/3 ) kms- 1 kpc. (16) 

Equations (p~4|) — (p~6l) are simple but powerful, allow- 
ing us to connect the visible properties of galaxies to 
their invisible DM halos, using some simple parameters 
and assumptions. They also provide robust observational 
constraints on some basic characteristics of galaxy forma- 
tion that are still far beyond the ability of raw theory to 

predict reliably. The average value of fj /* 2 ^ 3 for a pop- 



(15) fixed a = 2/3. We find that (fj fl 



2/3\ 



1.9 for Sb-Sm 



spirals and ~ 0.5 for ellipticals. For example, if we as- 
sumed an arbitrary (/*) = 0.2 for both types of galaxies, 
then we would infer (fj) ~ 0.65 for spirals and ~ 0.1 for 
ellipticals. This means a systematic difference in net an- 
gular momentum retention between the two galaxy types 
which, although there are many further details to work 
through below, will hold up as a basic result of this pa- 
per. 

To derive firmer constraints on fj, we need to break 
the f*—fj degeneracy by introducing well-motivated val- 
ues for for both spirals and ellipticals. We also need 
to consider the complication that /* cannot in reality 
have a simple, constant value, even on average. This 
is because the observed luminosity function of galax- 
ies has a dramatically different shape from the pre- 
dicte d mass function of DM halos (e.g.. iWhite fc Reesl 
1971 IWhite fc Frenkl fl99lt IMarinoni fc Hudson! I2002T 



Yang et al.H2003l; IMoster et aij|2010D . Below the charac- 
teristic "knee" luminosity L* , the galaxies are observed 
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to follow a shallower slope than the DM mass function 
dN/dM oc M -2 , while at higher luminosities, the obser- 
vations are steeper than the predictions. The implication 
is that the fraction of luminous-to-dark matter declines 
rapidly for galaxies fainter and brighter than L*; i.e., 
assuming a constant /b, the function (/*)(M v ; r ) has a 
characteristic inverted U shape. 

This empirical trend is thought to be caused physically 
by various feedback effects that inhibit star formation 
and become increasingly important in the low- and high- 
mass regimes (such as stella r and supermassive black hole 
feedb a ck, respectively; e.g..lLacev et al.lll993MCole et all 
1991 iSomeryille fc Primackl 119991: iBower et all 120061: 



Croton et alll2006[ ). Regardless of the explanation, any 



self-consistent ACDM-based model must incorporate a 
strong, systematic mass dependence on star formation 
efficiency, (/*)(M vir ). 

One might be concerned that such a mass dependence 
would transform an underlying j oc M 2 / 3 relation for DM 
halos into something very different for the stellar com- 
ponents of galaxies, and quite unlike our observational 
results. To check this, we will modify our simple model 
above to allow for a varying function (/*)(M v ; r ). Since 
this function is a tracer of undetermined baryonic physics 
during galaxy evolution, there is not yet any robust the- 
oretical prediction for it, but fortunately it can be es- 
timated empirically. This is done in an average sense 
through various techniques such as weak gravitational 
lensing, stacked satellite kinematics, and matching up 
the mass and luminosity functions mentioned above. 

There have been many studies that estimated 
(/*)(M v i r ), but few that did so separately for different 
galaxy types, which is important for our analysis. We 
therefore adopt the relations for (/*)(Mt) derived by 
iDutton et al.l (|2010h . For the spiral galaxies, we use their 
relation for "late- type" galaxies: 



</*>(M*) 



/o(M,/M ) 1/2 
[1 + (M*/M )] 1/2 ' 



(17) 



Below a characteristic mass log (Mq/Mq) ~ 10.8, this re- 

1 11 

lation has a dependence (/*) oc M* . At higher masses, 
it approaches a constant, /q ~ 0.33. Here we have con- 
verted the Dutton et al. results to our definition of the 
virial mass and to our adopted stellar IMF, while using 
ft = 0.72. 

For elliptical galaxies, we adopt the Dutton et al. rela- 
tion for "early-type" galaxies Fl 



</*)(M*) = 



f (M*/M ) 



0.15 



i + (M*/M y 



1/2 ' 



(18) 



where log (Mq/Mq) ~ 11.2, f ~ 0.14, and the asymp- 
totic behaviors at low and high masses are (/*) ~ Mj 315 
and (/*) ~ M~ 85 , respectively. One of the key features 
to notice here is that an elliptical galaxy typically has 



13 There has been very little work along these lines for el- 
liptical and lenticular galaxies separately, but there is some re- 
cent evidence that the halo masses for these types are the same 
(Truiillo-Gomcz et al. 2011). Note also that the Dutton et al. rela- 
tions were derived for somewhat smaller mass ranges than covered 
by our data, and that their stellar mass determinations may not 
be fully consistent with our methods. 



a much lower value of /* than a spiral with the same 
stellar mass: i.e., ellipticals inhabit systematically more 
massive DM halos, which in many cases extend up to 
"group" masses of My j r ~ 10 13 M Q and beyond (see also 
Ivan Uitert et alll2011[ ). 

These (/*)(M*) relations can be uniquely transformed 
to (/*)(Mvir), and taken together define an inverted U- 
shaped trend as discussed above. The relations were 
constructed using a compilation of different literature 
results, which showed an encouraging degree of mutual 
consistency, so we conclude that the average trends above 
are probably reliable at the ~ 50% (~ 0.2 dex) level 
There may also be non-zero galaxy-to-galaxy variations 
in /* at a fixed mass and type; the value of this scatter 
is less well established, but recent analyses s uggest that 
it ma y be ~ 0.15 dex (jBehroozi et al.l 120101: IMore et all 
I20T1 . We adopt this as our default value, which for- 
tunately is smaller than the expected dispersion in halo 
spin of ~ 0.23 dex and so will not have much impact on 
our conclusions. 

Using these variable (/*)(M v i r ) relations to construct 
mock j*-M* data sets as before, we plot the results in 
panels (d) and (e) of Figure [T7] For both spirals and 
ellipticals, we can see that the curvature in (/*)(M v j r ) 
translates to systematic deviations in the j + -M* relation 
from a simple a = 2/3 power-law. We will investigate 
how these deviations compare to real observations in the 
next subsection. 

Panels (d) and (e) of Figure [T7] also demonstrate that 
at masses of M* > 10 11 M@, the ellipticals are predicted 
to have higher j* than the spirals of the same mass, ow- 
ing to their differences in The more massive DM 
halos of ellipticals ought to provide larger virial-radius 
lever arms that lead to larger j v j r , and therefore larger 
j+ — j/they retain as much fractional angular momentum 
as spiral galaxies do. Therefore the observed offset in 
M* between spirals and ellipticals implies an even larger 
difference in (fj) than in the simple example above with 
fixed (/*) = 0.2. We will examine this apparent fj di- 
chotomy further in the next subsection. 

As a final illustrative exercise, we generate a mock 
data set for elliptical galaxies as in panel (e), then adopt 
(fj) = 0.1, with an assumed dispersion of ciog/* = 
0.15 dex. The results are plotted in panel (f), where 
we see that the galaxies have coincident ally returned to 
nearly the original j— M sequence for halos, modulo a 
little curvature and increased scatter. 

Figure [T7] thus shows how one could map the observed 
j^-M* properties of a population of galaxies (panel f) to 
a theoretical prediction for their halos (panel a), and re- 
cover some basic parameters describing galaxy formation 
(see Equation (jT6"]) l This formulation is closely related 
to a classic theoretical framework for the formation of 
spiral galaxy disks, whose observed sizes and rotation 
velocities are generally consistent with the approximate 
conservation o f primordial specific angular momentum 



(/,- ~ 1; e.g., 



Fall & Efstathiou 1980; Dalca nton et all 



ll^lMoetall IT9981 However, our formulation is more 
general by including also the early-type galaxies, as well 
as the bulge components within spiral galaxies (which we 
will discuss below). 

6.2. Investigating the spread in j+ 
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As just discussed, the observed dichotomy between the 
j+-M+ relations of spirals and ellipticals may imply dif- 
ferences in their specific angular momentum retention, 
expressed here by the factor fj. This interpretation is 
based on an implicit assumption that the parent halos 
of both galaxy types had the same average A. However, 
a natural halo-to-halo scatter in A is expected, and one 
could instead imagine the other extreme case, in which 
fj is the same for the two galaxy types, whil e their halo 
A val ues are systematically different (e.g., iKashlinskvl 
119821 iBlumenthal et al.lll984t iCatelan fc Theunsl 11996). 
In other words, spirals and ellipticals are drawn from 
the high- and low-spin tails of the A distribution, respec- 
tively. 

We call these two alternatives the "variable //' and 
"spin bias" scenarios. In reality, a mixture of both sce- 
narios may be present, which would be difficult to dis- 
entangle, but we can begin by investigating these two 
limiting cases in detail. Thus the aim of this section is to 
test how consistent each of these cases is with the data. 

The reason we can make headway on this issue is that 
there are predictions from ACDM not only for the aver- 
age value of A, but also for its probability distribution, 
i.e., a lognormal with a characteristic dispersion as dis- 
cussed in Section We continue to focus on the spirals 
and ellipticals as the two interesting extremes of the ob- 
served ji, range (at fixed Af*), and consider the lenticulars 
as intermediate either in fj or in A. 

We begin with the spin-bias scenario. If correct, adopt- 
ing a constant fj value for a complete, unbiased galaxy 
sample would allow us to work backwards to infer the un- 
derlying A distribution, which could then be compared to 
the theoretical prediction. One might think that we have 
already implicitly carried out this test by examining the 
residuals from the observed j+-M* relation in Section \h. 5 1 
and Figure 1161 However, that analysis did not account 
for the differences in /* between different galaxy types. 

We therefore proceed with a more direct comparison to 
theory by generating j^—M* model predictions for each 
galaxy type, and calculating the observed residuals with 
respect to these models. We use Equation (TT5]) with 
A = (A) = 0.035, along with the empirical (f ir )(M i ,) rela- 
tions (|T71) and (JTHJ) , and an ad-hoc (fj) — 0.55, to predict 
a mean j+-M+ relation for each galaxy type. We then 
derive the residuals Alogj* by subtracting the model 
from the observations as in Equation (|12l) . If the spin 
bias scenario is correct, then the properly reweighted dis- 
tribution of these residuals ought to follow a lognormal 
with dispersion eri og:)jl ~ 0.27 (which accounts for obser- 
vational errors and the intrinsic scatter in /*). 

Figure [18] presents histograms of these residuals, both 
by separate galaxy types (top panel), and in combina- 
tion (bottom panel), which uses a renormalization by 
frequency of g alaxy types from the ATLAS 3D survey, as 
in Section [531 We find that overall, the total distribution 
of Aj* has approximately the predicted width. However, 
the distribution in detail appears significantly different 
from a lognormal: there is an excess of low-Aj* galaxies, 
and a missing tail at high-Aj*. In particular, there are 
too many elliptical galaxies in the nea rby universe to be 
explained by the tail of low-spin halos Pi 
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Fig. 18. — Distributions of residuals in the observed stellar spe- 
cific angular momentum, with respect to the mean theoretical pre- 
diction for ACDM halos, after ass uming a fixed j'-retention param- 
eter, fj = 0.55. As in Figure [TBI red, green, and blue histograms 
in the top panels show the residuals for elliptical, lenticular, and 
spiral galaxies, respectively. The bottom panel shows the same 
distribution, renormalized for the relative frequencies of galaxies 
in the nearby universe. The curve shows a predicted lognormal 
distribution for comparison. The distribution of residuals for spi- 
ral galaxies is narrower than expected from the distribution of halo 
spins, while the overall galaxy distribution shows clear departures 
from the lognormal model (with an excess at low and a deficit 
at high >). 

This histogram analysis appears to exclude a simple 
spin-bias scenario, but there are some caveats, such as 
small sample sizes and the assumption of perfect lognor- 
mality for the distribution of halo spins. We can make 
further progress by recognizing that the scenario makes 
predictions for the residuals not only for all galaxies 
combined, but also as a function of mass. This is be- 
cause A is not predicted to depend on halo mass, while 
the relative frequencies of different galaxy types are ob- 
served to vary strongly. One can then immediately sec 
a serious problem with the spin-bias scenario: at high 
masses, almost all of the galaxies are ellipticals, which 
should thus be an unbias ed population representing the 
full range of halo spins (jDutton fc van den Boschl 120121 
made a similar point for low-mass disk galaxies). 

We investigate this issue in more detail by construct- 
ing a mock data set as in Figure [TT1 while this time in- 
corporating a schematic model for spin bias. We now 
assume that all galaxies have fj = 0.45, with the late- 
types inhabiting the high-spin halos, and the early-types 
the low-spin ones. Using the number densities of early- 
and late- types as a function of M* from ATLAS 3D , we 
use the (/*) (M+) relations to translate this to the relative 
fractions at fixed halo mass (which can be quite different 
from the fractions at fixed M*). We then randomly draw 



14 [Hernandez et al. (2007) also found in attempting to infer halo 
A values for spirals and ellipticals that an ad-hoc rescaling of the 



elliptical values was required in order to avoid a double-peaked A 
distribution. 
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Fig. 19. — Stellar specific angular momentum vs. stellar mass, comparing mock data generated from ACDM-based models (left- and 
right-hand panels) to real data (middle panel). The model on the left includes halo spin bias, while the model on the right assumes 
systematic differences in angular momentum retention between spirals and ellipticals. Blue open squares and red filled circles show spirals 
and ellipticals, respectively, with the solid blue and dotted red lines showing the best-fit power-laws for the real data. The relation for halos 
is also shown for reference as a gray dot-dashed line. The mock data sets include intrinsic scatter in the parameters A and /* at a given 
mass, but not observational errors. The simple variable- fj mock data on the right resemble the real data, while the spin-biased model does 
not. 



a distribution of biased spin parameters for each galaxy 
type; e.g., if spirals comprise 25% of galaxies at a given 
mass, we draw mock spirals from the top quarter of the 
spin distribution. We also adopt a similar mass range 
and total number of galaxies as in our real data sets. 

We show the resulting j+-M* mock data set in the left- 
hand panel of Figure [T5l which can be compared to the 
real data in the middle panel. We see that the low-mass 
ellipticals could indeed be drawn from only the low-spin 
tail because of their rarity. However, at high masses the 
ellipticals are common and their predicted values are 
similar to the spirals. To salvage the spin-bias scenario 
would thus seem to require a mass-dependent bias, which 
seems epicyclic and therefore not appealing^ 

The biasing idea can also be discredited by environ- 
mental considerations: there are strong observational 
correlations between environmental density and galaxy 
morphology, but as mentioned earlier, halo spins in 
theory depend only weakly on environment (which has 
some observational suppor t in the case of dis k galaxies; 
iCervantes-Sodi et all 120081 : iBerta etafl 120081 ) . In addi- 
tion, if we consider disks and bulges to be manifestations 
of the same j+—M+ trends as spiral and elliptical galax- 
ies, then the coexistence of these subcomponents within 
the same galaxies provides a clear argument against halo 
spin bias. 

We next turn to the variable- fj scenario, where spirals 
and ellipticals are drawn from the same underlying distri- 
bution of halo spins, but their baryonic components have 
systematic differences in retaining j. Given that we know 
(/*) for each galaxy type, Equation (IT6l) suggests that we 
can immediately use the observed jo normalization to in- 
fer (fj). However, the situation is more complicated since 
(/*) varies with mass and therefore one does not expect 



15 There may be reasons of stability for ellipti c als to be dom- 
inant at high masses (e.g., Dalcanton ct al. 1997; van de n Bosch] 
1998; Dutton & van de n Bosch 20121). but this ostensibly changes 
the morphology and not j*. 



an exact a = 2/3 for fixed fj (recall Figure UTTd.e)). 

As we did for the spin-bias scenario, we again construct 
mean j+-M+ relations for each galaxy type, while now 
leaving fj as a free parameter. Carrying out least-square 
fits to the data, we find values of (fj) = 0.56 ± 0.03 and 
(fj) = 0.12 ± 0.01 for the spiral and elliptical galaxies, 
respectively. The difference in (fj) of a factor of 4.7 ± 
0.8 is slightly larger than the observed j*-M* relative 
offset, as anticipated in the previous section because of 
the differences in (/*) (e.g., Equation (JTHJ) ) 

The next step is to verify that these best-fit models 
provide reasonable representations of the data. We again 
construct mock data sets, using the new fj models (with 
0.15 dex of scatter in /*), and show the results in the 
right-hand panel of Figure fT9l Here we see that, unlike 
the spin bias model, these variable- fj models provide a 
remarkably good match to the data. The curvature of the 
predicted j'*— M* relation turns out to be imperceptible, 
once we account for observational errors, small- number 
statistics, and a limited mass rangeFI Furthermore, the 
observed slope for the spirals is shallower than a = 2/3, 
which is predicted by the model. 

This comparison does not entirely succeed in account- 
ing for the scatter about the j+-M+ relations. As can 
be seen in Figure [T9l the real observations appear to fol- 
low tighter trends than predicted by our simple model, 
for both spirals and ellipticals. The model fits give rms 
scatters of a\ os f. = 0.18 dex and 0.25 dex for the spirals 
and ellipticals, which is already less than the expected 

16 Given the degeneracy between fj and /*, in principle the 
inferred fj dichotomy could be an artifact of errors in our adopted 
values for (/*). However, these errors would have to amount to a 
combined factor of ~ 5: e.g., with true (/*) ~ 0.1 for the spirals 
along with ~ 0.2 for the ellipticals, rather than ~ 0.25 and ~ 0.1. 

17 Future empirical estimates of and M* over a larger dy- 
namic range could provide a strong test of constant- fj scenarios. 
Given the observational difficulty of measuring at high masses 
where the underlying halos pertain to entire galaxy groups and 
clusters, the best prospect for improvement would be to study 
lower-mass galaxies, with log(M*/MQ) < 9. 
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scatter of 0.27 dex from A and /*, even without allow- 
ing for measurement errors, and scatter in fj (see also 
the histogram of spirals in the top panel of Figure [TBI 
compared to the curve in the lower panel). 

One possible explanation for this reduced scatter is 
that the baryonic processes responsible for j-loss could 
act as some kind of "attra ctor" to specific values of fj 
(cf. Ide Jong fc Lacevil2000l) . Alternatively, halo spin bias 
could be at work in a secondary role, even while fj vari- 
ation is the primary effect 

Our overall conclusion is that the variable- fj model re- 
produces the j^-M* observations well in general, is fairly 
insensitive to the exact trend of (/*) with mass, and does 
not require any additional variation of (fj) with mass. 
The spirals appear to have been fairly efficient in pre- 
serving the specific angular momentum imprints of their 
parent halos, while ellipticals have lost the vast majority 
of theirs. 

This is a plausible scenario from a physical standpoint 
if we return to our proposed framework where all galax- 
ies are composed of bulges and disks (Figure [2] and Sec- 
tion [53]). Unfortunately, we do not have (/*)(M+) re- 
lations for the bulges and disks themselves in order to 
directly derive their (fj) trends. However, given the sim- 
ilarities in j+-M+ that we found between these subcom- 
ponents and the galaxies overall, it seems reasonable to 
suppose that bulges and disks have (fj) ~ O.f and ~ 0.6, 
respectively, and that these values are characteristic of 
two distinct modes of galaxy evolutionpj We will return 
to this topic in the next section. 

Our conclusions about spiral galaxies echo 
similar findings in the literature, which have 
typically inferred (/,) ~ 0.5-0.6 overall (e.g 



Nav arro fe Steinmeta 12 000; Du tton et al.l 
20091: iDutton fc van den Boschl 120121 : 



2007; Burkert 
iKassin et al 



20121 ). In particular. IDutton fc van den Boschl $201% 



used a model parametrization similar to our 
and found that (fj) is fairly constant over a wide mass 
range. Note that these authors used a parametrized 
mass model to fit the Tully-Fisher relation, which was 
then converted to an average j'vir - -^vir relation. Our 
approach works instead in the space of observables, 
j+-M+, which is more direct and transparent while also 
allowing us to analyze galaxy-to-galaxy variations PI 

Our finding for the ellipticals is novel, as neither the 
predictions for j^-M* of ellipticals, nor their subsequent 
fj inferences, have been well-studied before now. We 

18 It has been suggested that later type g alaxies are biased to 
lower spin halos (D'Onehia & Burkert 2004). If correct, the net 
impact on the scatter is unclear, but one implication is that the 
fj dichotomy between spirals and ellipticals would be even larger 
than in our no-bias scenario. 

19 One concern here is that for more bulge-dominated galaxies, 
one might expect the disk-only {/*) to be relatively low, and thus 
the disk j* to appear relatively high. However, the observations 
are somewhat suggestive of the opposite trend, i.e., disk j* anti- 
correlating with B/T. 

20 As a consistency check, we also take a slightly different ap- 
proach and make a model prediction for the m ean relation be- 
tween size and rotation v elocity for spirals (cf. IMo et ahl 119981 ; 
IBurkert fc D'Onghiall200l 1. We adopt a value of (/*) = 0.56, and 
rather than assuming some function (f i ,)(M it ), we relate the disk 
rotation and the virial circular velocity by v s ~ 1.2v v i r . Given 
(A) = 0.035, there is a linear relation predicted b etwe en v s and 
a e , which we show in the right-hand panel of Figure 1121 To zeroth 
order, this prediction agrees well with the spiral data. 



have not carried out a comparable analysis on lenticu- 
lars since the constraints on them are less certain. Qual- 
itatively speaking, their observed logj* normalization is 
between the other two galaxy types, which for plausible 
values of (/*) implies (fj) values that are intermediate 
to those for the spirals and ellipticals. In addition, there 
may be two subpopulations of lenticulars as discussed in 
Section 15^41 with low and high (fj). 

There are two interesting implications about these 
findings. One is that that we now have a remarkably sim- 
ple and successful framework for describing and connect- 
ing some of the most fundamental properties of galaxies. 
The observable galaxies may be connected to their un- 
observable host halos using j+ and M* along with some 
relatively basic parameters fj and /*. Such a model may 
appear implausibly oversimplified in the light of our ever- 
expanding awareness of the complexities of galaxy forma- 
tion physics, but for some reason it seems to work. 

The other implication is that these parameters may 
give us insight into the formation of disks and bulges, and 
into the origins of the Hubble sequence. To illustrate this 
point, we use our modeling procedures as described above 
to work backwards and estimate /* and fj values for 
individual galaxies. The outcome is shown in Figure 1201 
where one should focus on the average results for each 
galaxy type, since no attempt was made to model the 
scatter in /* and A. 

The general picture that we obtain is that spiral and 
elliptical galaxies are clumped around two regions of pa- 
rameter space: (f*,fj) ~ (0.25,0.55), and ~ (0.1,0.1), 
respectively. Whatever processes formed and shaped 
these galaxies were efficient at both forming stars and re- 
taining total specific angular momentum for the spirals, 
and inefficient for the ellipticals. 

As discussed in Section [TJ early cosmologically-based 
simulations struggled to reproduce such high fj val- 
ues for spirals, finding typically fj ~ 0.01-0.1, which 
was later realized to be due in part to numerical arti- 
facts, and in part to inadequate feedback recipes. Feed- 
back could be particularly important for slowing down 
gas collapse and star formation so that the baryons are 
not affec ted by torque-dr i ven j transfer during early 
mergers dWeil et alj Il998t iSommer-Larsen et all 120031: 
iHummels fc Brvanll2012tlScannapieco et al.ll2012f) . How- 
ever, whatever physical processes are now invoked to 
explain the fj values of spirals must simultaneously al- 
low for much lowe r /,■ in ellipticals (e.g., by having less 
efficie nt feedback; iZavala et al.ll2008t iScannapieco et all 
[2001 . 

6.3. Physically-motivated models for galaxy evolution 

Now that we have derived a comprehensive frame- 
work for connecting j+-M+ observations with simulated 
ACDM halos, and thereby derived generic constraints on 
specific angular momentum retention, fj (Figure [2D]) , we 
will work through some case studies of plausible physi- 
cal processes in galaxy formation and evolution. These 
cases are not meant to be exhaustive, nor to provide im- 
mediate ammunition for current debates about galaxy 
formation, but to serve as practical examples of how the 
j-M diagram can be used as a tool to furnish physical 
insight. The models involved will treat fj and /* as co- 
variant parameters, unlike in the previous sections where 
for simplicity they were independent. 
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Fig. 20. — Specific angular momentum retention fraction plotted 
against stellar mass fract ion , as inferred for individual galaxies, 
with symbols as in Figure [19] The dotted diagonal line is the one- 
to-one relation, and the gray double-arrow shows the direction of 

2/3 

the uncertainties as driven by the fj oc f t ' degeneracy. The width 
of the shaded region around fj = 1 corresponds to the scatter in 
spin expected for ACDM halos. The black arrows show s chem atic 
vectors from 1:1 and 1:10 mergers, as discussed in Section [6.3l The 
spiral and elliptical galaxies occupy distinct regions of the diagram, 
while a simple model implies that converting spirals into ellipticals 
would require a very large amount of growth through ~ 1:3 mergers. 

A general constraint to keep in mind is that for each 
galaxy type, fj is approximately constant as a function of 
mass, including little additional scatter, which accounts 
for the observed j-M relations appearing so similar to 
those for theoretical DM halos. Any model for angular 
momentum evolution should explain why galaxies appear 
to remember so faithfully the overall initial conditions of 
their parent halos. 

The challenge of this fj constancy has been recog- 
nized previously for disk galaxies. There are a variety of 
physical mechanisms during galaxy evolution that could 
involve j transfer (e.g., gas cooling and feedback), but 
unlike gravitational clustering, these baryonic processes 
(and the resulting fj values) are expected to depend 
strongly on mass, which appears to require some degree 
of fine-tuning to reconcile wit h the observations (e.g., 
IDutton fc van den Bosch 2 0121) . Our inclusion of early- 
type galaxies in this framework, with near-constant fj, 
deepens the mystery: there are now two fine-tuning con- 
spiracies to explain. 

Here we emphasize again a distinction from com- 
parisons between internal distributions with radius 
of j for stars and DM halos (e.g.. | Bullock et all 
200lt Ivan den Bosch et~aLl [20011: IMaller et all 120021: 
Sharma fc Steinmetzll2005ll As mentioned in Section [TJ 
there is ample reason to expect redistribution of to 
occur within the baryonic component of a galaxy and 
thereby violate strong j conservation. However, this 
does not affect our examination of weak conservation, 



where the overall value of j may remain roughly the same 
(assuming negligible transfer of j between baryons and 
DM). 

We may reduce the potential explanations for the sys- 
tematic difference in fj between spirals and ellipticals 
into two basic scenarios, which we will examine before 
summarizing the overall picture. One general scenario 
is an internal angular momentum bias, where high- and 
low-j* galaxies were formed from parts of their available 
gas supply that had preferentially high or low j. The 
other is that these galaxies experienced systematic dif- 
ferences in angular momentum transport after star for- 
mation, and during subsequent galaxy assembly phases. 

Below, Section 16.3.11 discusses outflow and stripping 
scenarios, Section 16.3.21 considers biased collapse, and 
Section 16.3.31 examines mergers. Section 16.3.41 surveys 
the plausibility of these evolutionary modes in the light 
of the j-k-M* observations. 

6.3.1. Outflows and stripping 

One example of the first scenario involves gas outflows, 
whether caused by galactic winds or by some other mech- 
anism, let us assume that the baryons in a galaxy col- 
lapse into a thin disk while preserving the total specific 
angular momentum, i.e., fj — 1 (recall Figure nTT b)). 
The local specific angular momentum within the disk, 
j s (R) oc Rv ro t(R), is assumed to increase monotonically 
with galactocentric radius, which is unavoidable if the gas 
follows co-rotating circular orbits (the rotation-velocity 
profile cannot decrease any more rapidly than Keplerian, 
while the lever arm R in the j calculation increases lin- 
early) . 

Before many stars form, an outflow begins which we 
parameterize by a mass loss that is proportional to the 
gas surface density to some unknown power ft: 

AM g oc E£. (19) 

Because the gas is presumed to settle into a configuration 
where the density increases toward the center (e.g., an 
exponential profile), the parameter ft translates into a 
biased removal of gas from different disk radii, which in 
turn means depletion of gas parcels with systematically 
different j g . 

To analyze this scenario further, we now introduce Fig- 
ure [2U which like Figure [T71 illustrates schematic vectors 
of mass and angular momentum evolution, but now ex- 
tends to more specific, physically-motivated processes. 
In Figure l2"TT a). the horizontal arrow to the left illus- 
trates an outflow with ft = 0: the gas everywhere in the 
disk is depleted by an equal fraction, and its initial spe- 
cific angular momentum is preserved, while its mass de- 
creases. If ft > 0, then the outflows occur preferentially 
in the high-density, central regions that have relatively 
low j g , and so the overall j g for the galaxy increases (di- 
agonal arrow towar d upper left; cf. iBinnev etHI I200H: 
iSharma et all [20121 ). If ft < 0, then the mass loss is 
preferentially from the outer regions, and the overall j g 
decreases (diagonal arrows toward lower left) . Thus, out- 
flows could in principle produce either a net increase or 
decrease in fj. 

It should be kept in mind that these outflows repre- 
sent only material that is launched completely out of the 
galaxy, never to return. Other types of outflows may 
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Fig. 21. — Schematic evolution of galaxies in specifical angular momentum and mass, as in Figure [TT] but now considering evolution 
through gas outflows, stripping, and biased baryon collapse, and galaxy mergers. Panel (a) show initial conditions for pre-collapse gas 
(dots), and possible evolutionary vectors from outflows and stripping (arrows; see text for details). Panel (b) shows the collapse of gas 
and formation of stars at some initial redshift z\, preserving the j*—M* values until a final redshift zq (black arrow to the left, with dots 
illustrating a population of galaxies). The halo grows until redshift zq with no further star formation (black arrow to upper right). At zo, the 
expected trend with perfect j conservation is the dotted line, and net values for /* and fj would be inferred using the leftward and downward 
gray arrows, respectively. Panel (c) shows initial conditions for DM halos as gray dots, and schematic vectors of evolution through mergers 
(gray arrows): mass growth (to the right), specific angular momentum decrease through cancellation of the spin components (downwards), 
and increase through the orbital component (upwards). The net evolution is a black diagonal arrow to the upper right. The upper dotted 
track marks the initial conditions for stellar disks, and the blue dots show disks after having undergone four 1:1 mergers each. The upper 
black curved vector illustrates the typical evolution of a galaxy, with each black dot marking the beginning of a discrete merger event. The 
lower black curved vector shows the same for a series of 1 : 10 mergers (note that for clarity, the curved vectors are arbitrarily shifted relative 
to the fj = 1 starting point for the DM vector). In both cases, after the mass has grown by a factor of ~ 2, the orbital dominates the 
evolution, moving merger remnants along a j'*— M± track parallel to, but lower than, the initial disk trend. 

away from the initial j—M sequence. This conclusion 
is supported by iV-body simulations of ACDM halos, 
which find that the environmental dependencies of halo A 
are fairly weak (IZurek et al.lll988t ILemson &; Kauffmanr] 
[19991 iReed et alJl2005f T 

If instead the stripping occurs after the gas collapse, 
then j and M decrease for the DM but not for the 
baryons. This leads to elevated values of fj and /*, which 
could be investigated through observational constraints 
on M v ; r for field galaxies in comparison to satellite galax- 
ies in massive groups. 



also occur, where gas is expelled outward but remains 
boun d and falls inward again, as in a galactic fountain 
(e.g., iBrook et al.l I2012I ). However, such internal pro- 
cesses might alter only the detailed distribution with ra- 
dius of j, and not affect the overall value which concerns 
us here (see the discussion above of weak and strong j 
conservation) . More complex scenarios could also be con- 
sidered, where fountain material interacts with halo gas 
and e xchanges angular mo mentum (e.g., IMelioli et al.l 
120091: IMarinacci et al.l 1201 ID , leading to shifts in for 
the stellar disk that eventually forms. 

A mechanism related to gas outflows is galaxy stripping 
through gravitational interactions with other galaxies in 
a dense environment. Here the effects on and M* de- 
pend on whether the tidal stripping occurs before or after 
the gas collapses. If a galac tic halo is tidally st ripped be- 
fore the gas collapses fe.g. JLarson et aT1ll980f) . then the 
reservoir of M g and j s available for collapse is depleted 
in a manner that depends on the internal distribution 
of these quantities. F83 adopted some plausible distri- 
butions and worked out the resulting j—M changes: we 
will not repeat the analysis here, but merely show the 
equivalent evolutionary vectors as the three arrows in 
Figure [Ufa) pointing downwards to the left. 

There are two key features to notice with the gaseous 
stripping arrows. One is that unlike outflows, this strip- 
ping can only decrease fj {(i < 0) since it acts solely 
on the outer regions. The second is that plausible j- 
loss vectors are accompanied by substantial mass loss, 
which means that it is fairly difficult to move galaxies 



6.3.2. Biased collapse 

There is another scenario that is functionally equiv- 
alent in the j-M diagram to outflow or stripping, but 
which merits special attention. Here we consider a 
spatially-biased subcomponent of the initial gas which 
collapses and forms stars. Rather than our default as- 
sumption of uniform efficiencies /* and fj throughout 
the virial region, we assume that stars form preferen- 
tially in the inner regions of the halo, while the outer 
regions remain largely gaseous and form relatively few 
stars. 

This scenario was introduced bv lFalll (|2002f ) and is mo- 
tivated by the higher densities, and thus overall gas dis- 
sipation rates (through cooling and cloud collisions), in 
the inner regions. The consequent spatial bias in star 
formation can also be understood as a temporal bias, 
if one considers an idealized onion-shell model wherein 
galaxies form by inside-out collapse, with virialization 
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and star formation occurring first in th e central regions 
(cf. Ivan den Boschl 1998; K epnerl [1999D . Even in more 
realistic, hierarchical galaxy models, it is uncontroversial 
that a large fraction of the baryons within a galaxy halo 
at any given time will not yet have formed stars, and are 
located preferentially at larger radii. The stars observed 
in a galaxy at z = will have formed on average at higher 
redshifts, and from gas that was more centrally confined 
than the z = virial volume. 

Because j for a ACDM halo is expected to increase 
systematically with both internal radius and time, the 
above biasing scenario implies that j+ for a galaxy will 
be lower than its total j (including DM). Su ch a biasing 
framework was used by iKassin et al.l (|2012fl to connect 
observed disk galaxies with simulated ACDM halos, and 
thereby infer a radius of baryonic collapse. Here we out- 
line a generic toy model of collapse bias, to understand 
its implications in the context of j-M evolution vectors. 

For simplicity, we adopt a step-function model where 
at an initial redshift z- u all of the gas within the virial 
radius instantaneously collapses and forms stars with 
perfect efficiency and angular momentum conservation 
(/* = fj = 1), and subsequently no star formation occurs 
(/* = 0). This scenario is illustrated by Figure l2"TT b). 
where z\ marks the initial halo parameters. The leftward 
arrow shows the formation of the stars, with j*-M + pa- 
rameters that are preserved until zq = 0. The diagonal 
arrow to the upper-right shows the subsequent evolution 
of the halo. Because the halo continues to grow in M 
and j, the net values of /* and fj for the stars will de- 
crease with time, which is illustrated by the gray arrows 
which are the inferences made by connecting the final 
conditions of the halo and stars. 

This biasing scenario might seem to provide a tidy al- 
ternative for understanding galaxies that have apparently 
experienced baryonic angular momentum loss. However, 
it is important to realize that such biasing cannot ex- 
plain just any arbitrary set of j'^-M* observations. For 
example, the vectors in Figure [2"TT b) were constructed to 
represent a typical early- type galaxy with a net /* = 0.1 
at z = 0, which turns out to have a net fj = 0.22, i.e., 
not reproducing the apparent (fj) ~ 0.1 from observa- 
tions. Note that this model had an initial /* = 1, but in 
reality, we expect an initial /+ < 1, which would increase 
the discrepancy. We will discuss this scenario further in 
Section 16.3.4) for now, it serves as an important illus- 
tration of how constructing physically-motivated vectors 
in the j^—M* diagram can provide tight constraints on 
possible evolutionary scenarios. 

6.3.3. Mergers 

We next consider galaxy merging following star for- 
mation, which is likely to be more important for ellip- 
ticals than for spirals. The mass of a galaxy increases 
through a merger, while its final j is determined by the 
vector sum of three initial j components (the internal j 
for the two progenitor galaxies, and their relative orbital 
j), as well as by any exchange of j with the environment 
(e.g., between the stars and their surrounding DM halos). 
The random relative orientations of the first two compo- 
nents will cause them to partially cancel out, which con- 
tributes a net decrease to j. That is, after N equal-mass 
mergers, there will be average trends for the remnant of 
J oc N 1 / 2 and M oc N, and therefore j oc TV" 1 / 2 (jFalll 



[19791 lAarseth fc Fall[r980l ). The orbital j and the j ex- 
change processes are more difficult to model a priori. 

The effects of mergers on DM halos have been stud- 
ied extensively through numerical simulations, result- 
ing in a general picture where major mergers tend to 
"spin up" the halos, while minor merge rs and smooth 
accretion tend to spin them down (e.g., iGardnerl 120011: 
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Mailer et all l2002t IVitvitska et all ._ _ 
2004t iD'Onghia fc Burkertl 120041: iHetznecker fc Burkert 
2006). Given that the jvir^Afvir relation is scale- free 



and has a normalization that is ex pected to change only 
gradu ally, if at all, with time (e.g.. iNavarro fc Steinmeta 
11991 . we conclude that for individual halos, the co- 
addition of the above processes must amount to a 
random-walk that takes them on average along the jvir~ 
M vu - sequence. 

We illustrate this process in Figure l2lT c) with a 
schematic evolutionary vector for galaxy halos, broken 
down into subcomponents of j v i T and M V - 1T changes q 
Doubling the mass should typically increase j v - lr by a fac- 
tor of 2 2 / 3 = 1.6. 

The effects of mergers on the stellar components of 
galaxies, which have collapsed by large factors within 
their DM halos, are somewhat different. Qualitatively 
speaking, it is a generic dynamical requirement that the 
stars shed some of their orbital angular momentum, via 
tidal torques or dynamical friction, i n order to coalesce 
into a bound merger remnant (e.g.. | Frenk et all 119851 : 
IZurek et al.lfl988t IBarneiU98l ID'Onghia et al.ll20"06ll 

More quantitatively, we may make an initial, plausible 
guess that the "final pass" of the merger before coales- 
cence involves an impact parameter and relative velocity 
that are similar to the stellar scalelength and circular 
velocity of the larger progenitor. This would mean that 
the smaller progenitor would bring in an orbital of a 
similar magnitude to internal i of the larger progenitor 
(i.e., AJ* = j*, 2 M*, 2 ~ M*' 2 ). 

We sketch out some implications of this kind of merger 
evolution in Figure |2"T1c). Starting with galaxy disks 
randomly selected along the median trend as in 

Figure [17tc) (adopting a simple /* = 0.1 model with 
scatter included for halo A), we apply a sequence of four 
mergers to each disk. Each merger has a 1:1 mass ratio, 
and the relative vectors of internal and orbital are 
selected random ly (this is similar in spirit to the orbital- 
merger model of Ma iler et al.ll2002D . The blue dots show 
the end result after the merger sequence, and the upper 
arrow shows the median trend for a single galaxy, with 
black dots marking the discrete merger events. Note that 
at this point, the series of four 1:1 events is meant as a 
thought experiment and not necessarily as a likely merger 
history. 

After an initial decrease of in the first merger from 
cancellation of the internal spin vectors, the orbital 
dominates the evolu t ion of the merger remnant (e.g., 
lAarseth fc Fal ITiM IHetznecker fc Burkertl [2001 this 

21 In the merging of DM halos, the resulting angular momen- 
tum and mass are not the simple sum of those properties from 
the progenitors. The combination of the two virial regions in a 
merger increases the density within a fixed physical radius, but 
also increases the volume of the virial region, so that more of the 
surrounding material falls under the gravitational sway of the two 
galaxies together. A 1:1 merger typically increases M v i r by a factor 
of ~ 2.3; similar effects apply to j vlr . 
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Fig. 22. — Distributions of specific angular momentum residuals, 
relative to th e m ean trend for spiral disks, using the same analysis 
as in Figure Etc). The right histo gram shows the disk initial 
conditions. The middle and left histograms show merger remnants 
after having grown by a factor of 16 in mass, for 1:1 and 1:10 
mergers, respectively. The j* distribution has a smaller mean and 
dispersion for the 1:10 mergers than for the 1:1 mergers. 

also means that the results hardly change if the "ac- 
creted" galaxies are low-j* spheroids rather than disks 
as we have assumed here). Because the orbital j* term 
is assumed to be similar to the disk j*-M* trend, the 
final trend for the merger remnants parallels the disk 
trend, while being offset to lower by a factor of ~ 2 
(~ —0.3 dex). Referring back to Figure fTTl this corre- 
sponds to an effective angular momentum loss term of 
fj ~ 0.5. The distribution of the offset is also shown by 
a histogram in Figure [22l 

We have carried out the same exercise for a series of 
1:10 mergers, with a median trend shown by the lower 
vector in Figure r2"TT c). The result is similar to the 1:1 
case, with orbital j* dominating the evolution after the 
galaxy grows in mass by a factor of ~ 2. However, the 
final j* trend is now lower than the disks by a factor of 
~ 6 (~ —0.8 dex; fj ~ 0.15), with less scatter than in 
the 1:1 case (see Figure |2"21 again) . These differences arise 
because there is less stochasticity with the 1:10 mergers, 
where random walk effects tend both to wash out varia- 
tions and to dilute the orbital contributions to j*E3 A 
more realistic mixture of multiple mergers with varying 
mass ratios would presumably produce a j+ distribution 
with a peak intermediate to our 1:1 and 1:10 scenarios, 
and with a larger scatter. 

22 This scenario has some parallels to discussions in the literature 
about the systematic relations between angular momentum and 
merger histories, and the implications for the observed properties 
of galaxies (e.g . . ID'Onehia fc Burkertl 12004; IVitvitska et alj|2002t 
IPrimacM [2005t ID'Onehia fc Navarrd I2007h Isharma et alj 120121) . 
However, those studies did not always make a clear distinction 
between the differing merger dynamics of DM halos and of their 
embedded stellar components. 



These calculations are laden with simplifying assump- 
tions and could easily be wrong by a factor of 2 in j±. 
However, they are meant to illustrate some possible im- 
plications of merger activity in a hierarchical context. 
First of all, it is plausible that spheroids with a merger 
origin would follow a j*—M+ relation that is parallel to 
that of spiral disks, but offset to lower by a factor of 
a few0 Second, the scatter in introduced by random 
merging may be relatively small. 

These two results in our toy model are both driven by 
the dominant cont r ibutio ns of or bital 7*. Sim i lar po ints 
were made by EH (|19790 and by IZurek et aLI (jl988D . in 
the latter case based on the prediction that A would be 
fairly constant with radius inside DM halos. The stars 
that condense at the center of a halo, and then partici- 
pate collisionlessly in its merger history, would naturally 
follow the same j—M scaling relations as the overall ha- 
los, modulo a smaller scale-length in converting from A 
to j (in Equation (|13|) . \E\ is inversely proportional to 
the radius). 

6.3.4. Evaluating the possibilities 

We now step back and consider how well the preceding 
evolutionary scenarios (outflows, stripping, collapse bias, 
and mergers) mesh with the observational constraints 
(Figures ITU and |2"0")) . The idea is to find a vector (or com- 
bination of vectors) that connects up the well-established 
endpoints in the j-M diagram: the ACDM halo initial 
conditions and the z — galaxy observations. It should 
however be remembered that the focus of this paper is 
not to solve long-standing questions about galaxy evolu- 
tion which may require a detailed understanding of the 
physics involved. Instead, our more modest goals are to 
illustrate how the j-M diagram can be used in practical 
terms as a constraint on theory, while looking for any 
hints as to the viability of various scenarios. 

Recent work in numerical simulations of disk galaxy 
formation has emphasized how outflows might remove 
low-jg material, which counteracts j-loss through tidal 
torques during g alaxy collapse, and maintains a high ne t 
level of fj (e.g., IBrook et alJl20ll IGuedes et alll20ll . 
We could then imagine that the differences between spiral 
and elliptical galaxies originate from the spirals having 
much stronger outflows at early times. 

This outflow scenario implies more mass loss in spirals 
and so would initially seem to work the wrong way in 
explaining the /* differences — but there could be other 
factors besides gas-depletion that affect /*. It is beyond 
the scope of this paper to explore this scenario in detail, 
but we emphasize that the focus on reproducing fj and /* 
for spirals needs to expand to include simultaneously the 
constraints from ellipticals, beyond these being nuisance 
factors that represent failed disks. 

We have already discussed how stripping before bary- 
onic collapse is not expected to produce large changes 
in the observable relations, which may indeed be 

part of the reason that there is not more scatter in these 
relations There is also a more obvious constraint that 

23 More generally, a similar slope would presumably be driven 
by any merger history that involves a scale-free mass spectrum 
of progenitors. This is a basic property of ACDM halos, but is 
incorrect at some level for stellar galaxies, owing to the strong 
break in their luminosity function. 

24 There is one case where severe stripping has apparently led 
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both spirals and ellipticals exist in the field and in clus- 
ters, so present-day environment cannot be the unique 
driver of morphology and j evolution. 

Collapse bias is an appealing possibility because it 
would provide a natural explanation for the positive cor- 
relation between /* and fj as in Figure [20] In this sce- 
nario, elliptical galaxies would cease to build up both M* 
and j* at relatively early times, with the remaining bary- 
onic M and j at late times either residing in a hot gas 
halo or having been blown out into intergalactic space. 
Spiral galaxies would have more protracted star forma- 
tion histories that increase M* and monotonically with 
time. 

Besides explaining the relative positions of ellipticals 
and spirals in the j'^-M* diagram, this scenario also fits 
in naturally with the observation that the stars in spirals 
are on average much younger than those in ellipticals. 
There may be additional implications if one connects the 
baryon collapse to the overall halo collapse, which has 
a well-understood theoretical underpinning. At a given 
z = mass, some halos should have collapsed earlier than 
others, leading to their DM distributions being more cen- 
trally concentrated. Given a fixed A, the central DM and 
associated stars would then have relatively low j values. 
Since halo collapse time is correlated strongly with en- 
vironmental density, one would then expect the low-j* 
galaxies to reside preferentially in high-density environ- 
ments - which is indeed what is found observationally 
(through the traditional morphology-density relation). 

A potential problem with this scenario is that it does 
not appear by itself to be capable of explaining the appar- 
ent deficit of in ellipticals, as discussed in Section lB".3.2l 
More detailed analysis would be needed to see if halo 
concentration makes a difference, and to understand the 
baryonic physics of why early-collapsing galaxies would 
also shut down their star formation more drastically than 
late collapsers. In addition to collapse bias, other effects 
may also need to be involved, such as a bias to low spin 
for their halos, or a component of real j loss. 

The merger scenario is a common explanation for el- 
lipticals, since it accounts for spheroidal morphologies 
through violent relaxation (|Toomrelll977t ). and because 
there is strong observational evidence for some ellip- 
tical galaxies actively fo rming through mergers (e.g., 
iRothberg fc Jo seph 2006j). Our toy model analysis sug- 
gests that the overall effect of mergers is to reduce the 
of the remnant relative to an initial j+-M* trend for 
disks, while the combination of multiple mergers may 
move the remnants parallel to that trend (Figure [Hfc)). 
This might provide a natural explanation for the ob- 
served j+-M+ trend for ellipticals: the slope, scatter, and 
offset relative to disks. Note that it is not entirely clear 
in this context why the spiral bulges and the ellipticals 
would follow the same j*-M* trends. 

A more quantitative comparison of our model to the 
observations allows us not only to constrain the typical 
mass ratios in mergers (as Figure [22]) . but also to infer 
the amount of mass growth in ellipticals since their as- 
sumed primordial disk phase. We do so by mapping our 
toy model vectors for mergers in the key fj-f+ diagram 

to a larg e red uction in j±: NGC 448 6B, which is a low- j* outlier 
in Figure ITU and is discussed in Romanowskv ct al. (2012). This 
"compact elliptical" is a fairly rare type of galaxy. 



(Figure I20[) . starting from initial conditions similar to 
present-day spirals (/* = 0.25, fj — 0.6), and requiring 
that they terminate at (/* — 0.1, fj = 0.1). 

Recalling that M v i r growth slightly outpaces 
growth we find that reducing by a factor of 2.5 re- 
quires a very long series of mergers, with a final growth 
factor of ~ 100 in M* and ~ 300 in M vu .. Consideration 
of the fj constraint then suggests a typical merger mass 
ratio of ~ 1:3. Such "major mergers" seem like a reason- 
able pathway to forming elliptical galaxies, although re- 
cent work sugges ts a more dominant role for minor merg- 
ers (e.g., ~ l:10:lNaab et al.ll2009[ iBezanson et all 1200a 
Khochfar et al.l 120111 : lOser et al.l 120121 : Uohansson et al.l 
20121) . which is motivated in part by explaining trends 
in size evolution, and is also supported by the observed 
shapes of rotation- velocity profiles (see Section 15.21 and 
I Arnold et al.l [20111 ) 1^1 This apparent tension is not of 
great concern since our current results involve significant 
observational uncertainties and a crude model for the 
merging vectors in Figure [2lT c). while not taking proper 
account of the redshift-dependence of virial quantities. 
However, they are intended to illustrate conceptually the 
kinds of constraints that are possible with more careful 
modeling. 

A merger scenario may successfully explain the j*- 
M+ properties of ellipticals, but it should be remem- 
bered that in a cosmological context, all galaxies in- 
cluding spirals should experience a continuous rain of 
accreting objects. Even if spiral galaxies have sys- 
tematically avoided the most extreme merger events, 
they will ha ve still experienced events i n the ~ 1:10 
range Ce.g.. iKauffmann fc White! 119931 : iStewart et al.1 
2008; Fa khouri et al.ll2010D . which as shown in our toy 
models could significantly reduce j*. A more detailed 
analysis of j+-M+ evolution within a cosmological frame- 
work is needed in order to investigate the quantitative 
differences that might arise between spirals and ellip- 
ticals owing to varying merger histories. In particu- 
lar, an explanation for the observed bulge-disk bi- 
modality is needed, since a spectrum of merger histo- 
ries is more suggestive of a smooth distribution of j+. 
It should also be kept in mind that (/*(M*)) is obser- 
vationally constrained not onl y for present-day galaxies , 
but also at earlier times (e.g.. iConrov fc Wcchslcrl [20091 : 
iMoster et al.ll2012j ). which introduces additional "bound- 
ary conditions" to j-M evolution. 

Synthesizing the scenarios above, it seems plausible 
that ellipticals might be explained through a combina- 
tion of collapse bias and multiple mergers-which bears 
a notable resembl ance to recent dis cussions of two-phase 
galaxy formation (|Oser et al.ll2010( ). In this context, an 
early burst of star formation would both imprint a rela- 
tively low initial j'*, and allow more opportunity for sub- 
sequent mergers to reduce j* further. Spirals would be 

25 In more deta il, th e fast- and slow-rotator subcategories of 
ellipticals (Section 14.21) are often thought to originate in diffcr- 
ent merger histories, such as binary versus multiple mergers (e.g., 
Burkert ct al. 2008; Bois ct al. 2011). Our discussion concerns pri- 
marily the fast-rotators, since these represent the vast majority of 
ellipticals, and in addition, our j+ constraints for the slow -rot ators 
are less certain. However, as discussed in Sections 15.21 and 15.31 
we detect no systematic difference in j^—M* space between the 
two galaxy types, suggesting that they may have relatively similar 
merger histories after all. 
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those systems where late gas infall both brings in higher 
j, and avoids the most active merging period. 

There are of course other considerations besides angu- 
lar momentum when constructing models of galaxy evo- 
lution, which are beyond the scope of this paper to eval- 
uate. We have also been able to cover only a subset of 
possible scenarios. 

One significant omissio n is the disk-instability pathway 
for bulge formation (e.g.. iToomre 1964; Dalca nton et al.l 
119971: Ivan den Boschl 119981 iParrv et al I I2009T ). which is 
an internal process where the bulge and disk either form 
from high- and low-j material, or else exchange j through 
gravitational torques. While this pathway is usually 
considered in connection with pseudo-bulges, there are 
recent proposals that the special conditions in high- 
rcdshift galaxy disks can lead to the massive, classical 
bulge s of present-day spirals, lenticulars, and ellipticals 
(e.g.,INoguchilll999tUmmeli et aLll2004t lElmegreen et al.l 
[2001 IDekel et al.l l2009aU bl: ICeverino et al.l 12010ft . The 
filamentary nature of mass and j inflows at high red- 
shift may also re quire significant revisions to standard 
spherical models (IDanovich et al . 2012; S ales et aLll2012t 
IDubois et al~ll2012t IKimm et al.ll201lD . 

Our overarching emphasis here is that whatever the 
mechanisms for galaxy formation, they must reproduce 
the basic j±-M+ scaling relations observed for both spi- 
ral and elliptical galaxies. A combination of all the pro- 
cesses mentioned above, and more, could be operational 
in real galaxies, where each process must be associated 
with a vector of j'^-M* evolution that is not arbitrary but 
physically-motivated, as we have sketched in Figures [50] 
and [T7] The sum of these vectors over the lifetime of 
the galaxy must preserve the halo-like scaling relations, 
along with a relatively small scatter. These may be very 
challenging constraints to match in practice, particularly 
if one includes boundary conditions on /^(M*) evolu- 
tion with redshift, and requires that the j+—M* relations 
hold for both bulge and disk components simultaneously 
within the same galaxies. 

Thus, a fresh approach to j—M analysis appears to 
hold promise for providing new, powerful constraints on 
galaxy evolution. We would encourage numerical simula- 
tors to keep this approach in mind as part of their toolkit, 
tracking the evolution of their simulated galaxies in the 
j—M diagram, while refining our schematic estimates of 
Aj-AM vectors, and thereby gaining more insights into 
the underlying physical processes in the simulations. 

7. SUMMARY AND CONCLUSIONS 

We have revisited the pioneering study of F83 which 
derived observational estimates for the fundamental 
quantities M* and j+ (stellar mass and specific angular 
momentum) of spiral and elliptical galaxies, and com- 
pared these to theoretical expectations based on hierar- 
chical assembly. Although the amount and distribution 
of j* in late-type galaxies has been an intensively-studied 
topic in the intervening years, even the most basic trends 
for early-types have not been satisfactorily established. 
We have capitalized on the advent of radially-extended 
kinematic data for a large sample of early-type galaxies, 
to update and extend the analyses of F83. 

We focus first on detailed analysis of a small sample 
of galaxies with data extending to typically five effective 
radii, which is the distance one must reach for a high de- 



gree of confidence in the estimates. We derive various 
formulae for use in quantifying for pressure supported 
systems, including deprojection effects. In order to esti- 
mate j+ for a larger sample of galaxies without requiring 
detailed modeling and data to very large radii, we test a 
simple, heuristic j*-estimator. 

Based on the shapes of observed rotation- velocity pro- 
files for the detailed sample of galaxies, we find that a 
convenient metric for the characteristic rotation velocity 
v s of a galaxy is provided by the observed rotation at a 
semi-major-axis distance of two effective radii. This ap- 
proximation is accurate at the level of ~ 0.1 dex, which 
is suitable for studying galaxy-to-galaxy variations in j* . 

We next assemble a large sample of galaxies in the 
nearby universe with adequate photometric and kine- 
matic data for estimating j* and M*. This sample covers 
the full spectrum of bright galaxy types from bulgeless- 
spiral to diskless elliptical, as well as a wide range in 
M+, centered approximately at the characteristic mass 
M*. We use our simple formula for estimating j*, while 
adopting simple bulge+disk models for the spiral galax- 
ies. 

Along the way, we also introduce an important new 
observational scaling relation for galaxies of all types: 
v s versus M*. This relation is analogous to the well- 
known Tully-Fisher relation for disk galaxies, but is more 
closely related to angular momentum than to dynamical 
mass. Unlike the generalized Tully-Fisher relation, the 
mass-rotation velocity relation shows near-perpendicular 
rather than parallel trends for spiral and elliptical galax- 
ies. These rotation- velocity trends combine with size- 
mass trends to trace the more fundamental j^-M* trends. 

Our combined j+-M+ estimates confirm the basic result 
of F83 that late-type spiral and elliptical galaxies follow 
parallel sequences of roughly a ~ 2/3 log-slope, but with 
a large zeropoint difference (in our analysis, the ellipticals 
have a factor of ~ 3-4 lower j* at a fixed M + ). Although 
this conclusion has already been used in some theoretical 
analyses, now it has a much firmer observational basis. In 
particular, the data do not support previous suggestions 
that major mergers have transported large amounts of 
angular momentum into the outer regions of ellipticals. 

We confirm for the first time that lenticular galaxies on 
average lie intermediate to ellipticals and late-type spi- 
rals in the j+-M+ plane, with tentative indications for two 
families of lenticulars characterized by low and high j*. 
We see no indication of systematic, overall differences be- 
tween centrally fast- and slow-rotator ellipticals. We also 
find that spiral bulges are consistent with following the 
ji,-M+ sequence for ellipticals, despite having very differ- 
ent relations between mass, size, and rotation. Thus, as 
far as the fundamental parameters j* and are con- 
cerned, spiral bulges are essentially like mini-ellipticals. 

We examine the residuals of the combined galaxy j'*— 
M+ data with respect to the disk-only trend, and find 
that these correlate better with disk-to-bulge ratio than 
with Hubble type. They also deviate from a lognormal 
distribution, possibly suggesting instead a bimodality in 
j+. Considering all of these results together, we pro- 
pose an alternative framework to the Hubble sequence, 
based on more physically motivated parameters. In this 
picture, all galaxies are a combination of a bulge and 
a disk, which are distinct subcomponents with different 
characteristic amounts of j*. Galaxy morphology may 
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then be seen as a secondary manifestation of the mix of 
high- and low-j material, or equivalently, the position of 
a galaxy in j±-M+ parameter space is a reflection of its 
bulge-to-disk ratio. 

We next connect our observational results to a theo- 
retical framework based on the hierarchical assembly of 
galaxy halos in a ACDM cosmology. We use numerically- 
informed analytic methods that are much simpler than 
hydrodynamical simulations, but less susceptible to the 
large, lingering uncertainties about baryonic recipes, res- 
olution effects, and other numerical issues. We find that 
the predictions for universal mean values of halo spin 
translate into j v i T -M v i r relations with an a = 2/3 log- 
slope, which is remarkably similar to the observed j*— 
M+ relations. The zeropoint differences among these re- 
lations provide valuable clues to the formation processes 
of different galaxy types. 

Mapping between halo and stellar quantities involves 
two basic parameters: the net fraction of baryons turned 
into stars, and the fraction of specific j retained, fj. 
We find that realistic variations of /* with mass produce 
surprisingly mild deviations of the j^-M* relation from a 
simple a = 2/3 power-law. The most noticeable correc- 
tion is a slightly shallower predicted slope for the spirals, 
which turns out to agree well with the observations. 

We explore two simplified alternative scenarios for 
explaining the spiral-elliptical dichotomy in the j+—M+ 
plane: the formation of spiral and elliptical galaxies 
in low- and high-spin halos, respectively (spin-bias sce- 
nario); and a difference in j retention (variable- /j sce- 
nario). We find that spin-bias does not explain the tails 
of the observed distribution, nor does it agree with 
the observed trend as a function of mass for the elliptical 
galaxies. The variable- fj scenario, on the other hand, 
matches the data well and suggests universal values of 
fj ~ 0.55 and fj ~ 0.1 for spirals and ellipticals, or for 
disks and bulges, respectively. The near-constancy of 
these values is intriguing, and means that all the com- 
plexities of galaxy evolution somehow effectively reduce 
to a simple model, where galactic stars have preserved 
the "initial" conditions of their host halos, including the 
j v ir~M v - lr slope and scatter. This interpretation may be 
useful for semi-analytically populati ng DM halos w ith 
both spiral and elliptical galaxies (cf. IMo et all lT998). 

Our fj result for spirals confirms similar conclusions 
going back decades, that these galaxies have retained 
most of their primordial specific angular momentum. 
This is an unavoidable conclusion from basic comparisons 
between observational constraints and theoretical expec- 
tations, which for many years presented a major chal- 
lenge to numerical simulations of galaxy formation, as 
these apparently predicted very low values for fj . It has 
gradually been realized that such simulati ons overpre- 
dicte d angular momentum transport (e.g., IKeres et"aT] 
2012), with major uncertainties still lingering in the bary- 
onic physics included in the simulations. 

Our consolidation of the fj result for elliptical galaxies 
provides a new benchmark for models of galaxy evolu- 
tion, which to be fully credible must reproduce the ob- 
served fj (and /*) trends for both spirals and ellipticals 
in the same simulation. For example, any feedback pro- 
cesses that are invoked to prevent overcooling and j loss 
in spirals should be much less effective for ellipticals. 

We have explored a few toy models for galaxy evolu- 



tion that exploit the basic constraints provided by the 
parameter space of and M*, or equivalently of fj and 
/*. Galaxies cannot evolve arbitrarily in this parameter 
space, requiring physically-motivated diagonal vectors of 
change (Aj*, AM*). Thus we suggest the j-M diagram 
as a key tool for assessing and interpreting any model of 
galaxy formation. 

Our simplified models suggest that a combination 
of early baryonic collapse and repeated galaxy merg- 
ing (major or minor) might account for the parallel- 
but-offset trend of ellipticals relative to spirals. We 
have provided illustrative constraints on the numbers 
and mass-ratios of the mergers, which after refinement 
with more detailed modeling could be compared with 
cosmologically-based predictions for mass growth and 
merging. 

In summary, we have established a new synthesis of 
j+-M* trends from observations, whose general resem- 
blance to halos in ACDM cosmology provides important 
support for that theory, and in turn furnishes a valu- 
able framework for constraining baryonic processes as 
discussed above. Our course, the observations presented 
here must be relevant to any model of galaxy formation, 
even if ACDM theory eventually needs revision. More 
generally, we propose that the morphologies of galaxies 
are closely tied to the evolution of angular momentum 
during their assembly, with late- types being very efficient 
at retaining j, and early- types proficient at shedding j. 

In this context, there are several areas ripe for ad- 
ditional progress. First, clear predictions from high- 
resolution cosmological simulations are needed for 
versus M* as a function of galaxy type, to explore 
whether the dichotomy between spirals and ellipticals, 
or disks and bulges, can be explained by differences 
in their assembly histories. Second, the observational 
work on nearby galaxies should be extended to the next 
level, via a volume-limited, homogeneous survey of all 
types of galaxies including full two-dimensional, wide- 
field photometric-kinematic bulge-disk decompositions, 
with attention to stellar mass-to-light ratio variations. 
This would allow for more robust analysis of the de- 
viations of the residuals from lognormality, and for 
more secure treatments of SO/Sa galaxies and of bulge 
and disk subcomponents. Third, the extensions of our 
study to lower-mass (e.g.. iKormendv fc B ender 20 121) as 
well as to higher- re dshift galaxies (cf. lPuech et al.ll2007t 
IBurkert et al.ll2010|). to freshly accreted material within 
galaxies (cf. IStewart et alj|2011l). and to the orientations 
of th e vectors (e.g.. iCodis et al.l 120121 : iTempel et al.l 
2012), would provide additional, valuable diagnostics. 
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APPENDIX 

APPENDIX A: ANGULAR MOMENTUM FORMULAE 
A.l. General formulae 

We begin with the description of a galaxy as the six-dimensional phase-space distribution function of its particles (gas 
and stars) /(r, v), where r and v are the (vector) position and velocity coordinates relative to the galactic center, and 
/ is normalized to unity when integrated over all positions and velocities. Given a total mass M, the three-dimensional 
spatial mass density at position r is 



p(v) = M I f(r,v)d 6 v, (Al) 

and the mean velocity at that position is 

M 

v(r) = — / v/(r,v)d a v. (A2) 



The true (i.e., not projected) specific angular momentum is then 

. _J t _ Jrx v(r)p(r)ri 3 r 

Jt " M ~ J P {t)(Pt ' (A6) 

Given the loss of information in observed galaxies (one positional dimension lost in projection, and two velocity 
dimensions usually unmeasurable as proper motions), one must adopt some simplifying assumptions in order to re- 
cover jt from observations. Our main assumptions here are that galaxies are transparent, have axisymmetric density 
distributions, and rotate on cylinders that are aligned with the symmetry axis of the density (with no other net velocity 
component such as expansion or contraction). 

Adopting cylindrical galactic coordinates (R,z,<f>), our modeling assumptions imply that v(r) becomes a simple 

rotation- velocity profile v TO t,t(R) <fii independent of z and (where <f> is the unit vector in the azimuthal direction), 
and that the density p(R, z) is independent of <j>. Equation (|A3[) then reduces to a one-dimensional integral: 

. _ Jv Iot , t (R)Z(R)R 2 dR 

Ji JZ(R)RdR ' [ V 



where 



E(R) = J p(R,z)dz (A5) 



is the surface mass density when the galaxy is viewed pole-on. With these assumptions, all galaxies with the same 
v T0 t,t(R) rotational profile and the same pole-on E(iJ) have the same jt, e.g., whether they are thin disks or round 
spheroids. 

Even under these fairly restrictive assumptions, recovering the true angular momentum of an observed galaxy is 
a difficult inverse problem. Fortunately, there is a way to structure the problem that makes it conceptually and 
computationally simpler. We separate the calculation for j t into two factors: 

jt = a JP - (A6) 



Here the second factor on the right is the analogue of Equation (|A4[) constructed purely from observations along the 
projected semi-major axis x: 

■ _ ^ Vrot, v {x)T l {x)x 2 dx ^ 
^ p J T,(x) x dx 

Note that this "projected specific angular momentum" is not literally the projection of jt on th e pl ane of the sky 
(which we will discuss at the end of this subsection). Also, the two denominators in Equations (|A4j) and (|A7j) are 
closely related. In the case of spherical symmetry for the density, they are identical, and for now we will adopt this 
assumption for simplicity, retur ning to the more general axisymmetric case later. 

The first factor in Equation (|A6[) . Cj, is a numerical coefficient incorporating all of the additional deprojection 
effects that depend on inclination and on the shapes (but not the amplitudes) of the surface density profile and of the 
rotation- velocity curve. Substituting Equations (|A4[) and (|A7j) into Equation (|A6|) . we have 

= /i; ro t,t(fl)S(fl) R 2 dR 
J Urot,p(aO x 2 dx 

Before manipulating Equation (IA8|) further, we point out that the advantage of this formulation of the angular 
momentum problem is that the complicated expression d need not be evaluated for every individual galaxy — provided 
that it is not sensitive to the details of the density and rotation-velocity profiles, and instead depends primarily on 
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Fig. 23. — Geometry and coordinate systems for the calculation of the deprojection factor Ci. Here the galaxy is spheroidal and 
axisymmetric, with an inclination of i = it/2 (edge-on). The left-hand diagram shows the observer reference frame (x' , y' , «'), while the 
right-hand diagram shows the galaxy reference frame (x,y, z). An example line-of-sight observation is illustrated by the X symbol on the 
left, and by the dashed line on the right leading to the O symbol. The arrows tangential to the circles show sample rotation- velocity vectors. 
More general inclinations are similar, but with the line-of-sight components of the rotation reduced by factors of sini. 



the inclination. Using simple models below, we will verify that this is the case, so tha t we can treat Ci as a numerical 
coefficient (calib rate d by models) that we combine with the observables in Equation (|A7|) . and thereby reconstruct j t 
using Equation (|A6[) . 

Next, to derive a more specific expression for Ci, we begin by considering a galaxy observed from an arbitrary 
direction, which we characterize by the usual inclination i, the angle between the rotation axis z and the line-of-sight. 
The observer's reference frame is (x' ,y' , z'), where z' is the line-of-sight (measured away from the observer), and x' 
is aligned with x for convenience (see Figure 123(1 . Since we are adopting an approach of modeling rotation- velocity 
profiles along the observed semi-major axis, we simplify the problem by doing calculations only for the case y' = 0. 
The coordinate transformations between galaxy and observer frames are then 

x = x' (A9) 

y = -z'ami (A10) 

z — z'cosi (All) 

(see lRomanowskv fc KochanekHl997l for more general expressions). The azimuthal unit vector in the galaxy frame can 
be expressed as 

<fi = sin <f>x + cos <j>y, (A12) 
which after using unit-vector relations analogous to Equations (|A9|) and (j AlOjl . becomes 

4> = sin <j> x' — cos 4> sin i z' . (A13) 

For a given parcel of material, the observer can measure only the projection of the mean velocity vector onto the 
line-of-sight, which we express as a dot product: 

v p = v(r) -z' = v Iot , t (R)j>-z'. (A14) 

From Equation (|A13[) . we then find 

x 

v p = - v I0 t,t(R) coscj) sin i = - v T otfi(R) sini, (A15) 

H 

/ 2 1.1/2 

where we use Equations (|A9I) and (|A10I) to make the substitution R = ( x 2 + z' sin 2 i ) . We thus define a projected 
profile of rotation velocity along the semi-major axis: 

v TO t, P (x) = ~ v p (x), (A16) 

which we calculate by integrating Equation (|A15[) along the line-of-sight z' while weighting by the density p(r). 
Recalling also that for now, we are assuming spherical symmetry for the density, we find the projected rotation- 
velocity profile 

x sin i J p{r) v mt ,t{R) Rrhte* 

v ^ ix) = jtm^ ' (A17) 
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Fig. 24. — Relations between true and observed properties of idealized spherical galaxies with flat intrinsic rotation-velocity profiles. Left: 
Ratio of observed to true rotation velocity, vs. galactocentric distance, in units of the effective radius. Three different density profiles are 
used, with different line-styles and their Sersic n indices labeled in the plot. Three inclinations are also used, with different colors, as labeled 
in the plot. The horizontal dot-dashed lines show the resu lts if the dilution effect is ignored. We have verified that we reproduce related 
results from the literature (Binncy f978; Young ct al. 1 978) . Right: Ratio of the true specific angular momentum to the value estimated by 
using the observed rotation- velocity profile [Equation 1A7H . as a function of the cosine of the inclination. The black solid curve includes 
the dilution effect, while the green dotted one does not. 

where we can also substitute r — ( x 2 + z' j . The denominator is the surface density [Equation (|A5I) ]. and we have 



v ro t,p(x) 



x sin i 

Y(xJ 



P{r) Urot,t(-R) 



dz' 



(A18) 



There are a couple of notable features about Equations (|A15|) and (|A18|) . One is that the difference between true 
and observed rotation velocity is more than a simple matter of the "sini" inclination effect for a galaxy tilted away 
fro m edge-on. For a system of finite thickness, there is an additional cos <p t erm that represents a "dilution" effect 
(cf. lBinnevl[T978t lYoung et af] 119781: iFillmore et "ail 119861: iNeistein et all 119991) . corresponding to the projection of a 
circular orbit seen at varying azimuth (f> (see Figure [23]) . Even for an edge-on case, the rotation vector is not in the 
line-of-sight except at the tangent point for a given semi-major axis distance x. This effect also implies that even if 
the true rotation velocity (v ro t,t) is constant with radius, the projected rotation velocity (i> ro t,p) is generally not, and 
goes to zero toward the cente r of the galaxy (see left-hand panel of Figure [24] for examples). 

Next, using Equation (IA18|i . we find our final value for the numerical coefficient, 



Ci 



J Urot , t (fl)S(fl) R 2 dR 
sini J x 3 dx J p{r) v TO t t t(R) R~ x dz' 



(A19) 



where p{r) can in principle also be written in terms of the directly observed profile E(a;) via an Abel integral. 

Having arrived at our general expressions for calculating j t , we note that there is another, related approach found in 
the literature. This is to construct an observatio nal analogue to .it i using the data directly to calculate a two-dimensional 
cross-product in an expression similar to (|A3I) (Franx 1988; Emscll ern et al.|[2007t Uesseit et al.ll2009D . This technique 
has the merit of potentially providing a general characterization of the data, including non-axisymmetric rotation 
fields. We have experimented with this approach and found that the fairly small improvement in accuracy (assuming 
the availability of full two-dimensional data) comes with the price of added complexity, while it also becomes unclear 
how to connect the results to the true three-dimensional angular momentum j t , which is our concern here. 

A. 2. Thin disks 

Following the general framework above, we continue with more specific galaxy models. For context, we begin with 
the simple case of an infinitely thin disk that has a constant intrinsic rotation- velocity profile, which is the same as 
the circular velocity, v ro t(R) = v c (R). It also has an exponential pole-on surface density profile: 



cx exp (-R/R d ) , 



(A20) 
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where Ra is the disk scale- length. Using Equation (| A4|l . we find 

it =2v c R d = 1.19 v c R c , (A21) 

where R c is the radius enclosing half of the light of a face-on disk, and we have used the relation R e = 1.68i?d for an 
exponential profile (note that for other disk inclinations, the observed half-light radius is smaller). Because the disk is 
infinitesimally thick, there is no dilution, except for a perfectly edge-on case, and we have simply Cj = (sini)" 1 . This 
means that the observed (constant) rotation velocity is v s — w c sini, and we can map between observations and true 
specific angular momentum by the relation 

2 Ve Rd , . , 

it - — A22 
smj 

Using the notation of Equation (| A6|) . we may say alternatively that j p = 2 v s i?d and = (sin The exc epti on is 
the p erfectly edge-on case (i — 7r/2), where we find by constructing the appropriate analogues to Equations (|A4j) and 
flSf that Ci = 2. 

Real disks do have finite thickness and consequently an appreciable amount of observed rotation-velocity dilution 
toward their centers, even for i < 90°. Still, detailed formulae like (|A6I) - (|A19[) are not really necessary, and as we verify 
in the main text, one can instead arrive at a very good approximation to it by taking the maximum or asymptotic 
rotation velocity observed at large radii to be v s , which is then used in Equation (|A22|I . 

For a disk galaxy it is fairly straightforward to estimate the inclination i based on the apparent axial ratio q, using 
the formula 

8 ini=(-— , (A23) 

where we adopt an intrinsic axis ratio of qi = 0.11 (e.g., lGuthrid lT992'). However, in order to tie in with spheroidal 
galaxies, we first treat i as an unknown for individual galaxies, so that a "typical" it is recovered by statistical 
deprojection of an ensemble of galaxies. One approach is to consider the median case C me d, where for randomly 
oriented galaxies, half of them have C, < C mec i; this is equivalent to identifying the median value for i. Since random 
inclinations in an axisymmetric geometry are distributed uniformly in cosi, the median is i — 7r/3, which corresponds 
to C mod =2/V3_~ 1.15. 

One may ask instead what is the average inclination correction, C avg , which if applied uniformly to all observed 
galaxies would provide the correct ensemble average for j t . Interestingly, the answer appears to depend on the exact 
question being p osed. One may start in the reference fram e of a galaxy (or planetary system) and ask what the average 
projection is (cf. iCumming et al.ll2008t lLaw et al.l 120091 ) . Or one may start with the observations and estimate the 
average devroiection fe.g.. lErb et al.1120061 ). In the spirit of mapping from projected to intrinsic quantities, we take the 
latter approach as follows. 

The average value of it is calculated as an expectation value of a probability distribution: 



Ut)= Jjt^dj t . (A24) 



We consider a single observation j p with infinitesimally small uncertainties, and since the probability P is uniform in 
cosi, we recast this equation (see Equation (|A6p ) as 

(Ci) = J C l sinidi. (A25) 

Given Ci — (sini) -1 for a thin disk, we then find C avg = (Ci) = tt/2 ~ 1.57. The equivalent calculation using the 
average projection yields C avg = (C^ 1 )^ 1 = 4/7r ~ 1.27. These two alternative correction factors differ by 23% or 
0.09 dex, which will be a relatively minor concern in our overall it analysis. In the main text, we analyze a sample 
of spiral galaxies where the deprojections are known; because this sample is chosen to be avoid galaxies that are near 
face-on, the different options for averaging or medianing imply very similar values for C, and we cannot distinguish a 
best choice. In any case, the important point is to carry out relative comparisons of disk and spheroidal galaxies using 
the same deprojection approach, which we choose to do using median inclinations. 

A. 3. Spheroids 

We next consider rotating spheroids, with th e go als of (1) calculating the coefficient Ci to be used in combination 
with observational estimates of j p in Equation (|A6I); (2) deriving an even simpler, more approximate expression for j p 
that avoids t he detailed calculation of Equation (|A7[) . We assume spherical symmetry for now, and adopt the general 
ISersid (|1968[) law for surface density profiles, 

E(R) oc exp \-b n (R/R c ) 1 / n ] , (A26) 

where the shape index n determines the steepness of the outer density profile (higher values are shallower: n = 1 is 
an exponential profile, n = 4 is a de Vaucouleurs law), and the numerical coefficient is 

b n ~ In - 1/3 + 0.009876/n (A27) 
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(|Marguez et al.ll2000D . Approximate analytical deprojection formulae for p(r) are also provided in the preceding 
reference. 

If we somehow knew the intrinsic rotation- velocity profile i> ro t,t(-R), it would be straightforward to evaluate Equa- 
tion (|A4j) and find jt- For example, a de Vaucouleurs (n = 4) model with constant intrinsic rotation velocity vt would 
yield: 

j t =2.29v t R e (A28) 

(cf. [Zasov 1985]). To deal wi th th e projection effects, we begin by adopting a constant intrinsic rotation velocity, and 
calculate Cj using Equation (|A19I) . We find that C, depends only on i, and not on n (Figure^U right). The inclination 
dependence is weaker than C\ oc (sini) -1 , which means it is partially counteracted by the dilution effect. An example 
is the edge-on case (i = tt/2), for which we find Cj = 1.18, as compared with Cj = 1 if dilution were neg lectedEl We 
performed the same tests for rotation- velocity profiles that vary smoothly with radius and found very si milar r esults. 
For example, Cj has only a ~ 10% dependence on the details of the rotation- velocity profile in Equation (|A19|) . 

The implication is that given observations of spherical galaxies, we can recover an accurate estimate of jt using 
Equation (|A6I) and (IA7|) . where Cj is a numerical factor that depends only on inclination. A simple expression for Cj 
that is accurate to better than 3% everywhere is: 

C * °-" + - 14 *, (A29) 
sin i 

where i in the numerator is in radians. Note that this expression for Cj differs from the undiluted value by no more 
than 20%, which might suggest that we ignore the dilution effect in our studies of j t , but we include it because we 
want to avoid collecting multiple systematic errors of this level. 

In practice, we do not usually know the individual inclinations of spheroidal galaxies, and instead need to do 
a statistical deprojection as introduced above for disk galaxies. Conside ring t he median case (i — 7r/3), we find 
C me d = 1.32. For the inclination averaged value, we again use Equation (IA25[) but with the spheroidal expression 
for Ci [Equation (|A19[) . where in practice only the denominator is affected by the inclination averaging]. We find 
C aV g = (Cr 1 ) -1 = 1-45; the alternative value is C avg = (Ci) — 1.73. 

When modeling real early-type galaxies as in the main part of this paper, we do not use any of these choices for C. 
This is because our galaxy sample is not randomly selected in sub- type (E or SO), and these sub- types are known to 
have a systematic connection with inclination. Lenticular galaxies are difficult to identify when close to face-on, so 
the samples of lenticular and elliptical galaxies will be biased to high and low inclinations, respectively. To correct for 
this bias, we adopt a simplified picture where there is only one species of early-g alaxy, which gets classified as E or SO 
depending on whether its inclination is below or above a boundary of i = 7r/3 (cf. IJorgcns en fc Franxl[l994[ ). Therefore 
we apply median deprojection factors of C me d = 1-21 (i — 1.32) for the lenticulars, and C mc d = 1-65 (i = 0.72) for the 
ellipticals. 

When applying this approach for j t estimation to real data, one has to make allowances for the limited radial extent 
of the available v IotyP (R) profiles. One can still use Equation (|A7|) and extrapolate beyond the data, while taking care 
to quantify the uncertainties that this entails. We will provide examples of this procedure in the main text. However, 
a much easier approach is also possible — sacrificing some accuracy for the sake of speed and simplicity — which we 
outline below and test with detailed calculations. 



A. 4- Simple angular momentum estimator 

The general idea for simple j p estimation is that every galaxy can be characterized by a single observed rotation 
velocity v s , which if substituted as a constant value in Equation (|A7|) would give the s ame answer as using the full 
v X ot,p(%) profile. The rotation- velocity and surface brightness components in Equation (IA7|) are then separable, and 
we can reduce the calculation to a product of a numerical coefficient, a velocity scale, and a scale-length 

j p = k n v s R c , (A30) 

which is analogous to Equation (|A22|) . He re j p is a general approximation of j p , and the coefficient k n is a spatial 
weighting factor calculated using Equation (|A7[) with v I0 t,p(x) set to unity (i.e., corresponding to a weighting derived 
from the stellar density profile combined with a radius x lever-arm). k n is a function of n, with the following handy 
approximation, accurate to better than 4%: 

k n ~ 1.15 + 0.029 n + 0.062 n 2 . (A31) 

Thus, fci = 1.19 and fc 4 = 2.29, as in Equations (|A21|) and (|A28j) . 

Th e cru x of this approach is estimating v s from observations without having to evaluate the full integral in Equa- 
tion (|A7[) . nor requiring that the rotation- velocity profile be known to the outermost radii. The trick comes from 
realizing that every galaxy has at least one radius x s where the local projected rotation velocity is equal to v s . This 



° That is, given an observed rotation velocity v s , we would 
have jt = 2.70v B R c , which may be contrasted with the expression 



jt = 1.03 v B R c from Scorza & Bender (1995), who took dilution 
and non-cylindrical rotation into account but neglected all of the 
angular momentum outside of 5 R c - 
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(x. R)/R. x s /R e 

Fig. 25. — Relations between true and observed properties of idealized spherical galaxies with de Vaucouleurs profiles (n = 4) and 
inclination i = it/3. Different colors indicate models with distinct rotation-velocity profiles (red: constant; blue: t> ro t,t(-R) °c R 1 ^ 2 , orange: 
Viot,t(R) oc i? -1 / 2 ). Left: True (dotted curves) and projected (solid curves) rotation velocities vs. galactocentric distance, in units of the 
effective radius. The horizontal dashed lines show the characteristic global rotation velocities v a , with the large dots marking the radii x s 
where v s = v TO t, p (x s ). Right: Ratio of estimated to true specific angular momentum, vs. the chosen measurement location x s (see text for 
details). 



radius varies from one galaxy to another, but we might expect the variation to be modest if real galaxies have rotation- 
velocity profiles that are not too dissimilar. The implication then is that we can simply meas ure the rotation velocity 
at the same radius x s /R e for all galaxies and adopt v s = v rot , p (x s ) for use in Equation (|A30j) . 

There are two potentially different goals here: one is to pick the radius x s that yields on average the correct answer 
for j p for a variety of galaxies; the other is to pick a radius where the scatter in the j p approximation is minimized. 
The anticipated origin of this scatter is the variety of radial behaviors of early-type galaxy rotation- velocity profiles, 
e.g., rising or falling from central to outer regions, such that the radius x s which applies to one galaxy might not work 
well for another. It is the issue of scatter that is most important, since a systematic offset in j p can be calibrated out 
and subsumed in the value of k n . 

We test this heuristic concept via some simple model galaxies below and via some real observations in the main text. 
We find that x s values anywhere in the range of ~ (2-5) R c appear to work reasonably well, and in order to maximize 
the number of suitable galaxies, we adopt x s ~ 2 R c . 

To arrive at this point, we first experiment with galaxy models having simple power-law rotation- velocity profiles 
v m t,t(R) that ran ge b etween cx R~ x l 2 and cx R +1 / 2 . For each Sersic index n, and rotation- velocity profile, we compute 
j p via Equation (|A7j) . next determine the equivalent v s value, and then find the radial location x s where the local 
projected rotation velocity is equal to v s . Finally, we examine the ratio jt/jt, [where jt is known exactly from 
Equation (|A4j) ] for a range of x s , which allows us to diagnose the best values of x s to use in general. We make jt and 
jt comparisons here, rather than j p and j p , in order to incorporate the impact of neglecting the mild dependence of 
Ci on rotation- velocity profile [Equation (IA19|) ]. 

This procedure is illustrated for one case (n = 4, i = 7r/3) in Figure [25] The left-hand panel shows the intrinsic 
and projected rotation-velocity profiles for three different models. The horizontal lines show the values of v s which 
would yield the correct j p = j p when substituting fc 4 = 2.29 in Equation (IA30[) . The intercept of the value with the 
corresponding v ro t tP (x) profile gives the appropriate radius x s . In the right-hand panel, the accuracy of jt is plotted 
versus the chosen x s for each model. For the constant- rotation- velocity model, jt — jt for a wide range of x s , i.e., the 
results are insensitive to x s . For the other rotation- velocity models, the choice of x s is more critical and ranges from 
~ 3R C to ~ 6i? c . A compromise radius that works reasonably well is x s ~ 4.5 R c , which provides accurate jt estimates 
at the 15% level or better for all three models. 

Exploring other values of n and i, we find that inclination is not important, and that some aspects of the n = 4 
example are generally true for other n. For each n, the radius x s that minimizes the errors in jt among the different 
rotation-velocity profiles is close to the exact value of x s for the constant rotation velocity case. This "sweet spot" 
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also depends systematically on the density profile: for n = 2, 4, and 6, it occurs at x s ~ 2 i? c , 4.5 i? G , and 10 -RP 7 !. 

This x s trend demonstrates that the extended envelopes of higher- n galaxies contribute relatively more to jt, and 
require correspondingly extended kinematics observations. Given that the R e values for high-n galaxies are already 
large, the observations to ~ 10 R c become extremely challenging, and we will in general omit such galaxies from our 
current analysis. Observations at x s ~ 4.5 R c as needed for n ~ 4 galaxies are also relatively rare, but fortunately the 
systematic bias in j t when using smaller x s is small, and the scatter is still only ~ ± 0.1 dex when x s ~ 2.5 R c . 

The foregoing tests involve model galaxies whose intrinsic properties we fully understand, but we also need to carry 
out checks with real gala xies . In the main text (Section [3]), we present more detailed models of a half-dozen real 
galaxies, using Equation (|A7j) along with the full observed rotation- velocity profiles to estimate j p . There we show 
that the constant v s approach works fairly well for a broad range of radii, and we adopt x s ~ 2 i? c ; at this radius, the 
systematic offset in j p appears to be small, and the galaxy-to-galaxy scatter is ~ ± 0.1 dexF^I 

We have thus arrived at our adopted procedure for simplified j p estimation using Equation (|A30[) . This j p is then 
multiplied by the deprojection factor Ci to yield an estimate of the intrinsic jt . Of special interest is the de Vaucouleurs 
profile (n = 4), with the approximate result for a median inclination 

jt=3.03v s R c , (A32) 

which is related to Equation (|A28j) by a factor of C mc d in order to account for the rotation- velocity dilution (v s < Vt ; 
note that as discussed in the previous subsection, in practice, C mo( j may need to be corrected for additional inclination- 
selection effects). This can be compared to the expression used by F83: 

j t =2.5v m R*, (A33) 

where v m was the maximum value of the projected rotation velocity within the (central) observed region. This formula 
is now superseded by the more accurate expressions above. 



A. 5. Flattening corrections 

Real galaxies of course are not spheres. We could in principle set up a standard axisymmetric density model by 
making an ellipsoidal substitution for r as needed with p(r): 

? /2 2 • 

m = R -\- — — = x + z' sir i H - , (A34) 

q z q A 

where q is the intrinsic axis ratio, and again we are measuring rotation velocity only along the projected semi-major 
axis (y' = 0) . Substituting this expression into Equation (|A17[) along with additional modifications of Equation (IA19I) , 
we find from some test calculations that even strong flattening makes only a very small difference to u ro t )P and to Ci, 
which typically increase and decrease (respectively) by ~ 5-10% relative to the spherical case for q = 0.3. Given the 
mildness of these effects, and the unknown systematics of the cylindrical-rotation assumption, we will not attempt to 
make any correction based on the flattening. The extrapolation of the kinematic data to larger radii is in any case 
probably the dominant uncertainty for our final results. 

The other potential concern here is the effective radius. We have already mentioned for disk galaxies that the 
value of i?d used in Equation (|A21|) is not the ob serve d but instead the deprojected value. Similarly, for flattened 
spheroidal galaxies, we should not use in Equation (| A30|) the circularized value R c normally tabulated in catalogs, but 
the equivalent distance along the semi-major axis, a c = Rc/y/q 1 , where q' is the observed axis ratio (in the rest of the 
paper we simplify this to q). This is because even though we may not know the inclination of a galaxy, we do know 
(for an axisymmetric case) that if it were face-on, the true R c needed to calculate jt would be roughly the same as a c . 
Since early-type galaxies can be as flattened as q' ~ 0.3, this is a significant correction. 

A caveat here is that the observed flattening is correlated with the inclination: e.g., a galaxy with q' — 0.3 is probably 
nearly edge-on, and adopting the median inclination for random orientations would then cause us to overestimate it s 
j t . Also, the intrinsic R e value for a spheroidal galaxy is slightly smaller than Oc/^/q 7 (see also iWilliams et aLl l2009'). 
However, these two effects are much weaker in general than neglecting the q' correction entirely. Therefore in the main 
text we carry out all of the calculations using densities and kinematics along the semi-major axis: e.g., and a e 

rather than £(i?) and R c . 



27 These locations are very similar to the corresponding radius ■ , , . , , , ., , . . . . ... 

r .. , , . , J ,,, r j. r comparison because we do not know the detailed intrinsic rotation- 

ot gyration (where a point mass would have the same moment ot , c , , , , , , , „ 

. ™. r ii n, 7 , i ii r.r>r,ni \ velocity profiles needed to calculate C; exactly. 

inertia as the mil mass prohle; sec Coccato ct al. 2009). 

28 For the real galaxies, we do not attempt a full jt and jt 
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Fig. 26. — Examples of kinematic model fits to two-dimensional velocity data in early-type galaxies. On the left are subset s of stellar 
kinematics measurements in the galaxy NGC 2768 from S AURON using integral field spectroscopy (Emscllem ct al. 2004) and from 
SMEAGOL using the SKiMS technique ( Proctor ct al. 2003). These data sets are marked with small and large points, respectively, and 
are taken from radii x ~ 35" (~ 0.3 a e ) and x ~ 125" (~ 1.3 a e ), respectively. The horizontal dashed line marks the systemic velocity 
(the SAU RON data have been offset to match up with SMEAGOL), and the curves show flattened sinusoidal model fits to the data [see 
Equati on l IBll ll. The in crease of rotation velocity with radius is apparent here by eye. The right panel shows PN data for NGC 4494 
UNapol itano ct al. 2009) at a radius of x ~ 60" (~ 1.3 a e ). The solid curve shows the mean velocity model fit, while dotted curves show 
the intrinsic dispersion a of the model, marking ±2<r boundaries. For clarity, the measurement uncertainties of the individual data points 
(typically 20 km s _1 ) are not shown. 



APPENDIX B: DETAILED OBSERVATIONAL RESULTS FOR INDIVIDUAL GALAXIES 

Here we present in some detail the data and methods used to derive angular momentum profiles for individual 
early-type galaxies. Our general approach is to map all of the data to an equivalent semi-major axis rotation-velocity 
profile, as though making a standard long-slit observation. In some cases, the data are already available in this format 
while in others that involve discrete and semi-discrete data at various position angles, we must fit a two-dimensional 
kinematic model before mapping onto the semi-major axis. 

Ou r modeling me thod was recently devel oped for use w ith s parsely sampled t wo-dimensional data (|Proctor et al.l 
120091 lArnold et all 1201 It iFoster et ail 120111: iStrader et ail 120111: iPota et alJl2012T ). It is based on flattened sinusoidal 
curves of mean line-of-sight velocity versus azimuth <j>: 



,x) = ±v rot ,p(x) 



Q 



-1/2 



(Bl) 



where v TOti p(x) is the rotation velocity versus the semi-major axis distance x, and v p is evaluated along an ellipse with 
axis ratio q. The on-sky coordinates are the same as introduced in Appendix^ but without primes (') for the sake of 
simplicity; the position angle of maximum receding rotation velocity is at <fi = <pQ. This model is exactly equivalent both 
to a cos0 rotation velocity multip lied by the equation for an ellipse in polar coordinates, and to the classic tilted-ring 
approach used for gas disks (e.g., iGuhattiak urta et al.l 1988; Begeman 1989). Therefore it is a natural observational 
model to use in conjunction with our j t modeling scheme, which assumes a cylindrical rotation field in the reference 
frame of the galaxy (Appendix [Aj) . There is also a connection here to the "kinemetric" methods of mo deling data from 
integral-field spectrographs , where the two-dimensional kinematics are expanded in a Fourier series (jKrainovic et al.l 
200(1 iKrainovic et al.ll2008l) . Our model is equivalent to the first-order term of this expansion. 

In principle, the parameters q and <po should be determined by fitting to the kinematic data, but in practice, the 
azimuthal sampling is often too sparse for such constraints. For the sake of a uniform treatment of the data, we adopt 
the photometric values for q and 4>o in every case; these are generally held constant at a global value, e.g., neglecting 
kinematic twists with radius. The errors introduced through these simplifications are generally much smaller than the 
uncertainties in extrapolating v ro t iP (x) outwards to radii beyond the data. 

Figure [55] illustrates the kinematic fitting procedure for several types of data. One type tha t may be unfamiliar 
is the "stellar kinematics with multiple slits' (SKiMS) t echnique from the SMEAGOL project (jProctor et all 120091: 
lArnold et al.l 120111 IFoster et al.l 1201 ll ; iForbes et al.l 120121 ) . which uses a wide-field multi-slit spectrograph to pro vide 
sparse sampling of a two-dimensional field as needed for reconstructing the basic kinematical structure. Equation (IB1|) 
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is fitted to the data (left-hand panel), while occasionally disregarding outlying data points so as not to skew the results. 
The uncertainties on the rotation velocities are estimated via Monte Carlo simulations of mock data sets. 

This procedure is modified slightly for discrete velocity data such as from PNe and GCs, which are sparsely sampled 
not only in position but also in velocity. Here, each data point is not a measure of the local mean velocity, but is a 
random discrete velocity drawn from an underlying distribution. The kinematic modeling of rotation therefore involves 
simultaneous fitting for the velocity dispersion using a maximum likelihood method (which represents a refinement , 
along with the rotation-field flattening, of the commonly-us ed least-squares method; e.g. , iCohen fc Rvzhovl 11997ft . 
Note also that when using data from the PN. Spectrograph (|Napolitano et al.l 120091 120111 : iCoccato et al.l 12009ft . we 
ignore the results inside ~ R c since there are known systematic biases that can occur in the central regions, where PN 
detection may be velocity-dependent. 

To combine different data sets on the same galaxy, wc could in principle fit them all simultaneously to a single 
kinematic model. However, for the sake of simplicity and transparency, we create separate one-dimensional rotation- 
velocity profiles from each data set separately, and then combine these into a single profile. For this final combination, 
we set up a new grid of data bins that are smoothly distributed in radius, and in each bin average the weighted 
contribution of each data point according to its uncertainties and its fractional overlap with that bin. The price of 
generating these well-behaved profiles is that the values in each bin are not fully independent of those in adjacent bins. 

A significant complication in our analysis comes in extrapolating the rotation-velocity profile beyond the outermost 
measured data points. In a spiral galaxy, this can be done assuming a profile that becomes asymptotically constant, 
but we cannot assume any such model for an elliptical or lenticular. Instead, we fit a power-law model of rotation 
velocity with radius to the final few data points: v rottP (x) cx x 1 . 

The optimum range of points to use in this fit is not straightforward to determine: use too few, and the extrapolation 
becomes merely an amplification of noise in the outermost measurements; use too many, and the inner measurements 
with their small uncertainties dominate the fit, which may not even agree with the outer measurements. Therefore 
some degree of subjectivity is involved, where the radial region used in the fit is determined by the widest range that 
still appears consistent with a smooth power-law trend, which typically occurs outside ~ 2.5 R c for our sample. The 
permitted range of power-law exponents 7 is then determined by a A\ 2 — 1 criterion. In some cases the formal range 
of 7 is large and we adopt a plausibility prior that the rotation velocity increases no more rapidly than linearly with 
radius (i.e. 7 < 1). Note that although there may be lingering concerns about the robustness of these extrapolation 
procedures, half of the galaxies in our detailed sample have data out to ~ 8-10 R e , in which case the extrapolation is 
relatively unimportant for the j* calculation. 

The kinematic data for the eight early-type galaxies that we model in detail are described in the remainder of this 
section, followed by plots of the kinematic profiles in the left panels of Figures [25] [3T] Here the velocity dispersion 
profiles are included for comparison with the rotation-velocity profiles, but are not otherwise used in this paper. For 
completeness, we also include the two spiral galaxies that we model in detail as described in Section [3TT1 (Figure [27 ]) . 

We convert the rotation velocity data into profiles of projected specific angular momentum j p using the methods of 
Section [2] and Appendix lAl These are also presented in Figures l27l I3T1 fright-hand panels). We summarize the relevant 
observational and modeling parameters for the se galaxies in Tabl e ECS Note that in many cases from this sample, new 
GC kinematics data were recently obtained bv lPota et al.l (|2012ft . but are not incorporated here. 

NGC 821 — Isolated, disky elliptical. Stellar kinematics data to ~ 4 a c are combined in a sin gle kinematic model , 
drawing on a number of sources including SAURON, SMEAGOL, HET, GMOS and M DM (jPinknev et ail [2051 
Proctor et al.l 120051: jWeiimans et all 120091: IProctor et ail 120091: iForestell fc Gebhardtl 12010ft . PN data are taken from 
Coccato et al.l ( 2009ft . who used non-parametric kinematic mapping techniques and found a strong kinematic misalign- 



ment relative to the stellar kinematics. We find the same result, such that the PNe between ~ 1 and 2 a e are strongly 
counter- rotating with respect to the stars. At larger rad ii, the PNe and stars agree, and we suspect some kind of 
contamination in the intermediate-radius PN sample (cf. Sambhu s et al.ll2006ft . However, the PN rotation velocity 
along the major axis is still marginally consistent (given the errors) with the stellar rotation velocity outside 1 a e , so 
we include these PN data anyway since the overall results are dominated by the st ellar data in these re gions. The final 
composite profile has a remarkably low outer rotation velocity, as emphasized by IProctor et al.l (j2009) . 

NGC 1400 — Apparently round early-type galaxy, probably a member of the NGC 1407 group. It may be a lenticular 
seen near fa ce-on. The data inc lude major-axis long-slit spectra out to ~ 1 a e , and SKiMS covering ~ 1-3 a e from 
SMEAGOL (jProctor et al.l l2009). As Proctor et al. discussed, there is an unexplained velocity dispersion discrepancy 
between these two stellar-light data sets, but the rotation- velocity profiles are consistent. 

NGC 1407 — Bright, round elliptical at the center of the Eri danus A group. Ste llar kinematics data out to ~ 3 a c 
come from major-axis lo ng-slit spectra and from S MEAGOL (jProctor et al.ll2009ft . Kinematics data for 55 metal-rich 
GCs extend to ~ 12 a c ()Romanowskv et ail 12009ft . where we have removed objects bri ghter than i = 21.3 owing to 
their peculiar kinematics, which may imply contamination by ultra compact dwarfs (cf. iBrodie et all 120 lift . There is 
also a mismatch between the stellar and GC velocity dispersions around ~ 1 a c that could be caused by additional 
contamination from substructure or from the metal-poor GC subpopulation, but the overall rotation- velocity profile 
is not affected. We adopt q = 1 and a fixed (f>o — 250° for all of the kinematic modeling. 

NGC 2768 — Flattened E/S0 in the field; we have overruled the RC3 classification of E6 in favor of other classifications 
from the literature, including the RSA, as SO. Stellar kinematics data out to ~ 2 a c are drawn from SAURON and 
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SMEAGOL (see left-hand panel of Figure |2"6")) . From our detailed early-type galaxy sample, this is the case with 
the most strongly increasing v TO t tP (x) profile. However, the kinematic data are still somewhat limited, and there are 
preliminary indications from more radially extended data for a decreasing rotation-velocity profile (A. Romanowsky 
et al.j in preparation). Note that the extrapolation of the data outwards can reach an unphysical rotation velocity, 
exceeding a plausible circular velocity of ~ 300 km s , but restricting the maximum rotation velocity does not 
substantially alter the permitted range of j p . 

NGC 3311 — Disky, low-luminosity elliptical in a loose group ( Leo I). The kinemat ic data are from SAURON inside 
0.7 a e , dee p long-slit VLT/FOR S2 major axis spectra to 3.4 a c (jCoccato et al . 1 120091 ) . and 112 PN velocities between 1 
and 12 a e (|Coccato et al.1120091 ). In comparing data-sets in their overlap regions (including other long-slit kinematics 
from the literature: Bender et al. 1994; Kormendy et al. 1998; Halliday et al. 2001; Prugniel & Simien 2002) we 
discovered some fairly large discrepancies, particularly in the velocity dispersion. Assuming that the higher spectral 
resolution of the FORS2 results makes them more reliable, we therefore subtract 59 km s _1 in quadrature from the 
SAURON dispersions, and multiply the SAURON rotation velocities by a factor of 1.13. The impact of this issue on 
rotation velocity and j p is not large, but we note that characterizations of low-mass galaxies in general using quantities 
such as Wrot/cp should be viewed with caution since they could be strongly affected by systematic problems, as in 
NGC 3377. There is some indication of a ~ 20° kinematic twist in the outer parts of this galaxy, where the true 
rotation velocity could be ~ 60 km s _1 rather than ~ 20 km s _1 . 

NGC 4314 (M84) — Bright, slow- rotator elliptical, apparently in a subgroup along with M86 falling in to the Virgo 
cluster. T he kinematic data inc lude SAURON inside 0.2 a e , deep long-slit VLT/F ORS2 major and minor a xis spectra 
to 0.8 a e (|Coccato et alJ 1200 9). and 450 PN velocities between 0.1 and 2.3 a e (jNapolitano et all 120111 ). As with 
NGC 3377 (above), we found discrepancies between the SAURON and FORS2 data, and subtracted 85 km s _1 in 
quadrature from the SAURON velocity dispersions. Both the FORS2 and PN data show a sudden kinematic twist from 
near major axis alignment inside 65" to near minor axis alignment at larger radii. (We recovered this information from 
the FORS2 data by fitting a two-dimensional ki nematic model t o the major and minor axis data.) This minor-axis 
rotation was also apparent in the PN analysis of iCoccato et all (|2009f ) . although the isophotes also twist so that by 
~ 200", the rotation is actually along the major axis. There may be further kinematic twisting outside ~ 250" such 
that the rotation is again along the minor axis, with a velocity of ~ 40 km s _1 . Because of the strong twisting in this 
triaxial galaxy, and the fact that q = 0.9-1.0 in the outer regions, we adopt the circularized R c = 176" value for a c . 
There is one other complication with a high-n galaxy like NGC 4374, as discussed in Section 13.21 the cumulative j p 
converges slowly with radius, and so we choose the virial radius r v ; r ~ 6095" of the galaxy as the boundary for defining 
"total" enclosed j p . 

NGC 4494 — Round ellip tical in the Coma I cloud. The kinematic data include stellar l ong-slit spectroscopy to ~ 2 a c 
(jNapolitano et all 120091) . stellar kine matic s with multipl e slits (SKiMS) to ~ 4 a c ((Proctor et all 120091) . a nd PN 
velocities to ~ 9 a„ (jNapolitano et alll200l ICoccato et alll2009f ). The newest SKiMS dataset from lFoster et all ((20111 ) 
was not used, but would yield essentially the same rotation- velocity profile; note that this paper provided an inclination 
estimate of i ~ 7r/4. An example of the PN data is shown in the right-hand panel of Figure |2"61 

NGC 5128 — Peculiar early- type gal axy at only 4 Mpc dista nce. The kinematic data include pioneeri ng two-dimensiona l 
spectroscopy out to ~ 0.5 a e from iWilkinson et al.l ()1986l ) and 780 PN velocities to ~ 15 a c from iPeng et al.l (|2004f l. 
Analysis of the PN data does not show any strong kinematic twists with radius, so for simplicity we use a fixed 
<fio = 259° for the entire galaxy. Also, because of the lack of detailed photometry for the galaxy, we set q = 1. 
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TABLE 3 

Angular momenta of galaxies modeled in detail. 



NGC 


Type 


D 


n 


a e 


Omax 




v 3 


Op 




Tracers 


Ref. 






(Mpc) 




(kpc) 


Oe) 




(km s _1 ) 


(km s — 1 kpc) 






3054 


Sb 


34.4 


1 


7.8 


2.3 


0.62 


177 


16701^ 


11.12 


Ha 


P+04a 


3200 


Sc 


52.2 


1 


15.3 


1.8 


0.28 


261 


4750+ 330 
£ " OU _280 


11.37 


Ha 


P+04a 


821 


E6 


23.4 


3.4 


4.5 


10.7 


0.60 


23 


210±igg 


11.02 


IFU, LS, MS, PN 


E+04, C+09, P+09, 
























W+09a, FG10 


1400 


SO 


25.7 


1.9 


2.7 


4.4 


0.89 


48 


11.05 


LS, MS 


P+09 


1407 


E0 


28.1 


4.3 


7.8 


11.7 


1.0 


39 




11.66 


LS, MS, GC 


P+09, R+09 


2768 


SO 


21.8 


2.6 


10.7 


2.3 


0.40 


169 


3060^° 


11.33 


IFU, MS 


E+04, P+09 


3377 


E5 


10.9 


2.0 


2.5 


10.5 


0.67 


56 


200 ± 30 


10.50 


IFU, LS, PN 


E+04, C+09 


4374 


El 


18.5 


8.3 


15.8 


2.3 


1.0 


32 


1610+ 3860 

1D1U -1000 


11.68 


IFU, LS, PN 


E+04, C+09 


4494 


El 


16.6 


3.2 


4.5 


8.5 


0.84 


58 


4301^ 


11.08 


LS, MS, PN 


N+09, C+09, P+09 


5128 


SO 


4.1 


4 


6.0 


15.3 


1.0 


61 




11.28 


MS, PN 


W+86, P+04b 



Note. — Morphological Types are generally taken from the RC3 catalog ( dc Vaucoulcurs ct al. 1991). v s and j p are calculated 
from the full, detailed treatment of E quations f3i and (2J. The distances D are taken from surface brightness fluctuation analyses 
(Tonry ct al. 2001; Blakcslcc ct al. 2009) where available, and otherwise from redshifts. The Sersic parameters n and a e are taken from 
various literature sources. The tracers include Ha gas-emission rotation- velocity curves; long-slit stellar kinematics (LS); multi-slit stellar 
kinematics (MS); inte gral-field stellar kinematics (IFU ); plane tary nebulae (PN); and metal -rich globular clu sters (GC) . The references are 
IWilkinson et alj|1986l (W+86); lEmsellem et~aT][2004l (E+04); IPizzella et al.ll200l (P+04a); IPeng eF al. 2004 ( P+04b); ICoccato et alll2009l 
(C+09); INapolitano et~aLll2009l fN+09); IProctor et a l. 2009 (P+09); IRomanowskv ct al. 2001 (R+09); Wciimans ct al. 2009 (W+09a); 
ll'orestell fcTIebl iardt 2010 (FG10). 
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Fig. 27. — Kinematic results for the spiral galaxies NGC 3054 (top panels) and NGC 3200 (bottom panels). The left-hand panels show 
data for the profiles of projected rotation velocity (black filled circles) and velocity dispersion (purple open circles) vs. semi-major axis 
radius. The black curve with shaded region shows the outer power-law extrapolation of the rotation velocity and its estimated uncertainty. 
The right-hand panels show the data converted into a profile of cumulative projected specific angular momentum vs. radius (solid curve), 
with the shaded regions showing the uncertainties (barely visible for this galaxy). The dotted curve shows the expected distribution of 
j P (< x) for a constant rotation- velocity profile. Vertical dashed lines show the effective radius of the galaxy and the maximum extent 
probed by the kinematic data. 




Fig. 28. — Kinematic results for the early-type galaxies NGC 821 (top panels) and NGC 1400 (bottom panels). See Figure [27] for 
explanations. In addition, horizontal lines in the left-hand panels indicate the radial ranges spanned by the different data sets. Note also 
the dashed line at zero rotation velocity in the case of NGC 821. 
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Fig. 32. — Relation between Sersic index n and stellar mass M* for early-type galaxies. Points show results from the literature for 
individual galaxies in our sample. Symbol and colors indicate different galaxy types as in the legend; fE and sE mean 'fast rotator' 
and 'slow rota tor' ellipticals based on their central properties (see text for details). The black curve shows our adopted mean relation, 
Equation {CI) . 

APPENDIX C: SUPPLEMENTAL OBSERVATIONAL INFORMATION FOR FULL SAMPLE 

Here we describe some details of the full sample of observational data, with results reported in Tables CH] and CJSJ 
The methods for the spiral galaxies are described in Section I4.ll and do not require further elaboration. For the early- 
types, in addition to the generalities in Section l4~2l here we discuss details including estimates of the Sersic index n, 
an evaluation of bias in our galaxy sample selection, and some nuances of the size-mass relation. 

The issue with n is that, as discussed in the main text, its value is a significant factor in accurate j calculation, but 
is not available for some of the galaxies in our sample. Instead, we use a statistical estimate based on trends among 
other galaxies. Figure [32] shows the n-M* correlation for the early- type galaxies in our sample that do have Sersic fits 
available. There is a correlation betwe en M+ and n that echos oth er trends in the literature, and in particular matches 
the power-law slope of ~ 0.3 found by IGraham fc Guzman (2003]). The mean relation that we adopt is 

n = 5.4 x 10~ 3 x (M*/M Q ) - 25 . (CI) 

There is a scatter of a n ~ 1 about this relation, which is also similar to variations found in the literature for d ifferen t 
Sersic fits of the same galaxies. This uncertainty in n translates to an uncertainty of ~ 25% in k n [see Equation (|A31j) ]. 
and thus ~ 0.1 dex uncertainty in j [Equation (|6j)]. 

Next, to investigate whether or not our galaxy sample provides an unbiased, fair representation of the nearby 
universe, we could check a variety of different galaxy parameters, but given our focus on angular momentum, the most 
relevant parameter is the readily-measured central stellar rotation velocity. 

Our point of comparison is the ATLAS 3D survey, which presented photometric and kinematic pro perties for a com- 
plete, volume-limited sample of 260 nearby early-type galaxies, selected only by K-band luminosity (jCappellari et all 
l2011a|) . We use their tabulated central rotation metric, (v/cr) c /2, which is flux- weighted in two dimensions within 
an aperture of R e /2 (jEmsellem et al.l 1201 lh . From the galaxies in common between our samples, we calibrate this 
metric to our own major-axis kinematic data measured locally at a c /2. We find that the ATLAS 3D values of (u/cr) e / 2 
should be multiplied by a factor of ~ 2.2 to match our data (which, as expected, is a large difference because of the 
two-dimensional versus major-axis measurement locations). 

We next divide (v/a) from both data sets by a factor of [e/(l — e)] 1 ^ 2 to derive a rotation-dominance parameter 
(v/a)*. This is analogous to the standard parameter in Equation ©, but is constructed differently for the sake of 
convenient comparisons between the two data sets, and should not be used to make detailed inferences about the 
properties of these galaxies. 

We plot (v/a)* versus stellar mass in Figure [33] Our sample follows a similar trend to the unbiased sample, with 
central rotation becoming systematically less important for the more massive galaxies. Both samples span a similar 
range in mass. 

Another important issue here is the inclusion of different sub-classes of early-type galaxies. Our sample consists of 
15 lenticulars and 25 ellipticals, which implies a strong selection bias for ellipticals, since ~ 75% of nearby early-types 
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Fig. 33. — Rotation dominance parameter vs. stellar mass for early-type galaxies (see text for details). Our galaxy sample is shown in 
the top panel, and the ATLAS 3D sample in the bottom panel. Symbols show different galaxy types as in the legend. Our sample is biased 
toward ellipticals and more massive galaxies but otherwise appears to provide a fair representation of the central kinematics of nearby 
early-types. Note that the fraction of "slow rotators" is highest for the most massive systems. 



are SOs (jCappellari et al.ll2011al) . 

In Section B~B1 we compare the size-mass relation of our sample (both the early- and late- types) with the ATLAS 3D 
results. For ATLAS 3D , we adopt the power-law relations provided by iCappellari et all ()2011aD . after adapting them 
to the parameters that we plot for our sample. For the spirals, we convert from circularized R c to a c by assuming a 
median inclination of 60°. For the fast-rotator E/SOs, we use a median axis ratio of q = 0.63. 

Comparison with our sample (in Figure fTT)) is not straightforward since we have used disk sizes rather than total 
galaxy sizes (bulge+disk) for the spirals. In any case, we find that our results are consistent overall with ATLAS 3D 
but there appears to be a difference for the lower-mass ellipticals; at log (M±/Mq) ~ 10-10.5, these galaxies are ~ 30% 
more compact in our sample than the trend found by ATLAS 3D . Examining several galaxies in common between the 
samples, we find that t he discrepan cy persists for ca s es whe re we have made use of high-quality photometric analyses 
(|Kormendv et al. 2009). Figure 8 of lCappellari et all (|2011aD suggests that for the smaller galaxies, the ATLAS 3D sizes, 
which are partially based on NIR measurements, are overestimated relative to optically-based sizes. This illustrates 
the type of systematic uncertainties remaining in our analysis, and also of the comfortably low impact that they have 
on our analyses (only ~ ± 0.05 dex or ~ ± 15% in j'*). 

One additional issue is that our sample of spirals was selected to avoid strong bars which might in principle introduce 
a bias in j*. Inclusion of barred spirals could be an improvement made in a next generation, volume- limited j^—M* 
survey. 



58 



Romanowsky & Fall 



TABLE 4 

Angular momenta of spiral galaxies 



Name 


Type 
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(Mpc) 
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i 

(deg 
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JP 

(km s~ 1 
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it lo£ ( 
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NGC 224 


Sb 


3.0 


0.70 


0.19 


5.9 


0.17 


0.70 


77 


234 


73 


2230 


2290 


11.06 


NGC 247 


Sd 


6.9 


3.8 


0.00 


4.1 


— 


— 


76 


92 


— 


770 


790 


9.84 


NGC 300 


Sd 


6.9 


1.8 


0.00 


1.6 


— 


— 


47 


60 


— 


190 


260 


9.23 


NGC 701 


Sc 


5.0 


24.3 


0.00 


1.9 


— 


— 


61 


127 




480 


550 


10.46 


NGC 753 


Sbc 


4.9 


65.8 


0.04 


4.8 


0.02 


0.86 


50 


161 


31 


1480 


1930 


11.31 


NGC 801 


Sc 


5.3 


77.3 


0.21 


13.1 


1.1 


0.75 


86 


205 


56 


4270 


4280 


11.50 


NGC 1024 


Sab 


2.4 


46.8 


0.34 


6.2 


0.09 


0.56 


08 


252 


101 


2070 


2230 


11.34 


NGC 1035 


Sc 


5.0 


16.4 


0.00 


1.3 


— 


— 


71 


126 




340 


360 


10.13 


NGC 1085 


Sbc 


3.6 


89.9 


0.22 


6.0 


0.2 


0.90 


39 


194 


68 


1810 


2900 


11.42 


NGC 1087 


Sc 


5.2 


20.1 


0.00 


3.0 




— 


50 


103 




630 


820 


10.62 


NGC 1325 


Sbc 


4.0 


20.9 


0.03 


1.2 


0.3 


0.82 


70 


169 


56 


1360 


1450 


10.61 


NGC 1353 


Sb 


3.3 


20.0 


0.15 


2.8 


1.0 


0.74 


65 


205 


31 


990 


1090 


10.74 


NGC 1357 


Sab 


1.9 


26.4 


0.34 


4.7 


2.3 


0.85 


•15 


184 


34 


1210 


1720 


11.03 


NGC 1417 


Sb 


3.3 


54.7 


0.07 


6.4 


0.2 


1.0 


50 


218 





2580 


3370 


11.31 


NGC 1421 


Sbc 


4.1 


27.7 


0.05 


6.0 


1.0 


0.73 


78 


192 


79 


2230 


2280 


10.93 


NGC 1620 


Sbc 


4.5 


47.2 


0.07 


7.6 


0.7 


0.85 


69 


224 


34 


3150 


3380 


11.32 


NGC 2403 


Scd 


6.0 


3.2 


0.00 


2.0 


— 


— 


69 


116 


— 


470 


500 


10.01 


NGC 2590 


Sbc 


3.6 


71.8 


0.25 


5.9 


0.5 


0.54 


75 


246 


110 


2220 


2300 


11.34 


NGC 2608 


Sb 


3.3 


33.1 


0.07 


2.9 


0.4 


0.86 




107 


37 


580 


760 


10.67 


NGC 2639 


Sa 


0.9 


49.7 


0.61 


2.8 


1.8 


0.79 


65 
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69 


1110 


1300 


11.33 


NGC 2708 


Sb 


3.1 


30.2 


0.08 


3. 1 


0.3 


1.0 


61 
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1430 


10.81 


NGC 2715 


Sc 


5.2 


23.2 


0.01 


5.6 


0.03 


0.87 


70 


136 


37 


1520 


1620 


10.83 


NGC 2742 


Sc 


5.3 


22.6 


0.01 


3.5 


0.14 


0.89 


58 


148 


15 


1030 


1210 


10.65 


NGC 2775 


Sab 


1.7 


21.5 


0.20 


1.5 


0.4 


0.89 


44 
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11 


1400 


2020 


11.31 
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NGC 3898 


Sab 


1.7 


22.7 


0.66 


3.7 


2.1 


0.62 


67 


230 


107 


930 


1070 


11.07 


NGC 4062 


Sc 


5.3 


10.8 


0.03 


2.0 


0.50 


0.96 


61 


139 


13 


550 


610 


10.22 


NGC 4236 


Sdm 


8.0 


4.4 


0.00 


3.7 


— 


— 


72 


79 


— 


580 


610 


9.43 


NGC 4258 


Sbc 


4.0 


7.3 


0.00 


6.0 




— 


69 


186 


— 


2240 


2410 


11.08 


NGC 4378 


Sa 


1.0 


41.1 


0.52 


6.0 


2.8 


1.0 


35 


185 





1060 


1850 


11.42 


NGC 4419 


Sa 


1.1 


12.5 


0.10 


1.2 


0.07 


0.92 


71 


180 


19 


370 


400 


10.43 


NGC 4448 


Sab 


1.8 


7.0 


0.22 


1.0 


0.03 


0.90 


69 


182 


38 


300 


320 


10.03 


NGC 4594 


Sa 


1.1 


12.7 


0.86 


3.1 


8.4 


0.70 


81 


348 


101 


2030 


2380 


11.56 


NGC 4605 


Sc 


5.1 


4.5 


0.01 


0.8 


— 


— 


68 


79 




130 


140 


9.58 


NGC 4682 


Scd 


5.9 


34.7 


0.02 


4.4 


0.04 


0.77 


57 


143 


10 


1220 


1460 


10.73 


NGC 4698 


Sab 


1. 1 


13.7 


0.55 


3.6 


1.6 


0.96 


70 


233 


18 


790 


850 


10.96 


NGC 4736 


Sab 


2.1 


4.8 


0.35 


1.5 


0.14 


0.90 


30 


89 


23 


180 


360 


10.76 


NGC 4800 


Sb 


3.0 


18.4 


0.09 


1.1 


0.01 


1.0 


53 


136 





260 


330 


10.54 


NGC 4845 


Sab 


2.3 


24.9 


0.07 


5.2 


0.3 


0.80 


72 


171 


25 


1650 


1740 


11.10 


NGC 5033 


Sc 


5.1 


17.2 


0.10 


6.9 


0.10 


0.70 


69 


195 


57 


2440 


2620 


11.31 


NGC 5055 


Sbc 


4.0 


8.7 


0.05 


4.0 


0.2 


0.86 


51 


169 


27 


1280 


1590 


11.11 


NGC 6314 


Sa 


1.1 


98.2 


0.56 


6.2 


2.4 


0.70 


70 


215 


83 


1430 


1560 


11.50 


NGC 7171 


Sb 


3.1 


39.3 


0.04 


4.0 


0.04 


0.72 


52 


169 


35 


1300 


1650 


10.84 


INCiC 7217 


Sab 


2.5 


16.4 


0.25 


2.5 


2.9 


0.93 


35 


119 


23 


590 


1030 


11.07 


NGC 7331 


Sb 


3.9 


14.4 


0.26 


4.8 


0.2 


0.40 


67 


202 


111 


1450 


1580 


11.34 


7537 


Sbc 


3.6 


37.3 


0.18 


2.5 


0.2 


0.70 


79 


137 


34 


570 


580 


10.47 


NGC 7541 


Sbc 


4.7 


37.5 


0.01 


4.7 


0.3 


0.62 


72 


219 


34 


2050 


2150 


11.23 


NGC 7606 


Sb 


3.0 


31.3 


0.08 


5.8 


1.22 


0.60 


66 


238 


79 


2540 


2780 


11.34 


NGC 7664 


Sc 


5.0 


47.9 


0.03 


2.3 


0.2 


0.80 


58 


157 


36 


710 


830 


10.86 


IC 467 


Sc 


5.2 


32.1 


0.00 


3.7 






67 


129 




960 


1040 


10.44 


IC 724 


Sa 


1.0 


89.5 


0.41 


8.2 


1.0 


0.78 


55 


233 


86 


2360 


2890 


11.51 


UGC 2259 


Sdm 


7.8 


9.8 


0.00 


1.1 






22 


31 




70 


190 


9.14 


UGC 2885 


Sc 


5.2 


78.0 


0.14 


15.1 


10.4 


0.88 


65 


281 


76 


7560 


8380 


11.74 


UGC 11810 


Sbc 


3.7 


66.4 


0.05 


5.8 


0.4 


0.76 


74 


183 


22 


2010 


2090 


10.82 


UGC 12810 


Sbc 


3.8 


109.6 


0.07 


8.0 


0.7 


0.83 


70 


211 


67 


3140 


3350 


11.26 


Note. - 


Morphological Types are 


taken 


from the RC3 


catalog jde 


Vaucoulcurs ct al. I199UV 


The Hubble stage T[j u bble 


is from 



HyperLeda (Paturcl ct al. 2003). Distances are from NED, using the redshifts with a Hubble constant of h = 0.73 and correct ing the 
Hubble flow for Virgo, Shapley, and the Great Attractor. B/T is the r-band bulge-to-total luminosity ratio from Kent (1986, 1987, 1988); 



the B/T value for NGC 4736 was not tabulated, and we have estimated it from the rotation curve decomposition. The semi-major axis 
disk and bulge scale- lengths (R^ and a e b) are from parametric fits to the Kent surface brightness profiles. The bulge axis ratio % is also 
from K ent, except for NGC 4736 IIMendez-Abreu et al.ll2008l) . The inclination i is from lRubin et alj <[T980| . [1983. IT985D or lWhitmore et al.l 
(1984), except in some cases where it is derived directly from the observed disk flattening. The characteristic disk rotation velocity v s d 
is based on the (2—3) i?d regions of the Kent deprojected rotation curves, reprojected to the observed values. The bulge rotation velocity 
v s t, is estimated using a simplified model based on observed ellipticity (see main text). The projected and intrinsic specific angular 
momenta j p and jt are derived using Equations (|6jl and J2J, respectively. The stellar mass M* is based on 2MASS S'-band photometry 
(see main text). 
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Approximate angular momenta of early- type galaxies 
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NGC 821 


E6 


-4.8 


23.4 


3.4 


4.5 





60 


20 


180 


11.02 


LS, MS, PN 


C+09, P+09, 
W+09a, FG10 


NGC 1023 


SO 


-2.7 


11.1 


3.9 


3.2 





34 


175 


1270 


10.97 


LS, PN 


C+09 


NGC 1316 


SO 


-1.8 


21.0 


6.5 


13.5 





68 


70 


3640 


11.94 


PN 


M+12 


NGC 1339 


E3 


-4.3 


19.7 


4.3 


1.5 





71 


105 


400 


10.52 


LS 


G+98 


NGC 1344 


El 


-4.0 


20.9 


4.1 


5.1 





58 


20 


240 


11.13 


PN 


C+09 


NGC 1373 


E2 


-3.9 


19.3 


4.3 


1.0 





77 


10 


25 


9.79 


LS 


G+98 


NGC 1379 


E0 


-4.8 


19.5 


2.3 


2.3 





97 


20 


70 


10.68 


LS 


G+98 


NGC 1380 


SO 


-2.3 


21.2 


3.3 


5.4 





51 


220 


2330 


11.34 


LS 


B+06 


NGC 1381 


SO 


-2.2 


21.9 


3.1 


3.6 





26 


155 


1040 


10.75 


LS 


W+09b 


NGC 1400 


SO 


-3.7 


25.7 


1.9 


2.7 





89 


50 


3100 


11.05 


LS, MS 


P+09 


NGC 1404 


El 


-4.8 


20.2 


3.7 


2.6 





88 


115 


630 


11.26 


LS 


G+98 


NGC 1407 


E0 


-4.5 


28.1 


4.3 


7.8 


1 





20 


680 


11.66 


MS 


P+09 


NGC 1419 


E0 


-4.8 


22.9 


3.6 


1.1 





99 


3 


7 


10.18 


LS 


G+98 


NGC 1428 


SO 


-3.0 


20.7 


1.5 


1.4 





63 


65 


130 


9.97 


LS 


G+98 


NGC 2310 


SO 


-2.0 


14.9 


2.1* 


3.4 





19 


135 


680 


10.35 


LS 


W+09b 


NGC 2577 


SO 


-3.0 


31.7 


2.7* 


2.8 





62 


210 


1000 


10.79 


LS 


R+99 


NGC 2592 


E2 


-4.8 


31.7 


2.6* 


2.2 





83 


160 


610 


10.76 


LS 


R+99 


NGC 2699 


El 


-5.0 


27.9 


2.4* 


1.7 





93 


90 


240 


10.55 


LS 


R+99 


NGC 2768 


SO 


-4.4 


21.8 


2.6 


10.7 





10 


185 


3360 


11.33 


MS 


P+09 


NGC 2778 
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Note. — Mor phological Types ar e generally taken from the RC3 catalog (dc Vaucoulcurs ct al. 1991). The Hubble stage Th^mc is 
from HyperLeda llPat urcl ct al. 2003). The distances D are taken from the surface brightness fluctuation analyses {Tonry ct al. 2001; 
Blakcslcc ct al. 2009) where available, and other wise f rom rcdshifts. The parameters of the stellar mass profile (n,a c ,q) are taken from 
various literature sources as described in Section 14.21 Galaxies with a '*' by their Sersic index n have this value predicted from their 
stellar masses, rather than being measured directly (see text for details). v s is estimated from the observations at x s = 2 a e , and j p is 
calculated using Equation l[6j . For the subset of 8 galaxies with detailed j p calculations listed in Table [3] the j p values are listed here 
for reference, but are not used in our main analyses. The observational tracers are: long-slit stellar kinematics (LS); multi-slit stellar 
kine matics (MS); integ ral-field stella r kinematics (IFU); planetary ne bulae (PN); and metal-rich gl obular c lusters (GC). The r eferences 
are: IDavies et al.|[l983T(D+83'); ISimien fc Prugnicl 1997a. b (SP 97a,b) ; ISanchez-Blazquez et al.| | 2007| (S+97); IGraham et a.1.1119981 (G+98); 
Rix et al.1119 99 (R+99); Halliday ct al. 20 0j (H+01) ; IPeng et jO . 2004 fP+04b ^pedregaI et aUl2006l (B+06):IDouglas et al.lboO'l (D+07); 
Napolitano ct al. 2009 (N+09); Coccato ct al. 2009 (C+09); Proctor ct al. 2009 (P+09); Wciimans ct al. 2009 (W+09a); Williams ct al. 
2009 (W+09b); Forestcll & Gcbhardt 2010 (FG10); Spolaor et al. 2010 (S+10); Arnold ct al. 2011 (A+ll); McNcil-Movlan ct al. 2012 
(M+12). Other photometric references are: IMahdavi et al.H2005T 7M+05). 
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Fig. 34. — Bulge-fraction vs. morphological type for early-type galaxies (left panel), and specific angular momentum vs. stellar mass for 
the disk and bulge components of disky ellipticals (right panel). On the left, dotted curves show two alternative observational trends from 
Simicn & dc Vaucoulcurs (1986). Th e points show data from ATLAS 3D galaxies, using Equation IIDll l to estimate B/T with C = 1.35. 
The symbol typ es a re as in Figure [33] The dotted lines at right show the best-fit relations for spiral disks and ellipticals from the right-hand 
panel of Figure 1141 

APPENDIX D: DECOMPOSING EARLY- TYPE GALAXIES 

In Section 15.21 we sought to understand the j*-M+ trends for spiral galaxies by considering their bulge and disk 
components separately. Many early-type galaxies (both lenticulars and fast-rotating ellipticals) are also thought to 
consist of such subcomponents, but the decompositions tend to be more difficult. Here we develop a novel method to 
estimate the bulge fraction, /b = B/T, in order to study correlations involving this parameter, while we also analyze 
a small set of ellipticals with decompositions from the literature. 

First, we suppose for simplicity that the relative dominance of rotation observed in an early- type galaxy indicates the 
disk fraction, i.e., the bulge is assumed to be non- rotating, with all of the observed rotation attributable to the disk. 
This disk is assumed to rotate at the circular velocity v c , which we relate to the stellar rotation velocity and velocity 
dispersion, v and a, by v c ~ [v 2 + 2a- 2 ] 1 / 2 . The observed rotation is then v — C"(l — /b) x v c , where the parameter C 
accounts for projection and aperture effects (e.g., a central aperture probably encloses a larger disk fraction than the 
global value). After rearranging terms, we can solve for the bulge fraction based on the observed (v/a) parameter: 



C (v/a) 



2 + C' 2 (v/a)' 



1/2 



(Dl) 



Next, to estimate the numerical factor C , we calibrate aga i nst th e relation between morphology and bulge fraction 
derived from photometric data by ISimien fc de Vaucouleursl ([1986). As shown in the left-hand panel of Figure [Ml 
these authors derived two alternative relations with different extrapolations for the elliptical galaxies, and we select 
the version that extends to B/T = 1, i.e., with some ellipticals being truly disk-less. We then plot the ATLAS 3D 
data set, using their morphological types and (v/a) e values, and estimating B/T using Equation (|D1[) . Based on the 
lenticular galaxy data (Tjjubbic ~ — 2), we judge C ~ 1.35 to provide a reasonably good match; this is similar to our 
adopted model of C mc d = 1-21 for these galaxies (Appendix IA.3[) P^l As in Appendix [Cl there is then an additional 
factor (of 1/2.4 at 1 a c ) to convert our (v/a) observations to the ATLAS 3D values when those are not available. We 
estimate our final B/T values to be generally accurate at the ~ ±0.2 level. 

Next we investig ate the bulge and disk trends for a small sample of elliptical galaxies, making use of the 

pioneering work of iScorza fc Bender] (|1995f) , who carried out decompositions including the kinematics of the bulge 
and disk subcomponents. These authors reported for the subcomponents, and concluded that while the surface 
brightnesses and scale-lengths of these disks were consistent with an extension of scaling relations for lenticular and 



29 IKhochfar et al.l 1120111 ) adopted a similar convention for con- 
necting early- type galaxy rotation with B/T but did not quantify 
a relation except to estimate fy, ~ 0.85-0.9 as the division be- 
tween fast- and slow-rotators, based on semi-analytical modeling 



of galaxy formation. This is consistent with our relation, since 
the slow/fast boundary is typically at (v/a) a ~ 0.15, which our 
relation translate s to ~ 0.85. On th e other hand, the observa- 
tional study of Laurikaincn ct al. (2011) docs not support a strong 
correlation between (v/a) e and B/T. 
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spiral galaxies, the same may not be true for the disk va lues, with those of the ellipticals systematically lower than 
those of the disk galaxies. These data were further used by Ivan den Boschl ([1998) to support a scenario where disks in 
ellipticals are not related to spiral disks, but such conclusions neglected the mass dependence naturally expected for 

Revisiting the lScorza fc Benderl (|1995f) data set, we take their size and rotation decompositions as given, and update 
the distances, stellar masses, and calculations according to our current methodology. From a few galaxies in 
common, we find that these values tend to be higher than what we would derive from scratch (e.g., because of 
decreasing rotation velocity outside the central regions). Nevertheless, we plot the results in Figure [34] (right-hand 
panel), comparing them to the j'^-M* relations that we have found for spiral disks, and for ellipticals overall (which 
we recall is roughly similar to the trend for the massive bulges of spirals). 

There does appear to be a systematic separation between the relations for the bulges and disks of these ellipticals, 
but not as large as we would expect based on the results for spirals. Given the uncertainties in these decompositions, 
we conclude that it is not yet clear whether or not the disk and bulge subcomponents of ellipticals follow the same 
j*-M+ relations as in spirals. If they do, it would be natural to expect the total trend for ellipticals to be offset to 
slightly higher values from the bulge trend (owing to the disk contributions) , which from inspection of Figure [14] 
(right-hand panel) may indeed be true. 



